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THE ANALYTICAL DESIGN OF CONTROL SYSTEMS 


A. M. Letov 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 425-435, April, 1961 

Original article submitted November 4, 1960 


A method of dynamic programming is obtained for the problem of the analytical design of control 
systems that are optimum relative to the integral square error, 


1, The Method of Dynamic Programming 


In [1-13], various methods are described for solving the problem of constructing an optimal system. Solutions 
obtained in these articles give the law of control as a known function of the time and of the initial and final states of 
the system. These solutions can be realized by using digital computers operating according to certain algorithms ob- 
tained. Discrete time intervals, during which the switching of the controlling unit occurs, are used for the output. In 
[1-7], the so-called “synthesis problem” is also formulated, {.e., the problem of presenting the control law as a known 
function of the coordinates of the system. In [14, 15], the authors discuss the simplest case of the solution of this prob- 
lem, which is given the special name “analytical design." 


The problem of the analytical design of control systems was formulated as a problem of classical calculus of 
variations, in which the law of automatic control is obtained relative to some previously stipulated optimizing func- 
tional, The law of control was obtained analytically and written in the form commonly used in actual applications 
with the corresponding sensing elements and servomotor. 


The disadvantage of the methods applying the calculus of variations lies in the fact that they involve writing 
the equations of the variational problem, and these equations must be actually solved, 


In the present article, we use the fact that the dynamic- programming method developed by R. Bellman in [8] 
permits us to obtain results already known by a new method that is sometimes more conveniently applicable in numer- 
ical work, and also to investigate some other possible methods, For clarity, we will describe the basis of the method 
of dynamic programming by using a simple example. (The original exposition of the method is given in [8]). We 
will assume that the differential equation 


t= y) (1.1) 


is given, with the initial condition x(t) = u, ts = 0. 


We assume that G(x, y) is a bounded, continuous function of x and , differentiable with respect to y, defined 
in the open region N(x, y). Here y is a control function that minimizes functional 


= F (x, y) dt (1.2) 


in the class C of functions of x and y, Let y be a function minimizing (1.2), This minimum ts some function ¥ (u) 
of the initial state u which we denote by 


= minJ (y). (1.8) 


Let S be a positive number. We have 


(1.4) 


The essence of R, Bellman's optimization principle lies in the fact that if y = y(t) minimizes the functional 


foe) 
i F (x, y) dt, then it also minimizes the functional F (x, y) dt independently of the first integral. The prin- 
ciple expresses the natural condition of calculus of variations according to which any part of an extremal y(t) is also 
extremal, Using this fact as a basis, R. Bellman developed the following original method. The equality 1.4 can be 
written in the form 
8 
= min [| F(x, 
0 


the function y of the class C;. We will as- 


sume that the number S is sufficiently small and that y(tg) =v. If # has a derivative for xé[x(0), x(S)], then accord- 
ing to the mean value theorem we have 


plu + G5} = 
where 0<€<S. Then 
¥(u) = min| F (u, ») 5+ p(w) + 0(8)]. 
Here 0(S) is the remainder, and it is assumed that lim 0(S)/S~* 0 for S~* 0, Therefore when S ~ 0 we obtain 


min [F (u, v) +G(u, »v) (1.6) 


In order that the expression in brackets actually yield a minimum relative to y, its derivative with respect to 
¥ must be zero, Therefore zr 


F (u, v) + =0, 4 0. (1.7) 
After eliminating 5* , we arrive at the equation 
oF 
(1.8) 


relating the values of u and v for any instant of time ty = 0, The principal aim of the dynamic- programming method 
is to obtain the functional equations (1,7). We will apply this method in solving the problem of analytical design. 
However,we first make two observations, 


Observation 1, Since the equation (1,8) establishes the relation v = v(u) for any instant ts of time, and since as 
a result of the optimal principle it does not depend on this time, the relation v = v(u) is equivalent to y = y(x) with 
the equation (1,1), we obtain the complete system. For simplicity, in what follows we will therefore write the func- 
tional equation immediately in terms of the original coordinates, since this will not cause any confusion. 
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Observation 2, In the case when the function G is defined in a closed region N(x, y) and the optimal solution 
can be obtained from the boundary condition, the requirement that the function y(t) be continuous and differentiable 
can in general be relaxed. The difficulty in solving a problem by the dynamic-programming method lies in the fact 
that the second of the expressions (1.7) can, in general, not be obtained, since it is uncertain what the expressions 
> and 5° signify along the boundary of N. The originator of the dynamic- programming method proposes at this 
stage a very complex method for constructing the optimal solution, We have not applied this method. In the course 
of the solution of the analytical-design problem in the present work it will be shown that a simple optimal solution 
can be obtained, in the case of a closed region N, of the problem for continuous functions y(t), and in a subsequent 
article for continuous x(t) and discontinuous functions y(t), The value of the method we will describe lies in the fact 
that it makes it possible to obtain a closed solution for the problem of analytical design in all these cases, 


2. A First Order Object of Control 
We will start by considering the simplest problem of optimizing the functional (1,2) in the class of functions C,, 
We have the equation of the object of control 


= by + me, =n, (2.1) 


and the optimizing functional 


= \ (an? at. (2.2) 
0 


We write the dynamic programming equations (1,7) in the form 


ap 
an? + + (bn + mt) = 0, 
(2.3) 
+ mz = 0. 
We thus obtain 
mck? + 2benE — amr? = 0 (2.4) 
with the solution 


In the solution of (2,4) for — , we reject the solution 


which does not satisfy the requirements of control stability, The solution (2.5) was obtained in [14], For what fol- 


lows, it is important to obtain the form of the functional #(u). After eliminating & from the equations (2,3), we ob- 
tain 


ay? + = 7 (2.6) 
This nonlinear differential equation is satisfied by 


p(n) = BE), (2.7) 


i- 

n 
6) 
1) 
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where B(E) is an arbitrary function of & ; we should probably set B(E ) = const, 
3. The Optimal Solution in a Closed Region 


We turn again to the equations (2.1), (2.2), and we will investigate their solutions in the closed region N(t,) 
characterized by the relation 


(3.1) 


where € is a given positive number. We will assume that the comparison functions are of the class C, since the opti- 
mal solution of the problem can be obtained from the boundary of the region (3,1), The case when x(t) remains con- 
tinuous while y(t) is discontinuous at a finite number of points will be considered in a subsequent article. 


We therefore consider the single-valued nonlinear transformation [15] 


= (3,2) 
The function #(¢) js continuous, piecewise differentiable, and defined as °* 
for 
—§ for 


The number ¢° is determined by the form of the function ¢(¢). The transformation (3,2) maps the closed re- 
gion N(E ,”) into an open region N(¢,”) in which the role of the controlling function passes to ¢. 


The functional equations (2,3) in the present case will be 


an? + og? + (bn + mp) = 0, 


\ 29 (3.4) 
The last of the equations (3,4) breaks up into two equations: 
_9 
(3.6) 


It is evident that when c #0, the equation (3.5) taken together with the first equation of (3,4) yields the known 
solution (2.5). In view of (3,2) and (3,3), this solution exists for 


k+b 
 al<e’. (3.7) 
From Eq. (3.6) we also obtain 


(3.8) 


Then the first of the equations (3,4) becomes an ordinary differential equation for the functional ¥ with the so- 
lution 


= — — | — (by + mé)| + D. (3.9) 
*No specific form of the function #(¢) for |¢| < ¢* is given, since it isnot necessary in the present considerations. 
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For the purpose of analytical design, the expression for the functional (3,9) will be needed only to check on its 
continuity everywhere in N. Before doing this, however, we write the control equation 


We now assume that the initial deviation Np is sufficiently great in absolute value and that the system starts at 
the boundary of the region N. We further assume that the system ts stable [15]. The expression for the function in 
this case is obtained from the formula (3.9), where D is an arbitrary constant. Then the functional ¥ is a bounded, 
continuous function of 1, since the expression bn = mE is nowhere zero on the boundary. For |n|=|m/(k + b)| F, 
the function ¥ (n) of (2,7) must take the same value as the function # (1) of (3.9). This can always be arranged by 
the proper choice of the constants D and B. 


The solution (3,10), although it is the same as the solution obtained in [15] by using the calculus of variations, 
is interesting because it was obtained by applying a nonlinear transformation of the form (3.3) which avoided the es- 
sential difficulties involved in investigating the functional equations, connected with the impossibility of differenti- 
ating the expression (1.6) with respect to v = € for & =F, 


4. An gth Order Object of Control 


Case of an open region, We consider the general case in which the disturbed motion of the object of control is 
given by the equations 


(4.1) 
Tx (0) = No 
and it is desired to find the function € minimizing the functional 
oo 
I(t) = Vat, (4.2) 
where 
k 
is a positive definite quadratic form, 
The functional is to be minimized in the class of functions Cy. 
The basic functional equation of Bellman has the form 
8 
(tho) = min Vat + pine (4.4) 
0 
When we take the limit when S~ 0, we obtain 
O=V DeaNa 
+ + mat) (4.5) 
Since the function — must make the expression on the right-hand side of (4,5) a minimum, 
m = 0. (4.6) 
+ 2 
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The function ¥ must therefore satisfy the partial differential equation (4.5), where the integral # must be cho- 
sen so that it satisfies the relation (4,6). 


In order to find such an integral, we simplify Eq. (4.5) by using the equality (4.6), As a result we obtain 


4 2 
The solution of this nonlinear partial differential’ equation must be sought in the form of a quadratic form in 
the variables ,: 
= 
ka (4.8) 
Here Ayq ate completely definite numbers which can be found by comparing coefficients after the substitution 


of the expression (4.8) in the equation (4.7). The substitution of the expression (4,8) in (4.6) gives the equation of a 
perfect control with a fixed negative feedback [14]. 


Pete: (4.9) 
a 
We will show that (4.9) is the same equation for a perfect control (3.7) obtained in [14]. With this aim in view, 


we consider the equations (4,5) and (4.6). We will denote the right hand side of Eq. (4.5) by U. Taking (4.9) into ac- 
count, we must assum that 


U =U (m,..-, tn) = 0. (4.10) 


dU = 


Since the equations (4,1) and (4.9) determine a system of independent solutions 1, then (4,11) will be satisfied 
only when the equations 


(4.11) 


= (k =4,..., (4,12) 


and write the equations (4,12) as 


— — (> beat + mit) = 0. 


But from (4,1), 


(2 + mit) = in, 


Hence 
are satisfied. 
We now use the notation 
368 


and so along with (4.6) we obtain 


he = — + 


4.14 
k 


It follows from what has been said that the systern of equations (4,14) is equivalent to the system (4.5), (4.6). 
But the equations (4.14) together with (4,1) are the equations of the variational problem of Lagrange, solved in sec- 
tions 2 and 3 of the work [14]. The conclusions drawn concerning Eq. (4.9) follow from this. A second, no less im- 
portant conclusion can now be drawn. Namely, as was shown in [14], the equations (4,1), (4.9) describe a stable con- 
trol system having a control function that optimizes the functional (4,2) in all nonspecial cases as noted in [14).° 


Case of a Closed Region, We again turn to the equations (4,1) and (4,2), and assume that the coordinate — is 
restricted by the inequality (3.1) and that the functional (4,2) is minimized in the class of functions C. If we repeat 
the reasoning in part 3, we obtain the basic functional equation in the form 


V+ + mee) = 0; 


As before, we obtain two solutions; the first coincides with the solution (4.9) for the points | £ pat. | <E, the 
second solution @ = + € exists for points | Zp.na|>§. In order to check that in both cases the first equation of 
(4,15) defines a continuous functional, it is necessary, in addition to the solution (4,8), to find the solution of the par- 
tial differential equation 


(4.16) 


+ + (Toran + mb) = 0. 


We will, however, not consider this further, 


It is obvious, from comparison, that the dynamic-programming method yields the same result as was obtained 
in [14, 15] by the methods of calculus of variations with the application of the functional (4,2), This makes it possi- 
ble to obtain directly the law of control in the form in which it is usually realized with a known selection of sensing 
elements, transformations, amplifiers, and servomotors, 


5. The Case of a Limited Servomotor Velocity 


In this case (see [15]), addition to the equations (4,1), we must use the equation 


E=/(5), (5.1) 


where is a function of the class A [16]. The function f (9) and its argument . . . %p,€) must be chosen to 
minimize the functional 


in the class C of functions of 1), Mp, and €. 


*In this connection, see [5, 18]. 
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Reasoning similar to that used above leads to the Bellman equations 


Vt + 1% + + mb) 55, = 0 


+) = 0. 


The desired solution consists of two parts, The first part is determined by the limitations imposed on the veloc- 
ity, and is obtained from the condition 9f/do = 0, This means that 


(5.3) 


= +f, (5.4) 
where 0° is a positive number. 
For obtaining the second part of the solution, we have 


In combination with the first equation of (5.3), this yields the partial differential equation 


++ mab) aa (%) 
which can be satisfied by the function 
= +8 Bette + RE. (5.7) 


_ The coefficients of the form (5.7) are obtained by the method already described, The equality (5.5) yields a 
linear function f(9) for |o|< 0°, i.e., 


The quantity 0° is obtained from the condition ho* = f, 


As in the previous case, we can show that this solution of the problem coincides with the solution obtained in 
(15). 
6. An Example 

As an example, we will consider the problem of an aircraft flying at a specified height. With certain assump- 


tions concerning the rate of damping of natural oscillations of the plane relative to its transverse axis, the equations 
for the disturbed motion can be written as 


= bum + + maf, 
= baits + + mak, (6.1) 
Ns = 


The coefficients of the equations are given in terms of the initial parameters of the aircraft by the relations [17] 
bu = —c,,, bia = my == — per, , 


(6.2) 
bu = bes = — 2kyc,,, ms = pet, 


Here t is a parameter of the aircraft. In the equations (6.1), 1, is the inclination of the trajectory, 1, a dimen- 
sionless velocity disturbance, and Ns a dimensionless disturbance of the height, 


The optimizing functional is given by the integral 


oo 
1) =| (ng + (6.3) 
0 


The basic equation for the determination of the functional is 


and the solution of this equation will be sought in the form 


@ = 4, 2, 3). 
ka@ 


The equations for determining the coefficients are 


budnt + Ara = + 


= (myAj3 + 


(6.8) 
+ (ben + bus) Asa + + Aas + mE mms (An Aart 


+ + Ase + m3 + mamta + 


The equations (6.6) can be solved by the method of successive approximations, and as a result we obtain 


i 


1. A Relation between Dynamic Programming and Lyapunov's Functions 


In a discussion of the contents of the articles [14,15], N. N. Krasovskii drew my attention to the fact that the 
functional (4.2), interpreted in these articles as an integral square error, can be considered as a Lyapunov function, 
N. N. Krasovskii pointed out this relation between Lyapunov functions and optimal systems in [7]: and also in [18, 19}. 
R, E. Kalman noted it later in [13, 20}, 


We will show here that the relation of Lyapunov functions to the problem of optimizing the integral square error 
arises naturally from the formulas in section 4, Actually, using (4.9), we will calculate the total derivative of the 
function defined in (4.8). This function must be positive everywhere that In; + 0. It will have a meaning 
when the system (4,1) and (4,9) is stable. 


Then the basic functional equation (4,5) becomes the partial differential equation we have used for construct 
ing a Lyapunov function with a total derivative given by the formula 


dtp 
(1.2) 


so that the functional ¥ is the Lyapunov function for our problem, Conversely, we assume that by applying dynamic- 
programming methods we have found tha functional (4.8). If there is some way of verifying that the function # is 
positive everywhere for =nj = 0, then it is a Lyapunov function for the problem. The condition that # be positive 


is the stability criterion for the optimal system, This is of special interest, since the dynamic-programming method 
can be used to solve simultaneously both the problem of optimizing a system according to the criterion (4.2) and the 
problem of the stability of the optimal system, An optimal and stable system corresponds to a functional ¥ that is 
positive for Int + 0. 


This latter fact is especially important in the case when the system to be optimized is nonlinear 


= Ex Na» §). (7,2) 
In this case the basic functional equation will be 


k 


and the relation (7,1) remains valid, 


The author wishes to thank E, A, Barbashin and N, N, Krasovskii sincerely for their very useful advice and com- 
ments, 
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ON NEW STABILITY INEQUALITIES 


I, S, Arzhanykh 


(Tashkent) 
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Original article submitted October 25, 1960 


New inequalities based on Shur's theorem are given; these are useful in the analysis of the stability of 
systems described by differential equations with constant or periodic coefficients, 


Suppose the variational equations 


Pots (s == 1, 2, 8) (1) 


ive constant or periodically varying coefficients, Then, as is known, the question of asymptotic stability may be 
solved by inequalities guaranteeing that the roots of the characteristic equation are distributed in the left half-plane 
or in the unit circle,respectively, With the help of a linear fractional transformation with real coefficientsone may 
transform the left half-plane onto the unit circle and thus consider in both cases the question of whether the roots lie 
inside this circle, The criterion that the roots of the characteristic equation lie in the left half-plane is given by the 
Raus-Hurwitz inequalities, However this criterion is very inconvenient in application, since it is first necessary to 
write the equation in clear polynomial form, and then, from the coefficients of this polynomial, form the Hurwitz de- 
terminant, The question naturally arises: Is it possible to obtain an equivalent inequality without first reducing the 
equations to polynomial form? 


A positive answer to this question is given in the present note, We shall seek inequalities guaranteeing that the 
roots be distributed inside the unit circle, and at the end we shall indicate the slight changes necessary in order to ob- 
tain inequalities guaranteeing that they lie in the left half-plane, 


Thus suppose the characteristic equation has the form 


According to Shur's theorem [1] the necessary and sufficient conditions for the roots of equation (2) to lie inside 
the unit circle consist in the satisfaction of the inequality 


| ¢o| >| en| (3) 
and the condition that the polynomial 


hi (2) +...+c 


n—1’ 


(4) 


defined by the equation 
xfs (2) = Cof (x) — (2), (5) 
where 
f (2) = 24 (6) 
possess the same property. 


~ Using these algorithmic rules, it is easy to construct stability inequalities without writing out f(x) explicitly, 
We have 


f (2) = [A + Bz), (7) 
Qin 


Onn 
Cn = [Ax + = (b"). 


Cy = [A + = = (8) 


(9) 
We calculate f,(x) by the formula 


af, (2) = Az + B]—cnlA + Ba] 
and {* (x) by the formula 
file) = 4+ 8+))} = 
= ¢)[A + Bx] —ca[Az + B]. 


fala) Aa + (A + Ba}) = ... + ,, 
fil) +... + 


= Lim fi (2) = [A] — B] = (6") = (10) 
d 
ef), = lim f, (2) = + B) — cal + 
Denoting 
[Az + B] = [ax + 6, ax+b,..., az + 6}, 


(i1) 
[A + Bx] = [a+ ba, a+ ba,..., a+ bz), 


—en({b, @)+[a, +[a, a,..-, N= (12) 
= — cn (ba). 
After this we calculate 


= [42 + B)— ent A + Bz}) — (¢,[A + Bz] — [Az + By} 


a*f, (x) = [Ax + B] p; (x) —[A + Bz] qi (2). (13) 


Ps (2) = py + chs (14) 


qi = + ap 


3) 


4) 


Po = = (b”). (15) 
We show that the formula 


(xz) = [Ax + B] pe_s(z) — [A + Bz] (#), (16) 
holds, where 


pr(x) = c™ pps (a) + 29,_, 
gr (2) = (2) +, (2) 


(17) 
for all values of k= 1,2,...,-—1, With this goal we compute 
= (A+ Bz} — [Az + B)q,_, (+) 
and form the expression 
af k++ (x) = f k (z) — k = 
{{Az + B) p,_, (@) —(4 + Bal 
— [A + Bal + BI 
From this we find 
(2) = (Az + Bl (2) + (>) 2)— 
= [Az + p, + Bz] q, (2), 
where the polynomials py and q, are determined by formulas (17). For k = 1 the formulas (16) and (17) have the 
form (13), (14), and (15), We extend by induction formulas (16) and (17) to all values of k= 1,2,...,0—1, 
Suppose 
fy (2) = ++ + (19) 
Then we have 
a® 
ct), = lim f(x) = + Bl (2) + Bel (20) 
Moreover, 
fe, (x) = ++ +++ (21) 
(x) = (A + Br] (x) + (2)- (22) 
Consequently, 
= lim f; (2) = { ([A + Bz] — [Ax + Bl 


Thus we obtain the following theorem. 


Theorem. Necessary and sufficient conditions that the roots of equation (2) lie inside the unit circle are that 
the following inequalities be satisfied: 


> lenl, 
| | >| | (k=1,2,..., (25) 
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where the numbers ef” ana f, are defined by formulas (20) and (23), and the polynomials p; and q; are calcu- 
lated recursively according to formulas (15) and (17). 


We examine the application of the theorem to the cases n = 2 andn = 3, For n= 2 we have: 


[Az + B] = [ax+ 6, ax+b], [A+ Br] = [a+ bz, a+ bz}. 
= (a, a] = (a*), c,=[b, = (0°), 


(cy [Ax + + Bz\)| 
= ({ab} + (ba]}) — ({ba} + [ab}) = (¢, — (ab). 
The stability condition is 
lel > lel |. 
By way of a numerical example we consider the equation 


f(x) B] = 0, 
ay, = 1, by, 
aq, = — 4, ag, = 1, ba, = 1, bes = B. 


For the roots to lie inside the unit circle it is necessary and sufficient that the inequalities 


Jeo” 
= c8— = (a? (a? +.B + 2), 
be satisfied. 
~ Consequently, one should have 


For x = 3 we have 


c() = (ci)? = (elt) — (¢, (a*b) —c, (ab*)). 
The stability inequalities have the form: 
If the variational equations (1) have constant coefficients, then it is necessary to transform the ) -plane into the 


x~plane so that the left half-plane is mapped onto the unit circle of the x-plane. This is accomplished by the substi- 
tution x = (1 + A)/(1— A), as a result of which the equation 


goes into equation (2), where 


A = (p,, —5,,); B= (p,,+6,,) 
(Srs is the Kronecker symbol). 


Pa Pog — Pan a: 0 

Pug Peg 
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In those cases then the system of variational equations consists of equations of order higher than the first, by 
the introduction of new unknowns it may be brought into the form of a system in which each equation is of order one; 
then the inequalities indicated may be applied. But such a case may be investigated directly by computing the Shur 
polynomial in the above manner. We examine, for example, the variational equations of the form 

+ disks + disks + Pinks + + Crake = 0, 
+ + + + + Carbs = 0 


with constant coefficients. The characteristic equation will be 


+ +- Cor, + + Cos 


We transform the left half-plane (Re A < 0) onto the unit circle (|x| < 1) by the substitution 


=x Q, 


We examine the characteristic polynomial 


where 
=A, +b,+¢,, By, Vin = — + Cy. 


in the unit circle, To apply Shur's algorithmic rule it is neces- 
sary to compute the coefficients cy, cy, We have: 


Then we compute the third degree polynomial A, (x) by the formula 
(2) = (2) —c4A° (2) = + 28 + 1] lo + Ba + 
and find 


= limA,(z), ef) = lim (2) | 
x0 


= + 8+ t+ Be+ 
= (Bt) — (Ba) = + — + = 


| 
= + + — + eB + = = — 


The second degree polynomial A,(x) is calculated by the formula 
A, (x) = cA, (2) — = 


= (ey [22 + 28 + 11— + Be + — 
= {fata + + (c,c + — [ar + Bar + + 


= Vim (2) + + (2) — + + 9,(2))}_ = 
= + (BY) + (Br) — (rat) + (BD) — (Ba) 
= lim lim {{a + fie + + — 


— fae? + Be +11 + 
Finally, we compute the first degree polynomial A(x) by the formula 
x, (2) = oA, (2) — (2) = + 28 + + — 
— + Ba + (cc + — 
— fs + B+ + -- 
—[2+B- + + = 
= (lata + 2B + Ps (z) — + Ba + (z)}. 


= 2 (c, (aB) — (1B) + (ce, ((av) + (B*)) + 
+ ((1*) — 
P= lim lim {fa + Be + P,(-)— [ae* + Be + 7] q,(+) 
The stability inequalities have the form: 
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The paper studies a method for correcting automatic control systems which is based on changing the 
structure of the system as a function of the signs for the deviation of the controlled parameter before 
and after those control system elements whose properties are to be corrected. The method is analyz- 
ed using the example of a relay automatic control system developed by the author for controlling the 
electric power supplied to an ore furnace. A comparative theoretical and experimental investigation of 
the system is performed for the cases of correction and no correction. 


1, Statement of the Problem, A reduction of the effect of noise in automatic control systems (ACS) in those 
cases where the useful signal has a spectrum that consists of lower frequencies than the noise is achieved by introduc- 
ing a series filter [1]. However, when this is done the narrowing of the frequency band causes a reduction in the sta- 
bility margin of the system, and the duration of the transient response increases. Varying the gain in such a system 

can often lead to a worsening of the static or dynamic properties. In relay ACS the introduction of such a filter re- 

duces the number of switching operations [2] and therefore increases the operational reliability, but (as in other sys- 
tems) it worsens their static or dynamic properties. This contradiction could be eliminated if we were able to design 


an ACS with a narrow frequency passband relative to the perturbation and a wide passband relative to the useful sig- 
nal, 


In this paper we show that this problem can be solved by modifying the structure of a relay ACS . For this pur- 
pose it is necessary to cut in a filter, depending on the signs of the quantities at the input and outputof the filter. This 
can be achieved by means of a special corrective network. Such a system will have good dynamic and static proper- 
ties for a minimum number of switching operations for the relay elements. 


2. The Essence of the Correction Method. The idea for the correction method treated in this paper consists of 
the following. Assume that after introduction of the filter the transient response at the input of the ACS is character- 
ized by decaying oscillations 1 (Fig. 1). The variation of the quantity at the output of the filter which is connected 
in series with the detector element has the response given by curve 2 under these conditions. This same ACS with no 
filter does not cut off the high-frequency component of the perturbation. Its transient response is characterized by the 
curve 1". The problem of making the transient response 1 approximate 1° while simultaneously filtering the perturba - 
tions can be solved as follows. 


We shall introduce a network which compares the signs of the quantities at the input and output of the filter but 
worsens the dynamic properties of the ACS. Assume that at instants when the sign of the ordinate of curve 1 changes 
relative to the sign of the ordinate for curve 2,the structure of the original system will cahnge in such a way that dur- 
ing these periods the filter will be disconnected. 


The nature of the transient response in such a system with a variable structure can be represented by the curves 
3 and 4, respectively. 


A sign-sensitive element is connected to point a at the instant the sign of curve 1 changes relative to curve 2; 
this element disconnects the filter. The curve ab is then characterized by a steeper slope than that of the correspond - 
ing curve bc. Then the sign sensitive element which changes the structure is again disconnected, and this restores the 
original structure for a smaller error magnitude. 


It is evident from Fig. 1 that the transient response in the system with correction is completed during @ shorter 
time and approaches the transient response in the systern with no filter. Thus, any deviation which has appeared ini- 
tially and has been caused by the perturbation passes through the network elements which assure the appropriate filt- 
ering of the high frequencies; the signal passing through the feedback loop is received in such a way that it bypasses 
the filtering unit. 
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The general statment of the problem involved in 
changing the system structure as a function of its operating 
conditions was made by V. V. Solodovnikov [3], and for 
pulse and relay systems it was formulated by Ya. Z. Tsyp- 
kin (4, 5}. The literature indicates the expediency of us- 
ing quantized devices for this purpose [5]. In [6] a study 
is made of the introduction of nonlinear feedback into a 
linear ACS for the purpose of achieving the optimal tran- 
sient response. The paper describes the principle of step- 
wise variation of the speed feedback of a servosystem as a 
function of the product obtained by multiplying the sumof 
the error and its derivative by the error derivative. In [7] 
the possibility of controlling a neutral object with an asta- 
tic controller via the introduction of a quantized element 
is demonstrated. 


In [8] the author proposes the introduction of a non- 
linear element which varies the damping coefficient as a 
Fig. 1. 1) Transient response at the input to function of the error magnitude into a linear servosystem. 
the filter in the ACS with no correction; 

1’) transient response in the ACS with no fil- 
ter; 2) transient response at the output of 
the filter in the ACS with no correction; 3) 
transient response at the input to the filter 
in the ACS with correction; 4) transient re- 
sponse at the output of the filter in the ACS 
with correction; 5) connection of the sign- 
sensitive element 24 @p is the dead band 
for the sign-sensitive element; § is the trans- 
ient response time in the ACS with no filter; 
Tis the transient response time in the ACS 
with a filter when correction is present; T,, 
Tz, and T; represent the periodicity for con- 
nection of the sign-sensitive element. 


The block diagram for designing the ACS is shown 
in Fig. 2. Here 1 is the controlled object, 2 is the detect- 

. ing element, 2a is the comparator element, 3 is the series 
corrective network (a filter that filters out the higher har- 
monics of the perturbations, 4 is the control network (for 
example, a relay system), 3' and 4‘, respectively, represent 
the corrective and control networks in the auxiliary control 
loop, 5 is the actuating organ, 5° is the corresponding ac- 
tuating organ whose presence is not always mandatory when 
tme element 5 cau be connected to tne circuit 6-3'-4-1, 

6 is the network which responds to the sign of the controlled 
parameter deviation relative to its preceding value (Xs). 
In the relay system under study (Fig. 2) the total transient 
response time will be a minimum (other conditions being 
equal) if the pulse in the auxiliary circuit causes the ac- 
tuating organ to be disconnected. 


In a number of cases it is expedient to use quantities differing from the above-mentioned ones for correction 
purposes. In the general case we assume that the filter consists of a set of certain elements and sections in the ACS 
which has definite dynamic properties that are to be ccrrected. 


In order to exclude the effect of noise on the corrective network it is important to match the static and dynam- 
ic characteristics of the sign-sensitive element 6 to the amplitude-frequency or spectral responses for the noise; these 
are not difficult to obtain on the basis of investigating the object and the operating conditions of the ACS*. The 
majority of industrial ACS are characterized by the effect of noise signals which have a relatively small amplitude 
compared to the amplitude of the basic perturbations, and a wider frequency spectrum. In order to suppress noise of 
this type we use a filter in the basic circuit. In that case the effect ot noise on the corrective network can be prac- 
tically eliminated if the sign-sensitive element 6 also responds to the absolute magnitude of the quantity xs and per- 
forms the switching when xs and xs are commensurate. 


In the general case,eliminating the effects of noise on the corrective network is possible if the static and dy- 
namic characteristics of the elements in this network are determined as functions of the amplitude-frequency respon- 
ses for the noise. When this condition is not satisfied we will have a situation where sectors are cut out of the noise 
spectrum cutoff by the basic circuit; these sectors correspond to the passband for the corrective network frequencies 
at instants when the sign of the noise xg is opposite to the sign of xs after the transient response has arisen. This pro- 
duces an appreciable lowering of the degree to which noise affects the ACS. 


*~Such characteristics are treated below on the basis of a specific example. 


The method treated above for changing the structure of the system in order to obtain specific dynamic or static 
characteristics relative to the useful signal and the pertrubations (as a function of the signs of thesignals at the inputs 
and outputs of the elements whose properties are to be corrected) can be used for improving the dynamic properties of 
a broad class of industrial ACS in which such a problem arises. 


We shall study the application of this correction method to a specific relay system. 


Fig. 2 


3. Description of the Network. We shall analyze the application of 
the method above to ACS which control the electric power delivered to 
an ore furance [10,11], The control network (Fig. 3)° is designed on 
the basis of measuring and controlling the phase admittance of the fur- 
nace by means of an electromechanical device for obtaining the quo- 
tient resulting from a division of the phase current by the correspond- 
ing voltage. In an electric ore furnace the variations in the phase ad- 
mittances of the furnace are independent [12] when the proper distances 
between electrodes are chosen. 


If under these conditions the dynamic characteristics of the elec- 
tromechanical computing device are matched with the corresponding 
characteristics of the network, then the autonomous operation of the 
controllers in individual phases is achieved in both the static and dy- 
namic modes. In order to ensure the reliability of the network and cor- 
respondence between the conditions governing the operation of the ore 


furnace it is important to obtain cutoff of the high-frequency noise components and a minimum displacement of the 
electrodes. This is achieved by choosing the frequency responses of the electromechanical device that measures the 
ratio and by assuring autonomous operation of the controllers in neighboring phases. 


| 
! 
! 
! 
! 
! 


7 


L 


Fig. 3. TrF is the furnace transformer; E are the electrodes; F 

is the fumace. 1) Detecting block; 2) compensation servomotor 
(CM) with control contacts (Ryy); 3) auxiliary relay units; 4) con- 
trol network. EM is the motor for the electrode drive; H, L are 
the starting coils for the electrode-drive motor; S. N, is the set- 
ting network; FSH and FSL are the final cutouts for the electrode 


drive. 


: The network is shown for one phase of the furnace; the networks for.the remaining two phases are analogous. 
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Figure 4 shows the amplitude-frequency responses for the noise and fundamental perturbations, obtained on the 
basis of investigating the ore furnace as an object for automatic control of electric power. The curves are plotted in 
percent of nominal value of furnace phase power for a lead-smelting electric furnace; the data for the furnace is given 
in the Appendix. From Fig. 4 it is evident that the noise has a spectrum with higher frequencies and an amplitude- 
frequency response that lies appreciably below that of the fundamental perturbations; this verifies the reasoning cited 
above and demonstrates the expediency of introducing a filter for the purpose of excluding the effect ofnoise on the 
system and eliminating excess oscillations of the controlled parameter at a frequency close to the noise frequency. 
From these curves it is also evident that the superposition of noise on the fundamental curve has practically no effect 
on the average value of the controlled parameter. 


Our network is sufficiently characteristic, since it 
takes into account the special operating features of a ma- 
jority of industrial control systems that include electronic 
automatic potentiometers, bridges and other sections with 
compensation detecting or control instruments ("EPD,” 
“EMD,” "IR-130," "EPP," etc.). Just as in many industrial 
systems, we are dealing with a relay system that incorpor- 
ates a constant-speed servomotor. 


Ale) The network operates as follows. When there is no 
deviation of the controlled quantity the contacts Rjy: and 

: Rye (Fig. 3) of the electromechanical device are open, 

5 Ay(w) <+-- N = the armature of the null-relay R which controls the con- 
TT | nection of the compensation motor CM is in the interme- 
0 Qi 05  1.5w.sec™ diate position, and both contacts R;, Rp are open. The re- 
Fig. 4. A» (w) is the characteristi for the fun- lays 1SR, 2SR, which operate when the contacts of R areen- 
damental perturbations; A; (w) is the character- ergized are engaged, and the normally open contacts 1SR 
istic for the noise during the smelting period; and 2SR in the circuits of the starting coils for the electrode 
Ay (w) is the characteristic for the noise during drive motors are closed. The relays 1SR and 2SR together 
the smelting period —the period ir which with the control contacts Ry: and Ry form an element 
physicochemical ‘transformations occur in the fur- which is sensitive to the sign of the variation in the control 
‘A Calf for parameter relative to its previous value. When the admit- 
the perturbations; A,(w) is the characteristic of tance of the furnace changes, one of the contacts of the re- 


the basic ACS circuit. , 


tion of deviation for the control parameter corresponds to 
closure of Ry and Ry; the other direction corresponds to closure of Rz and Ry. The pulse for connection of the ele- 
trode drive motor is transmitted through the electromechanical device from the contacts Ryy; the cutout pulse is sent 
from the contact 2SR and SR. The actuating motor which displaces the electrode is disconnected before the instant 
at which it would be disconnected due to opening of the control contacts R,;. This makes it possible to eliminate the 
effect of the inertia and the time required for movement of the compensating servomotor in the detecting instrument. 


It should be noted that in the nctwork being analyzed the limited frequency range of the instrument used for 
frequency measurement is directly associated with the principle governing the operation of the electromechanical 

device; the output quantity of the device reaches a steady state only after the compensation process has been com- 
pleted. 


Analysis of the control System. Based on the analysis of the automatic control network, its block diagram( which 
takes into account the dynamic properties of the individual elements) can be represented as in Fig. 5; in this figure 
1 is the controlled object (the electrical ore furnace); ¢ is the deviation of the controlled quantity in relative units 
¢=G/G,, where G is the running value of the phase admittance for the furnace, Gp is the admittance of the furnace 
when it is used as the base quantity. The input quantity for the object is the displacement of the electrodes p = h,/Neby 
where he and hgp,.respectively,denote the running and base values for the electrode displacement. Studying small de- 
viations p and ¢, we obtain the linearized equation for the object in the form 


TT? 


al 

2s T 


g=ky (1) 


1) 


The transfer coefficient k for the object determines the slope of the static characteristic and depends substan- 
tially on the engineering conditions governing the operation of the furnace (temperature, the state of the melt); this 
coefficient varies as the smelting process proceeds. The detecting organ includes the relay amplifier 2 and the inte- 
grating section 2a (Fig. 5) in accordance with the network shown in Fig. 3. For the condition that the frequency re- 
sponses of the instrument and the network are matched, the output quantity of the instrument reflects the input quan- 
tity with practically no inertia: the presence of the integrating element in this case is equivalent to introducing a 
constant lag corresponding to the maximum transient response time for the instrument. Then the equation for the de- 
tecting instrument is written as: 


o=0 for 
S=NP for Peo (2) 


where gp is the sensitivity threshold of the instrument; %gm is the upper limit for the measurement, and” is a pro- 
portionality coefficient n =do Aig. The lag introduced by the instrument in the other network elements is taken into 
account by section 4. 


In a case where section 2a has the characteristic of an integrator® , ¢' (Fig. 5) corresponds to a pulse represent- 
ing the derivative of 0. 


The function for the transducing element 3 (Fig. 5) is ful- 


f f filled by the main contact of the power controller. The equa- 
L. a ial tion for element 3 corresponds to the equation for a relay witha 
as 2 | | 22 dead band (240) which is determined by the controller adjust- 

ment 
0=0 for (3) 
for |5| >| As}. 
4, 


The equation for the delay section 4 in the basic circuit 
1-2-2a-3-4-5 is written as 


Fig. 5. e=/(0—0,) for |8|>0, (4) 


where 6, = 6(t°— 1), where ¢ is the time corresponding to the instant at which the quantity % changes (d9io » 0). 


The servomotor equations which take into account the difference between the dynamic torques of the drive for 
raising and lowering the electrodes and the dead band which accompanies connection and cutout are written as follows: 


0 e<0 
0<e<e, 
0<e<es p< 0 pO 


0< e< p>o (5) 
<e<0 


& > 

& < Los 


* In the particular case involving the correction under study. 
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Here €o4 and €o3,respectively, determine the dead band for connecting the servomotor (for movement up and 
down); €, and €, analogously determine the dead band for disconnection of the servomotor when it moves up and 
down; Myo; and Myg. are the dynamic torques of the servodrive when it is connected for upward and downward dis- 
Placement; my, and Mye are the analogous dynamic torques for cutout: 


where M, gy, and M) 4, are the average starting and braking torques for the motor which are determined on the ba- 
sis of its mechanical characteristics; My, is the nominal torque; T is the electromechanical time constant 


tm 


and ¥,,, is the maximum rate of change for the controlled parameter (¢) of the object. 


The auxiliary relay unit 3'(Fig. 5) takes into account the sign of the variation in g'relative to 6, Taking this 
into account, we write the equation for the element 3° in the form 


for sign® = signg’, (6) 
for signO--signg. 


The lag introduced by the elements of the corrective circuit 1-2-3'-4'-5 is taken into account (just as in the 
basic circuit) by the delay section 4° whose equation is written as equation (4): 


e’ = (0’ —0,,) for |6’|>0; (1) 


= 6’ (t’* — +), t” is the time corresponding to the instant at which the quantity 9'changes. 


The equations (5) for the servomotor when it operates in the corrective circuit will include ¢€' instead of ¢ 
while keeping the same form. 


The block diagram under study differs substantially on the networks used in industrial systems with quantized 
stopping of the actuating mechanism. In these systems the cutout signal is not associated with a variation in the pro- 
perties of the ACS elements and their static and dynamic characteristics due to the effect of ac perturbations; rather, 
it is determined by a constant standard input introduced from the outside. In that case, for example, the variation of 
the system gain can lead either to the appearance of oscillations or to the prolongation of the transient response (i.e., 
in the final analysis it leads to a deviation of the static or dynamic characteristics of the ACS from the desired ones). 
In the block diagram under study both the basic and the correcting pulses pass through the object; this permits us to 
associate the introduction of the correcting pulse with the variation in the properties of the stabilized object (or has 
its output quantities varied according to some definite law) by the ACS. 


In order to analyze the method of correction under study and to determine its effectiveness we shall perform a 
comparative investigation of two networks: a) with no corrective circuit (the control system 1-2-2a-4-5-1 (Fig. 5)), 
and b) with a corrective circuit (the control system which depends on the sign of g'— 0:1-2-2a-3-4-5-1 or 1-2-3'- 
4'-5-1). We analyze each of these systems on the basis of investigating the corresponding phase trajectories. 


It follows from eqs. (1)—(7) and the block diagram in Fig. 5 that the investigated system can be treated from 
the point of view of dynamics as the combination of a relay servomechanism which is derived by a second order 
equation (in the general case with variable coefficients) and a delay section. Under these conditions the characteris- 
tic of the servomotor has a dead band and is asymmetrically arranged relative to the axis of abscissas; the time de- 
lay may be different, depending on the direction in which the servomotor is made to rotate, An analysis of such sys~ 
tems in the general case requires an examination of the possible movements on a many-sheet phase surface [9]. 


GD 
see 


Fig. 6. Phase trajectories for the ACS: a) with no correction; b, c) for 
introduction of the corrective circuit. The broken lines show the tra- 


jectories for an analogous network without correction. 


Based on a simultaneous analysis of equations (1)—(7), it is not difficult to eliminate the intermediate variables 


and to obtain the equations of motion for the ACS on the phase plane*. For a system with no corrective circuit the 
phase representation is shown in Fig. 6a. It follows from an examination of the phase trajectory that the limiting cy- 


cles will be absent if we have 


for all curves of the form 4° (Fig. 6a); here 24 ¢is the dead band for the controller (4 g=40/n). From this the ne- 
cessary and sufficient condition for an absence of periodic modes and for a transient response without oscillations is 


Written in the form 


* The derivations of the equations of motion for the ACS are given in the Appendix. 


the system may also be stable, but the transient response has one fluctuation. 


The main feature of the networks under study from the point of view of dynamics resides in the fact that in a 
network with correction any connection occurs with a certain lag when the pulse passes through the integrating block 
of the instrument; disconnection occurs in such a way that the integrating block is bypassed. However, in the net- 
work without correction the connect and disconnect pulses pass through the integrating unit in both cases. The filter- 
ing properties of the circuit 1-2-2-3-4-5-1 (Fig. 1) are taken into account by the delay section 4, The time delay 
T, determined by the section 4° characterizes the delay of the relays R, 1SR or 2SR and the servomotor; it is much 
smaller than T. In accordance with this we determine the positions of the switching lines on the phase plane for the 
network with correction. For connection of the system the switching lines for the phase trajectories of both networks 
are analogous, and for disconnection the phase trajectories for networks with correction are shifted in the direction of 
an earlier transfer of the trajectories for the image points onto the plane ®g. The nature of the phase trajectories for 
the analyzed network with correction is shown in Fig. 6b, where the switching lines for the network without correction 
are denoted by broken lines, The graphsin Figs. 6a, 6b,and 6c were plotted for the case where all analogous sections 
of the compared networks have identical parameters. This permits convenient determination of the effect achieved 
in improving the dynamic properties of the system when it is designed according to the described method. 


The equations for the phase trajectories of the network with correction do not differ in form from the original 

equations (cf. Appendix I), The condition for stability and the absence of oscillations is written analogously with (9): 
T 

Taking the relationship T,;« T into account, we come to the conclusion that for equal object and controll- 
er parameters the network with correction permits the limitations imposed on the magnitude of the dead band and 
the characteristics for individual elements to be reduced. In particular, for the same dynamic properties a network 
with correction can be used in conjunction with an object that has an appreciably higher transfer coefficient than the 
value which can be handled by the original network. Here the quality of the transient response in both networks may 
be identical. The relationship which determines the possible increase in the transfer coefficient for comparable net- 
works is not difficult to write on the basis of (9),(11) if we assume ¢ const for both of the compared networks and 
an identical rate of electrode displacement (up and down, respectively). Then after we substitute ¢ pp = kil pp in (9), 
% me = Kulp in(11) and perform the simultaneous solution of (9) and (11), we have 


(12) 


t+ 
where ky and k are the maximum allowable transfer coefficients in the network with and without correction. For the 


condition 7 ,> = we have approximately 


(13) 
ky Thy 


Based on the resulting relationships we determine the boundary value of the transfer coefficient and the possible 
limits for the improvement of the dynamic properties of the network under study. The use of the network with corree- 
tion for the specified object permits us to reduce the dead band of the controller® or to narrow the range over which 
oscillations exist for possible deviations of the system gain (the variation of the transfer coefficients of the object and 
other ACS elements during the process of operation). 


The phase trajectories cited in Fig. 6 illustrate the use of the described correction method in a specific non- 
linear‘system. Here the phase trajectories are plotted for the case where introduction of the corrected circuit permits 


. For an equal degree of stability. 


us to narrow the range in which oscillations exist while si- 
multaneously assurirg cutoff of the high-frequency noise 
component. For specified initial conditions and for main- 
tenance of the dynamic characteristics of the system, the 
introduction of the corrective circuit permits reduction of 
the static system error 24.9(Fig. 6c). Assuming that we 
have the same static error for networks with and without 
correction, it is not difficult to plot the phase trajectories 
which illustrate the expansion of the region for the initial 
conditions (%p, %p) or which assure that the transient re- 
sponse approximates the optimal response. Here the sys- 
tem with correction which is under study can assure the 
minimum displacements of the actuating organ for a mini- 
mum number of relay-element switching operations as com- 
pared to the analogous network without correction. 


We eliminate the effect of noise on the corrective 


circuit in the system under study by satisfying the condi - 
tion (Fig. 4): 


2Aq, > [A, ()) max: (14) 


where 2Aggis the dead band for the auxiliary relay unit in 

the corrective circuit. This does not affect the static char- 

acteristics of the ACS which are determined by the quan- 

| tity 4¢ in the basic circuit. However, an appreciable in- 
crease in Ag» may lead to a reduction in the stability mar- 

Fig. 7. 1) Electrodes; 2) frame of the model; 3) elec- That 

wolyte; 4) metal trough which stimulates the molten 4." ftece of noise will be reduced in the high-frequency 

metal. AM is the actuating mechanism; CT are cur- which is abl : Pg oid fin 

rent transformers; EA is an electronic amplifier; H is in this is 

the electrolyte level; h, is the immersion depth of the 8 


achieved by varying the time constant in the circuit of the 
coils 1SR and 2SR for rel 


Such a compromise solution permits us to assure the 
desired degree of reduction for the effect of noise on the corrective network in both systems. On the other hand, the 


possible effect of noise on the corrective circuit may manifest itself in terms of an earlier cutout of the system; from 
the point of view of reducing the possibility of the appearance of oscillations this is useful, although it may lead to a 
certain extension of the transient response without causing any deterioration in system stability. 


For the analyzed system the effect of noise on the corrected circuit is not substantial*®, The introduction of a 
corrective network, even when the analyzed effect of noise on the corrective network proper is taken into account, 
solves the problem on an over-all increase in the stability region for the ACS. 


The analysis above illustrates a particular case where the correction method is applied. In general form the in- 
vestigation and the choice of parameters for ACS with correction will be studied separately. 


5. The Experimental Characteristics of the Network. An experimental investigation of the compared networks 
was performed in an electrolytic bath (an electrical model of the electrical ore furnace*) by connecting both con- 
tollers to it. The scheme illustrating the tests and the method of obtaining oscillograms is shown in Fig. 7. For in- 
stantaneous variation of the setting input to the regulator we fixed the following variable quantities: the phase current 
in the model (¢), the voltage across the rheochord Rrh (0), the operation of the null-relay contacts R(¢"), the opera- 
tion of the controller contacts (®) and the contacts of the intermediate relays 1SR and 2SR(®'), and the displacement 


Taking into account the imposition of a limitation in the sense of the sign, amplitude, and instant of appearance 
of the noise (cf. Section 2). 


** The experimental characteristics for the basic perturbations and the noise are shown in Fig. 4 for an actual ‘ 
electrical ore furnace for smelting lead. 
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t of the electrode motor, The oscillograms for the variations of the basic quantities in the compared networks are 
shown in Fig. 8 for identical conditions. 


It is evident that for a network without correction 


y— ee a ae (Fig. 8a) oscillations occur for the specified transfer co- 
A Bie efficient of the object. In a network with correction (Fig, 
eet eer?” ae ey 8b) the transient response is close to an aperiodic response, 


CONCLUSION 


mt ad San Bee The method studied in this paper for correcting the 
b dynamic properties of ACS is sufficiently simple and ef.- 
fective and can be used for many industrial control cir- 
Fig. 8. Oscillograms for the transient response: a) in the cuits. The equipment which is used can easily be com- 
network without correction; b) in the network with correc- bined with standard instruments and controllers. In par- 
tion ticular, the network in Fig. 3 can easily be constructed 
on the basis of the “IR-130" proportional -integral con- 


troller. 
The adjustment of the network is also simple. 


Engineer A. N. Ryndin participated in the experimental research. In conclusion the author thanks Ya. Z. Tsy- 
pkin for his valuable remarks and his help in preparing the paper for publication. 


APPENDIX I 


In order to determine the equations of motion for the ACS on the phase plane we shall study the equations for 
the individual elements (1)—(7) simultaneously while eliminating the intermediate variables in accordance with the 
method described in [1]. For an input deviation | ¢ | >] Ag] we obtain 


km, (t — T) 
Sout \— dat. 
at the output of the system. 


In order to fix the deviation 9% we have 


km,,t? + km,,tt, (15) 
cout = —— gp + 


after integration. 


If the ACS under study were not to have lag, then the cut-in time would be t, In an ACS with a delay 7 the 
variation of y, during the cut-in time t, is smaller by an amount km, te /T; however, for switching it is necessary 
that the cut-in time for the system be the same in both cases in order to liquidate the perturbation. Therefore, the in- 
crease in speed in both cases must be. determined by the equation 


km,,t, 


(16) 


which is derived on the basis of (5), (1). From this it is not difficult to derive the equation which determines the 
switching line for the ACS on the phase plane when lag and a dead band are present. For an input deviation 9, — Ag 
we have the equation for the switching line 2 (Fig. 6a) after substituting (16) into (15): 


— + Oo. (17) 
Analogously, we determine the equations for the switching lines 2', 3, 3'(Fig. 6a): 


The equations for the phase trajectories on the planes %,, &, 4, that characterize the acceleration and brak- 


ing modes when the servodrive is connected in different directions are also determined on the basis of the simultan- 
eous solution of equations (1)—{7). After the simplest substitutions [1] we obtain the equations for the phase trajec- 
tories on the planes ,, ®,, and 4 respectively, in analogous fashion: 


The phase trajectories in Fig. 6a are plotted for the following values of the basic parameters corresponding to 
the ACS for controlling the electric power to a lead smelting cannery 
k = ko 


2Ag = 0.2; k= 0.8; myo, = 1; 


heb = 20 cm ; T = 7.6 sec ; Mogg = 1.2; Gp = 400 cm ; Hing = 0-75; my, = 1.5; 
r 03 0 0 ke = 
t=0.3 sec; m .75; 
v2 bd cm 
hp, is the base value for the electrode depth (cf. Fig. 7). 


The drive motor for the electrode is an induction motor with a short-circuited rotor (type MA~-142 2/6 ); the 
rate of displacement in raising the electrode is 80 cm/min, and it is 90 cm/min for lowering the electrode; the brake 
is of the KMT type; the furnace transformer is rated at Py; = 3600 kva, k = 14,000 amp, short-circuit voltage 1, .. = 
=9.7%; the line Z,o, = 0.0012 ohm; the electrodes are graphite with a diameter of 700 mm. 


For Figs. 6b and 6c, %=0.1 sec., K= 1.4K, Ky=1.12, The remaining quantities correspond to Fig. 6a. 
APPENDIX II 


The Data for the Model and the Controller. We used an electrolytic bath (of the three-phase type; cf. Fig. 7) 
which was constructed as a model of a multiple-slag electrical ore furnace on scale 1: 10 with respect to an indus- 
trial lead-smelting furnace. The electrolyte consisted of an aqueous solution of table salt. The electrodes were graph- 
ite with a diameter of 5.5 cm, 


A transformer of the type TOS rated at 1.5 kw with three legs and 1, . = 4.8% was used; the nominal voltages 
for the primary and secondary windings were U,/U, = 220/21 v, and the nominal secondary current was ly) = 23 amp. 


The electrode drive motors are actuating mechanisms of the type IM 2/120 with a pulley; D=40cm, Gp = 
=1.50hm™, =4cm, ky =0,2 ohm™/cm, k=0.55. The controller dead band corresponded to AG =0,030hm™, = 
=0,0245 sec, T =0,08 sec, T =0.2sec, Ty = 0,04 sec, fi = 0.36; = 0.43. 


The limiting value of k for which there are no oscillations in the network with correction was equal to kj =1.2 
on the basis of the experimental data. Computations according to formula (12) yield ky *2k, i.e., ky = 1.1,which is 
close to k}. 


In our experiments the maximum perturbations at the controller input were equal to AG,,,, = 0.50hm™~ (or 
~33%Gp), She max =2-5 cm. 
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ANALYTICAL SOLUTIONS OF LINEAR EQUATIONS DESCRIBING ONE 
CLASS OF DYNAMICAL SYSTEMS WITH VARIABLE PARAMETERS 


Ts. G. Litovchenko 


(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 457-465, April, 1961 


Operational calculus is applied to obtain the solution of linear equations with linearly varying co- 
efficients in general form. Conditions under which this solution may be represented in closed form, 
i.e.,expressed in terms of elementary functions and tabulated special functions, are established. 


The mathematical investigation of the behavior of nonlinear automatic control systems in the neighborhood of 
stable motions, with first order accuracy, reduces to the analysis of linear equations with variable coefficients, the 
equations of first approximation for certain unstable motions of the system. A general solution to these equations, 
as is known, may not be obtained;and therefore it will be valuable to single out certain classes of equations which, 
on the one hand, have a practical value, and on the other hand, may be solved in closed form. This will permit a 
simple analysis of the automatic control systems corresponding to them. 


In this paper we examine one class of equations of the form 


v=N 
(t) = u(t), 


where Fy , is a certain linear stationary operator defined by 
t 
F,, z(t) = 2 (x) de. (2) 
0 
The structure of the system corresponding to equation (1) is represented in Fig. 1. In the particular case when 
fy (t —<) = &*+1) (t—+), equation (1) takes on the form of a differential equation: 


One deals with equations of type (1) or (3) in planning various linear automatic control systems whose variable 
parameters may be considered to the linearly varying. 


Examples of such automatic control systems are gyroscopes in those stages of motion when the kinetic moment 
changes approximately linearly in time, tracking systems with linearly varying driving moment of inertia, radioloca- 
tion systems in conjunction with automatic escort with variable band pass, and also the control systems of the type 
examined, for example, in [1, 2]. 


9. 
10, 
11, 

12. 


of 


Equations (1) or (3) are equations with singular points, and their solutions 


ult xt) may be represented by infinite power series (3], which, however, are not con- 
oad d 7 ut venient for the analysis of the properties of automatic control systems. There- 
; ay “Sy fore it is of interest to distinguish those cases when these solutions may be put 
1 into closed form, i.e. when they may be expressed in terms of elementary func- 
. tions and tabulated special functions. The fact that many particular types of 
equation (3) for N =2 have such solutions (Bessel 's, Laguerre's, the degenerate 
' | hypergeometric equation, etc. [3]) lends itself to this idea. 
4 The problem posed in conveniently solved by application of the Laplace 
(0,448) : transform. Applying it to the equation (1) or (2), we obtain, by virtue of the 
at linearity of the variable coefficients, an equation of first order, the solution of 
which is easily obtained by quadratures, Then we determine, firstly, the con- 
Fig. 1. ditions under which the transform obtained has an inverse, and secondly, we ob- 


tain the inverse itself. As a result we find all particular cases in which it is 
possible to solve in finite( analytic) form, and also find the solutions themselves of equations of type (1) independent- 
ly of the form of the operators Fy» and of equations of type (3) independently of the order N. 


We apply the Laplace transform to equation (1). Considering that inverses of the form (2) correspond to the pro- 
duct of transforms, and also the fact that the product of a certain function of time by t corresponds to differentiating 
its transform with opposite sign [4], we obtain the equation 


N 
— Fo (p) X (p) = (4) 


in whichE | (p) is the transform of an impulse transient function fy(t). In the case when equation (1) takes the form 
of differential equation (3), /',(p) = p’. 


After transformation, 
N 


N N 
— (p) X (p)= X (p)— x’ (P) (p) 
v=0 


v=0 


equation (4) may be written in the form 
(p) + P (p) X (p) = Q(p), 


(5) 
(P) Se F, (Pp) 
> (Pp) > (P) > (Pp) (6) 
v=0 v=0 v=0 


The functions F,(p), forcing functions of linear stead systems possessing together with variable parameters (a ,+ 
+b,t)the same structure as the system examined may be represented in the form of rational functions of the argument 
p. Therefore the functions P(,) as well is a rational function, with the order of the numerator for all @,=-0 and 
b,=£0 equal to the order of the denoininator. The latter fact follows from the obvious fact that the second term of 
the function P(p) for p-»*© approaches a constant limit equal to 


N N 
= lim] — >) a.F.(p)/ >} (p) 


v=) 


and the first term approaches zero, 


On this basis, the function P(p) may be represented in the form of a whole part and the sum of proper fractions: 
n Ma 
c 
P => Ss (7) 

where p,, are the poles of the function, N is the number of poles, and My the order of the nth pole, and c, », expan- 

sion coefficients defined by the formula 

My-—m 


enm = — Pa) (PN (8) 


dp 


uk, 
| 
2) 
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In the case when the functional equation (1) takes the form of a differential equation (3), the function P(p) is 
equal to 


N 
vb, — ap” 
P (p) = = N (9) 
> 
v=0 
In this,the whole part is equal to ¢, = lim P(p) = — >". 
N 


The general solution of equation (5) 


X (p) =e + dp | 
after substitution of the expression for P(p) from formula (7) takes the form 


X (p) R*(p) [C (p) dp, 
in which the functions Kp) and V (p) are equal to: 


nm d 
V (p) = R(p) Q(p). ag 

Simple transformations lead to the following formula for the function R(p): 

N My nN Ma 
n,m n,m 1 
(P) = xP i) 2 a P,)™ ap} a (P— Pp)™(m — 1) 
+ ena In (p— = TT (p— pa) I 
n=] n=1 n=] m=2 


We further determine the inverse x(t) of the function X(p). Using the convolution rule [4] and operational re- 
lations (27) established in Appendix I, from formula (10) we obtain 
t 


z(t) = (t — t) dt, (14) 
0 


where the functions r® (t) and v(t) are the inverses of the functions R™(p) and V(p). 


Thus formula (14) gives the general solution of the original equation (1) independently of the form of the opera- 
torsE,- In this the formula (13) shows in which cases the integral (14) may be taken in closed form, i.e.,expressed 
in terms of elementary and tabulated special functions. This is possible when the coefficients c, ,, of series (7) for 
m= 2 are equal to zero, and for m=1 they are actually whole numbers. In this case the function 


N 
R(p) = I] (19 

is a rational function; the functions R~*( p) and V(p), also rational functions, have inverses expressable by sums of terms 
of the form t*e*", Then integral (14) in taken in the general form (see Appendix III) and the solution x (t) is expressed 
in the form of some combination of functions te*'E;(t). If, however, the numbers c, ;, for m>2 are not equal to 
zero and at least one number cp y is not whole, then the inverse r(t) of the function R(p) may not be obtained in 
closed form, and the solution x(t) of the original equation (1) or (3) may not be obtained in such a form either. This 
follows from the fact that for cp m+ 9 and m= 2 the function R(p) contains a factor of the form e*/P Pa for which 
in inverse, on the basis of Appendix II and the displacement theorem [4], is expressed by Bessel functions of the form 


yee! (+) I, (2 V at) a However, the convolution integrals of these functions, defining r(t), are not expressible 
in terms of known tabulated functions [5, 6]. 


As a result one may formulate the following conditions under which it is possible to solve equations (1) or (3) in 
closed form. 


1, The function P(p), defined by formula (6) or (9), does not contain multiple poles; this guarantees that the 
numbers Cy m for m= 2 are zero. 


(9) 


(10) 


(11) 


(12) 


(13) 


1€ 


2, The function P(p) does not contain complex poles and their residues cy , are whole numbers. 


Assuming that the indicated conditions are fulfilled, we find the solution of equation (1) or (3) in finite form. 
If the transform of the perturbing function U(p) is a rational function, then the function V(p) defined by formulas (12) 


and (6), and also the function R™(p), are rational. This allows one to represent R*(p) and V(p) in the form of the 
following summation : 


K Ly 
k=1 l=1 (mt 


where pj, are the poles of the function R“*(p), K is the number of its poles, L, the order of the: ki kth pole, gy) the co- 
efficients of expansion defined by the formula 


S is the excess of the order of the numerator of the function R™*(p) over that of the denominator, B, are the expansion 
coefficients defined by the formula 


4 8 


P; are the poles of the function V(p), I the number of its poles, ©; the order of ith pole, qjg the expansion coeffi- 
clents defined by the formula 


gio = lim 


J the excess of the order of the numerator of the function V(p) over that of its denominator ,and A; the expansion co- 
efficients defined by the formula 


J 


The inverses r*(t) and v(t) of the functions R *(p) and V(p) are obtained, after applying to formula (16) the 
inverse transform formula, in the form: 


v(t) + 3 


Finally, the last step is to take the integral (14), in which we substitute r®*(t) and v(t) defined by (21): 


™—Co 
j=0 


k=1 l=1 * 


rit e1 t x) | 


is 
14) 
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(Pi 


(t—*)™ dt, 


Jin (t, m, n) = 


determined in Appendix III for p; — p, =a, we obtain the final expression for x(t): 


J 8 K ly 
1 

j=0 s=0 k=1 l=1 
Thus to obtain the general solution of equations (1) or (3) in closed form it is necessary : 
1) to compose the functions P(p) and Q(p) in formula (6) or (9), 


2) expand the function P(p) according to formula (7), (8) in a series and establish the possibility of obtaining 
the solution in finite form according to the conditions introduced above; 


3) compose the functions R™'(p) and V(p) according to formulas (15) and (12) and calculate the coefficients of 
their expansion in the series of the form (16) according to formulas (17), (18), (19), (20): 


4) obtain the final solution of the original equation (1) or (3) by formula (22). 


The method presented for solving equation (1) or (3) 
with linearly varying coefficients differs in that it depends 
neither on the order of the differential equation (3), nor on 
the complexity of the operators F) ; entering into equation 
(1). Since for a large number of terms of the series (for 
large values Cp) calculations by formula (22) are some- 
what complicated, in these cases an approximate method 

Fig. 2. of obtaining the solution may be recommended. This ap- 

proximate method is based on graphically calculating the 

integral convolution (14), determining the functions rt) and v(t) beforehand by the corresponding transformations( for 
example, by the known method of trapezoidal frequency characteristics [7)). 


Example. We shall examine a control system for the velocity of rotation of a body of variable mass. The sys- 
tem (Fig. 2) contains an amplifier of the mismatch signal 1, a rotary amplifier 2, a motor 3, and tachogenerator 4. 
The moment of inertia of the load 5 changes according to the linear law J =Jg+jt. The equations of motion of the 
system have the form 


(t) = u(t) —Q(t), M, (t) =A(t), JQ(t)+ /Q(t) = M,(t), J=Jo+it, (23) 


where the operator A "e5 ay to the transformation of the signal 6 (t) at the controlled moment My(t) according 
to the transfer function My (P) =! = the circuit “one and the transformation 


of the signal in the amplifier and motor. Initial data: K =0 A2kgm/ f= 0.01 kgm/ i= 0.06 kgm sec¥sec; 
Jo = 0.06 kgm sec’; G=1sec; =0.5sec. 


(Jo + jt) 2) + (F + A) Q(t) = Ayu (t). 


Comparing equation (26) with equation (1), we obtain 


@=1, b=0, Fo(p)=f+A(p), 
b=j, Fi(p)=p. 


Introducing the function 


24) 


We determine the reaction of the system to the 


signal of the form u(t) =1(1). 
1. According to formula (6) we compare the func- 
sec 
A (p) K 
Fig. 3. 2. We determine the expansion coefficients by 
formula (8) 
for 
K i i 
= =4 =—i, 
T; 
for 
For m=2c, 
3. We compare the functions R “(p) and V(p) according to formulas (15) and (12) : 
2 i 
We calculate the coefficient according to(17). In this,K=1, py Ly .,=4. 
619 = lim + 2)! + 2p(p + 2))——1, 


Flim [2(p +2) +2(p +2) +27] lim (6 +2] 


The coefficients B,=0. 
We calculate the coefficients g; according to formula (19). Here I=2, p;_,=0, O,.,=3; = -2, 


The coefficients Aj=0. 

4. We obtain the final solution according to formula (22) 


= (-4-# t+ (t + 1) — Bi ay + 


The corresponding graph is given in Fig. 3. 
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Appendix I 


oo 
We multiply the left and right sides of the relation F(p)=( e~?' f(t)dt by e~®P and integrate with respect 
to p. We obtain 0 


0 0 
From this follows the operational relation 
The quantity C may be chosen arbitrarily, since it corresponds to adjoining to the function f(t)/(t+a), having a 


discontinuity at t = -a, the term c5(t+a), defined only at the point of discontinuity. In the particular case when the 
function f(t)/(t=a) does not contain a discontinuity (for f(-a) + 0 and f"(-a) # ©), one must take C = 0. 


Appendix II 
We apply the inverse Laplace transform to the function on3 5 We obtain the operational relation 


J, 
for «<0, 
-—(4)"neva for a>0 


for a<0, 
(+)"nevay for 


Appendix III 


We calculated the function J (t, m, n) = 
pansion, we obtain 


e** (¢ — 


dx. Representing (t- T)™ in the form of an ex- 


ower 


m t 
J(t, m,n) = >} (27) 


The integral entering into (29) may be brought, after a series of transformations, into the form 


= ta) - “de = 


(28) 


—at, o=1 


eP > 
| for which the transformations 
lead to the operational relation 
| 


We calculate the definite integrals entering into (28): 


— dz = Ei (at — ac) — Bi (— ac) , 
—at 


v= 


Substituting them into (28), and afterwards into (27),we obtain finally 
m 
J (t, m,n) = (— | [Bi (at — acy) — Bi (— + 
»=0 
e=p+n v=e—1 
e=1 v=0 
Tables for integral exponents of this function are given, for example in [5]. 
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PROBLEMS OF NONLINEAR FILTRATION 
I. THE CASE OF ONE PARAMETER 
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Original article submitted October 11, 1960 


The problem is solved of the filtration of a random process (parameter) from a nonadditive mixture 
with another random process (the carrier signal for the parameter) under certain assumptions concern- 
ing both random processes. Operators are found that carry out the filtration; the applicability of the 

principle of feedback is shown in this realization. An example from radio technology is described. 


There has been a great deal of work done on the linear filtration of a random process from an additive mixture 
with another random process, starting with the fundamental articles by Kolmogorov and Wiener. In practice, however, 
cases often occur when the function to be determined is in @ nonadditive mixture with other signals. Thus the radio- 
technological methods of measuring various quantities (frequency, signal delay, etc.) are essentially to separate aran- 
dom variable from a misture with a signal of which it is the modulation parameter and with various types of interfer- 
ence, In a more common case, the elements of control systems are interconnected by channels in which the desired 
information is coded in the parameters of certain carrier signals. In the transmission process, the latter signals 


fluctuate and are mixed with noise. Since these random disturbances are nonadditively mixed with the desired infor- 

mation, methods of reproducing the information can hardly be considered to be in the realm of linear theory. Methods 
for constructing filters in the general case will be discovered in the theory of statistical solutions, where it is not nec- 

essary to make any previous assumptions concerning the structure of these filters, Complete statistical information on 
the carrier signal and possibly on the parameter is needed, however. 


In what follows, we will consider filtration to be the continuous measurement (estimation) of some random quan- 
tity which is a parameter of the received signal; in the general case the received signal is also random. In the present 
series of articles, we will solve the Bayes filtration problem when the a priori statistics of the parameter are known. 
Here we will consider the case when the signal depends on a single parameter with a Gaussian distribution. In subse- 
quent articles we will consider the cases of many parameters, and cases of sampling depending on the estimate. Since 
we are mainly interested in practical applications, we will not aim to give the rigorour mathematical basis of all the 
results obtained. 


1, The Design of a Nonlinear Filter 


Let the function y(t) and one of its parameters A(t) at times ty, tz, . . ., ty from the interval (to, t) take,respec- 
tively,the values ., Ynt = y and fAy, Az, .. .,A,} =A. Let their joint probability density distribution be 


P Yor Yq | Aas An) Po (Aa, An) == P(y|A) P, (A), 


where P(y|A )is the conditional density distribution of y for given A or the probability function [1,2] and P)@) is the 
a priori density distribution of the parameter. The filtration problem is to obtain a function A(t) from the selected 
values of y(t), where this function is the approximation to A(t), and is in some sense the best possible approximation 
to the real values of the parameter. In the theory of the solution, the closeness of this approximation is characterized 
by the “function loss,” which is a minimum for agreement between the estimate and the real value of the parameter. 
The most common is the square function loss, which leads to the requirement of a minimum mean square error. We 
can see that in this case, A” satisfies the vector equation 


i=1 
the known formal solution of which is the conventional mathematical expectationA for a given y 
(y| A) Po 
J Py Po(d) dh 


However, except in the additive case (y = A +n, where n = n(t) is the interference), the integration in (2) cannot 
be carried out, and so the structure of the optimal operator remains unknown. 


(2) 


We note that in practice, only the case of a sufficiently accurate parameter measurement is of interest. It is 
natural, that for a relatively low interference level, the error made in measuring the parameter when an optimal op- 
erator is used will actually be small compared to the width of the a priori distribution. If we also take into account 
the fact that the rate of variation of the parameter A(t) is usually small in comparison with the rate of variation of 
‘the signal y(t) in which this parameter is coded, then Pty |A)can be represented approximately in the form 


P(y|a) ~exp{—L(’) +2 (3) 


AL (2° 


is a square matrix, and a" denotes the transpose. Thus the 


is a columa vector, A =|aul =| 
probability function, as a function of A, is approximated by a many-dimension Gaussian distribution near the point 

2» with coefficients depending on the realization y(t). The basis for such a substitution is the convergence of the loga- 
rithm of the probability function as a function of the parameter A(t) to a parabola, near the real value of Ag(t) for time 
intervals considerably longer than the correlation interval of y(t), but shorter than the time of variation of A(t). By 
assumption, such time intervals actually exist. It is clear that the expansion of the logarithm in a series can be car~- 
ried out to a point close to A9(t) that gives the same parabola. 


The approximation (3) makes it possible to construct optimal filters for a whole series of important cases of 
parameter statistics. Of prime interest is the Gaussian distribution for A: 


P, (2) = det” R exp 


where R =IIR(ty, ty)ll is the correlation matrix of A(t), V = the inverse matrix,and ={X,, ., An} 
the column vector of mean values, 


If we substitute (3) and (4) in (1), we obtain the result: 


n n 

= ene am (di — ] + (5) 
k=) t=} 

where Cp} is an element of the matrix C , which is the inverse of the matrix W+A in the interval (to, t,). The 

equality (5) related the estimate at the last instant of time t = t, with its values at all the preceding instants. In (3), 

(5), etc., the estimate of %,, for the times t,{k = 1,. . .,n) is formed according to the realization of y(t) in the inter- 

val (tp, t,,), which corresponds to the requirement of physical realizability of the operators sought. If the choice of 


y(t) is discrete then (5) is already the solution to the problem, while is the realization of y(t) is continuous, we intro- 
duce new functions by the relations 


tyt Ata y+ At 
C (tn, te) = Z(t) dt = 2, dt ds A(t, 8) = (6) 
ty— At/a 


and transform (5) from a summation to an integral by carrying out the limiting processes tj,,—t; = At 0, n+ @ 
(i =1,2,.. ., nm; nAt = t-te), which yields 


a(t) = (C(t, +\ A(x, 8) {4° (s) —K(s)]ds} de + (1) 


The physical interpretation of the relation (7) is shown in the block diagram in Fig. 1. The given realization 
is led to the element 1, which forms the quantity Z(t) equal to the mismatch between A(t) and %(t), The element 1 
can thus be appropriately named the discriminator. Simultaneously, y(t) is lead to the elements 2 and 3 for the for- 
mation of the impulse reactions C(t, T) and A(t, T) of the linear variable of the filters. The output of the element 
2 is added to the quantity A(t), and the result led to the elements 1, 2, and 3 for the adjustment of their parameters. 
The output of 2 also passes through the filter 3, is combined with the output of the discriminator 1, and the result is 
again led to the filter 2. This system is essentially a self-regulating system, the characteristics of the elements of 
which depend on the form of the realization y(t). It can be seen that the variations of the characteristics of the fit 
ters 2 and 3 tend to decrease the influence of the filtration component at the output of the discriminator, the strength 
of which increases with increasing mismatch €= A - %. Such a component occurs both in basic radiotechnical 
measurements of a parameter and in the system we are considering. When the noise is strong enough, however, this 
component does not play an important part. We can therefore reject the formation of the filters 2 and 3 by using the 
given realization, and using instead of A(t, T) its mean value obtained from the over-all properties of y(t), known to 
us a priori. In this case, the dotted line in Fig. 1 can be omitted, and we obtain a considerable simplification of the 


scheme of solution. If we use the assumption concerning the smallness of the mismatch and the rate of variation of 
the signal, we can set 


_ LG) = K8(t —*), 


where the bar denotes the mean over the set of values of y(t), and K is a coefficient which is determined from the 


coding of the parameter in the signal and the statistics of y(t), and which, under rather general conditions, is indepen- 
dent of the time and the measured parameter (see section 2). 


When we solve (7) for % (t), we obtain 
t 
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where the impulse reaction G(t, T) from (7) and (8) satisfies 


t 
C(t, d+ s) G(s, t)ds = G(t, (10) 


Fig. 1. Block diagram of a nonlinear filter: Fig. 2. Simplified single-loop system: 
1) discriminator; 2) filter with impulse reaction 1) discriminator; 2) smoothing circuit. 


C(t, T); 3) filter with impulse reaction A(t, T) 


Since © =[V +AJ“=[V+KAd J", where I is the unit matrix, then C(t, T), in its turn satisfies the equation 
t 
C(t, s) R(s, t)ds = R(t, (11) 


A similar relation also holds for the discrete case. 


The expression (9) shows that under the assumptions made, an optimum nonlinear filter is a single-loop servo 
system (Fig. 2) with a nonlinear discriminator which isolates the derivative from the logarithm of the probability func 
tion for A(t) at the point %(t). The system is closed by linear smoothing units with impulse reaction G(t, T) depend- 
ing on the correlation function R(t, T) of the parameter, and through the coefficient K on the statistics of y(t). It will 
be shown, in what follow:. that the smoothing units are linear only for Gaussian parameter distributions. 


~ In general, Z(t) depends nonlinearly on € = A - %, but if the probability is high that the quantity € (t) issmall 
then for the analysis of the accuracy of measurement we can set 


1 ab 
= E(t) + Ke (ty), 


where € (t) is the fluctuation component, which is not dependent on K5(t,-t,). Here the spectral density of the fluc- 
tuations € (t)/K, converted to a value of the parameter, is evidently equal to 1/K, so that K characterizes the signal/ 
interference ratio in the discriminator. We also note that the quantity 1/KT is the dispersion of the estimate of the 
probability maximum of the parameter A, which can be obtained from the realization y(t) after a time T, if the para- 
meter is constant during this time. 


From (9) and (12) it is easy to see that KC(t,7) is the impulse reaction of some linear system which is equiva- 
lent to the system in Fig. 2 for small mismatch. To the input of the equivalent system is lead the total “signal” 
A(t) - A(t) and the “interference” £(t)/K converted to a value of the parameter. Thus with the neglect of a fluctua- 
tion component proportional to the mismatch, for small €(t) the parameter and the interference ate additive and in- 
dependent to an equivalent degree, If we take into account the fact that the parameter has a Gaussian distribution, 
then in the light of what has been said, it becomes clear why the equation (11) coincides with the equation for a 
Wiener filter, for the isolation of a parameter with correlation function R(t, T) from a mixture with white noise witha 
spectral density 1/K for a finite time of observation (to, t). The main difference from the Wiener filter is the presence 
in the scheme in Fig. 2 of a doubly nonlinear discriminator element (for the input signal y(t)) and the necessity in 
principle of the feedback, which shifts that part of the discriminator characteristic that is linear with respect to the 
parameter, immediately following a variation in the parameter. 


When we use (9)-(12), we obtain the following expressions for the mean value and the correlation function of 
® (t): 


W(t) (tr) — (42) — = R(t, —€ (4, 4). 
Finally, the error of measurement of the parameter has the variance 


(t) = = C(t, t). (13) 


2. Discriminator Characteristics 


The results obtained above determine the operations carried out by the discriminator in the formation of the 
function Z(t), this function being a measure of the current mismatch and also of its basic characteristic, i.e.,the quan- 
tity K; but in practice these operations depend on the statistics of y(t) and the method by which they are coded in the 
parameter A(t). Certain general considerations can, however, be introduced, that make it possible to use a model for 


a quick determination of the operation of the discriminator and of its characteristics in the most interesting practical 
cases, 


If we ignore difficulties concerning the limiting transitions, the logarithm of the probability functional can al- 
ways be written in the form 


Lin. (y, = A; t) dt, (14) 


wherel(y, A, T) depends in principle on the realizations y(t) and A(t) in the interval (ts, 7), but in actuality only on 
the values of y(t) and A(t) in an internal of the order of the correlation time 1, of the signal y(t) preceding the time 

tT, According to the assumption concerning the relation between the rate of dhenge of the signal and the parameter, 
during this interval A(t) maintains an almost constant value which determines the function I(y, A; t) in the interval 

(t- T,, t), but has no influence on its value for other time intervals. On the basis of these considerations, it can be 
seen that Z(t) is statistically equivalent to the derivative of -2(y, A; t) with respect to A for a constant value of A 

equal to the value of the estimate X(t) for the given time; i.e.,the operation of the discriminator for optimal measure- 
ment reduces to the formation of the quantity 


Z (t) al (y, t) 3 (15) 


from the values of the signal realization y(t) that have been taken, and from the current estimate * (t). 
The quantity K characterizing the optimal discriminator is here given by the formula 


T 
4 A; t) 
K = lim \ at. (16) 
0 
The operation (15) can be realized physically in several ways. First of all, in certain cases, the derviative of 
Uy,A;t) can be formed exactly. Secondly, in practice we can replace the derivative by the difference 


i,e.,in the measuring device we use two channels which form Uy, A; t) at distances of #AA/2 from the value of the 
parameter being considered, We note that each of these channels is an ordinary optimal receiver for the signal y(t, A) 
without final integration [2]. It is also possible to use only one channel which is successively shifted in both directions 
from %*(t). All approximate methods must be evaluated relative to the characteristics they give compared with the 
characteristics of an optimal discriminator that forms Z(t) accurately. 


Since we have not yet chosen any definite parameter form or any system by which it may be coded, we will now 
give some examples of the representation of the logarithm of the probability functional in the form (14). The realize 
tion of the operation for the formation of Z(t) by using (15) or one of the approximate methods will be almost obvious. 


A. A regular signal v(t, 4) in white noise with a spectral density Ny. Here 
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T 
1_( (a(t, 
L(y, A; t) = sy, ly(t)—v(t, K = | 


There is no difficulty in generalizing this to the case of correlated noise. 


B. A noncoherent signal in white noise. The case of a periodically modulated signal with independent, random 
phases of a high-frequency filler in each modulation period T, is sufficiently well known [2]. For sufficiently 
large values of the signal/noise ratio, the quantity [(y; A; t) coincides with the voltage envelope at the output of the 
adjusted filter [2], and for a symmetric signal spectrum the quantity K is given by 


where u(t, A) is the amplitude and #(t, A) the phase-modulation of the signal. 

C. A fluctuating signal in normal noise. In practice, the most interesting case is that of a normal, fluctuating, 
high-frequency signal, with arbitrary amplitude modulation u(t, 4) and phase modulation (At, A), and with a correla- 
tion function equal to Ret, te) = u(ty, 4) (te, A) cos — te) + P(t, A) — Plte, p(t, — ts), 
where p(t, — t,) = x { S,(@)cos@(t,;—t,)d@ is the correlation-signal envelope in the absence of modula- 


tion and Sc(w) is the corresponding spectral density. If the function p(t) varies slowly in comparison with u(t, A)and 
¥ (t, A), then neglecting terms that are independent of A, we can write 
L(y, 4; = 


Here the function h(t), representing the impulse reaction of a filter, has the square of the absolute value of its 
frequency characteristic given by the expression 
/ 2No 


8 
and an arbitrary phase characteristic, u(t, 4) is normalized so that Jim | ut (t, A)dt = 4, and Warh is an 
0 


arbitrary intermediate frequency. The quantity K for this case is given by 
T 


T +00 


—oo 


3. The Characteristics of Smoothing Circuits 


Although the methods of designing Wiener filter are known for several cases (see, for example [3]), the repre- 
sentation of the filters in the form of servo systems can be used to obtain the known results by new methods, and also 
to obtain some new results, 


A. The parameter a linear combination of known functions with random coefficients Let A(t) = p f(t) + Ar), 
where f(t), X(t) are known functions and p is a normally distributed function for which fi = 0, y"=0" so that 


R (ty, te) = 3°f (t1) f 
The relation (11) is now an integral equation with a degenerate kernel, having the solution 


C(t, = SEO 
1+ /*(s)ds 
i, 


Thus the error in measuring the parameter 95, (t) always tends to zero for increasing t: 


1+ Ko jr (s) ds 


The impulse reaction G(t, T) of the smoothing unit, according to (11) is given by 
G(t, = (t) 
1+ f*(s) ds 


The smoothing consists of multiplying the output of the discriminator Z(T) by the function f(T)/ 


I[1 + Ker p(s) ds] (Fig. 3) and then integrating. As a result, we obtain the value of the coefficient », which 


can be used as a multiplier for the law of variation of the parameter f(t) generated by the system. After adding A(t), 
we obtain an estimate of X€). 


The gradual “switching off" of the discriminator is typical of this setup, and the 
larger the value of K, the more rapidly this takes place. The law of variation 


gud 3 HEH A(t) is known except for a multiplying factor, and so a finite realization interval 
is sufficient to predict the behavior of the parameter at all succeeding times. 
Fig. 3. The smoothing unit for We now consider the general case where i (t) = > phe, 
with f(t), . .. fn(t) known linearly independent functions and Hy, random coeffi- 
2) ator; 3) multiplier. cients for which = 0, (i, j= 1,... , m), and we find that 


R(t, t) = =F OSHC), 


where f(t) = {f,(t), . . . .f_(t)} is a column vector and $= | S4j| is a square matrix. The solution (11) has the form 
where Then from the equation (10) we obtain 


G(t, t) (t) [I + ST (21) 


It follows from (21) that the smoothing unit consists of m parallel channels. In the kth channel the function 
Z(t) is multiplied by a function depending on all the functions §(T) (i = 1,... , m) and on K, the resulting function 
is integrated, and the result of this integration is multiplied by f(t). All m outputs are then added, together with A(t). 


B, The parameter a stationary process. This case has been more thoroughly investigated than any other, be- 
cause of its great theoretical and practical interest. For a finite time of observation, C(t, 7)is obtained by methods 
described in [3]. If the time of observation exceeds the correlation interval of the parameter, then if we let t, tend 
to and try to find C(t, T) = C(t-T), the relation (11) yield the well known Wiener-Hopf equation 


The physical realization of the solution (22) has the form [3] 


+00 


C (iw) ={ C(t) (@) = R(t)e dt, ¥ (iw) ¥(—io) =14+KS(@) (24) 
0 


is a special representation (factorization) as a product of factors with zeros and poles in the upper and lower complex 
half-plane,respectively. Such a representation is possible everywhere, for example if Sq) is a rational function of w* 
The operation [ ], is the separation of that part of the expression in brackets that has a pole in the upper w half-plane; 
analytically, this is expressed by using the two integrals in (23). 


C. The parameter is a Wiener process, i.e. the integral of white noise, Its correlation function is 
R(t, t) = Bmin(t, t) (B = const). 
In this case, when we differentiate (11) with respect to T we obtain 
| 
t)— KBC (t, rt) = 0, 
the solution of which is 


Cit, 1) = (t) = BR t (to = 0). 


in the region T< t. 


After a period of (t<(BK) 44), of the regime with o6,¢ (t) increasing like the variance for the parameter Bt, a re- 
gime is established in which 


Sout (t) lim = (25) 
From (10), the impulse reaction of the smoothing units is 
G(t, = V B/ (26) 


Because of the previous knowledge of the continuity of A(t) and the zero initial error of measurement, the out- 
put of the discriminator for T< (BK) ¥2 is suppressed. For T > (BK) ” the smoothing becomes a constant amplifica- 
tion and integration, and C(t, T) tends to 


V B/K exp {(— V BK (t — »)}. (27) 
Thus the use of one integrator in the smoothing circuit of a servo system, often occurring in practice, is optimal 


for a parameter in the form of a Wiener process. 


D. The parameter a process with stationary growth. In C (para, 4), for a nonstationary process with a long ob- 
servation time, the filter appeared to be particularly simple in the case of the class of circuits with constant parame- 
ters. 


Since a Wiener process is a particular case of processes with stationary parameters [4], it is natural in the case 
of stationary parameters to seek general algorithms for obtaining limiting operators. 


For a process with stationary parameters of the kth order, the statistical characteristics of the kth difference 


k 
= 
are independent of the time. The spectral density S(w) of such a process is defined by the relation 
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[A¥, (ty) — (4) (4) — Aba = = (28) 
-+00 


The transition in(27) to the limit for 1,, T, +0 gives the relation between S(w) and the correlation function 
P(t, - t,) of the kth derivative of-x(t) 


Jim wrth) = P(t —t) = (an 


In (11),we replace the function R(t, T) by its expression in terms of P (t; - t,): 
t 
Rt, = jan dt, (t —t,)"* ty)" p(t, — 


We further set tp = -”,and seek C(t, T) in the form C(t - T). When we take the Fourier transform and use (29), 
we find that C(iw) is expressed by the same formula (23), oe factorization (24) involves the representation of 
the factor w* in the denominator of S(w) in the form (iw)*(-ia)*,where 1/(i w) is formally assumed to be a pole in 
the upper half-plane and 1/(-iw) a pole of the lower w half-plane, For the justification of this method we can pro- 
ceed as in [3], and extend the class of functions differing from zero on the half-axis (0 = t < + ©), by the addition of 
the functions 


t>0 
and introduce similar function on the half-axis (-~ <t= 0), For example, in the case of the Wiener process, 


(m = 1, 2,...) 


1+KS(o)=1+4 VBK — iw 


—iwo 
and so from (23) we have 


B i 
VBR + io |, 7) K VBR+io 
and this agrees with the Fourier transform (27). 


For another example, we can consider the process obtained by taking the double integral of white noise with a 
spectral density D(S(w) = D/w), It can be seen that 


iw + (KD 4)" 
(iw) = (io)? + (AKD)" 


(30) 


and 80 o% yr is equal to (4D/K*) 4. In this case, the smoothing units have the transfer function 
V 1+(4/ KDy"p 
K P 


ie, they consist of a double integrator and a correlating RC-circuit with constant time (4/KD)%. Such smoothing 
filters are very often applied in practice, but the conditions for their optimization are still unknown, 


EB. The mixed case, It is sometimes important to represent the parameter as a combination of two terms~a 
linear combination of known functions with random coefficients and a certain nondegenerate process (for example, 
with stationary growth) with correlation function r(t,t). Then 

m 


Rit, 1) = >> Suh) +r 
4, 


where we will assume that r(t, t) <2 Sisfi (t) 15(2)- To a first approximation, the equation (11) takes the 
same form as in A (para, 4), with the solution C,(t,7) given by (20). For a second approximation C,(t,T) we have 


t t 
C3(t, *) + K\¢, (t, s)r(s, t)ds = r(t, t) — (t, s)r (s, t)ds. (31) 


When t— te ~ ©, the term containing C; (t,t) vanishes in (31) as well as in all subsequent approximations to 
C(t,T). The smoothing circuits in this case can be approximately realized by two parallel channels: a channel for 
estimating the linear combination, which is gradually disconnected from the discriminator, and a channel for main- 
taining the nondegenerate process, The outputs of these channels are combined. For processes with stationary growth 
in a steady regime, the second channel contains circuits with constant parameters. 


4, A Practical Example 


We consider a radio transmission of random information coded in the delay of a harmonic frequency modula- 
tion, When there are fluctuations, the correlation function of the signal (see B of para, 3) is 


Ry (ti, ts) = p (ti — te) cos + (sin —2) —sin (44 — 9) | 


where w mand Q are the amplitude and frequency of the modulation, and \ is the measured delay, According to 
(15) and (17), the operation of the discriminator is described by the formula 


t 
Z(t)= \ h (t — c08 @in¢(t — y (t)cos (on + | ar} x 


t 
x{ h 608 v(x) c08 2 — 2") sin | (09 + + | ax} 


and an exact realization can be obtained in the form of that part of the block diagram in Fig. 4 within the dotted 
lines, The signal is combined in two mixers 1 and 2 with the heterodyne signals from the heterodyne generator 16 
with a phase shift of #/2, obtained by the phase-changer 15, and with identical frequency modulation with delay 
time corresponding to the measured value of \* and imposed by the control delay line 17, In one of the channels 
there is also a mixer 3 with frequency ChM, which transmits to the intermediate frequency wjnterm the ChM side 
bands that arise at the output of the mixer 2 for the mismatch between \ and d° with a strength proportional to 

X— 2X. There are also filters 4 and 5 with frequency characteristics (18) tuned to the intermediate frequency. In 
the phase detector 6, the error Z(t) is obtained from the filter outputs. The quantity K which characterizes the oper- 
ation of the discriminator, in the case of a signal with a rectangular fluctuation spectrum of width Af, is, according 
to (19), 


(32) 


where h = P,/2NeAf, is the signal/noise ratio in the signal band and P, = p (0) (w) do. 
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Fig. 4. Block diagram of the delay-measuring unit. 
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Fig. 5. Dependence of the delay- measurement error on the time 


Let the measured quantity vary according to the law 4(t)=X(t)+pt+€(t), where ACO is the known 
mean value, ) a random, normally distributed variable (yu = 0, es o”), and & (t) a Wiener process with variance Bt, 
Then the impulse reaction of the smoothing units, according to(10) and (11) is 


sinh Y exp{ \ 


V BK \~ 


and it is realized approximately by the remaining part of the block diagram in Fig, 4. The discriminator output is 
led to the tracking channel and the channel for evaluating the velocity. The channels have controllable amplifiers 
7 and 8, and integrators 9, 10, and 11, 


The tracking channel is gradually triggered, while the estimation channel is gradually disconnected from the 
discriminator, The element 13 controls the amplifier, and the programming mechanism 14 introduces the mean 
value which is summed in 12 with the outputs of the two channels, The output of 12 is the estimate of A(t), and is 
led to the discriminator, thus closing the system. 


By using (11) and (13), we obtain the variance of the measurement error 


sinh_sds 


(t) = VB) eothy BKt — 1 BRe + (34) 


where VB/K =o3y, (00) is the steady variance, and v = (B°K)/o* 
is the ratio of 0% yt (e) to the variance of the linear term pt through 

the interval (BK), equal to the time for the establishment of a 
steady state in the system and calculated only for the Wiener process. 
In the steady smoothing regime, the system becomes a simple inte- 
grator with the introduction of a mean value, and for Gpy,;(), from 
(32) and (34) we obtain the formula 


Fig. 6, The dependence of delay- measure- 
ment error on the parameter variance and 
the signal/noise ratio, = | 


B i+h 
Ajo, 


Figure 5 shows the time dependence of 6%,,1(H / 6%, (0o))for some values of v=(B*K}Yo?, obtained from (34). 
The dotted curves show the geometrical location of the maxima] values and the asymptotic value. Figure 6 shows the 
relation Of OpyrwW;, @nd_ the ratio of the growth B/4f, of the a priori variance in the correlation interval of the 


(35) 
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signal to the square of the ChM time constant 1/w{p for various values of h, the signal/noise ratio, 


Many other examples could be given, showing the feasibility of designing synthesizing, nonlinear measuring 
and filtering apparatus. 
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THE ERRORS CAUSED BY A MEASURING APPARATUS, 
CHECKING THE VARIATION OF A VARIABLE QUANTITY 
AT EQUAL TIME INTERVALS 


Vv. A. Bolkonskii 


Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 479-482, April, 1961 
Original article submitted June 23, 1960 


We consider an apparatus for measuring the values of a variable quantity at discrete instants of time 
and signalling when this quantity has exceeded a certain allowable level, An estimate is obtained 
for the mean sum of the durations of those time intervals during which the controlled quantity is 
beyond the permitted level, that are not observed by the discrete control, 


During the course of a technological process, it is often necessary to ensure that certain variables characteriz- 
ing this process (for example, the supply of heating gas for a blast furnace, or its temperature and pressure) do not pass 
beyond certain limits. This result is obtained by determining the value of the variable at certain fixed intervals of 
time (we will call these intervals measurement cycles), and if it is observed that the measured value passes outside 
the allowable range, a signal is transmitted indicating the necessity for regulation. 


These excursions beyond the allowable limits (we will call them overshoots, taking over the terminology in 
radio- technical literature) can occur, however, in the intervals between measurements: and they will then not be ob- 
served by the controlling apparatus, It would therefore be interesting to find out how often it is necessary to take 
measurements, to ensure that the time of duration of the undetected overshoot is sufficiently small. It is also natural 
to take into account the fact that for a high frequency of measurement, the undetected overshoots will have only short 
durations and can not strongly influence the course of the technological process, We will therefore pose the problem 
in the following way. 


We will assume that the variation of interest to us is that of a normally distributed random process with zero 
mathematical expectation and known correlation function rr. 


We will find the mean sum of the durations of the undetected overshoots beyond a level 5 for a long time of 
observation T, with a given time A between successive measurements? We denote the mathematical expectation of 


*In a similar way, we can take into account the fact that undetected overshoots can only exceed the level 5 by small 
amounts, and it is possible, for example, to derive the dependence of the mean number of undetected overshoots, ex- 
ceeding not only the level 5, but also the level 5(1+ €) on €(€> 0), In this form, the problem can be applied in an- 


swering a question that is important in metrology: At how many points of a scale must a measuring apparatus be 
checked, in order to ensure the correctness of readings taken by it? 


al 


the total duration of all overshoots (detected or not detected) by M,, The effectiveness of the control for a given fre- 
quency of measurement can then be measured by the ratio M,/Mj. 


It is obvious that with the above assumptions, we have 


M, 


where P(x, > 5) is the probability that a value x, taken by the process exceeds the level 5 , tp is the variance of the 
process (the value of the correlation function r, for t = 0), and 


We denote by » the mathematical expectation of the number of overshoots starting in a unit time, and by m(x) 
the conditional (with the condition that the overshoot begin at the zeroth instant of time) mathematical expectation 
of the quantity Tx which is equal to the duration of an overshoot if it does not exceed x, and equal to zero if it does 
exceed x. It is evident that the mathematical expectation of the total duration of overshoots beginning during the 
time interval [t, t+ dt], 0<t<A, that are undetected by the measurement made at the instant A, is pdtm(4- 0). 
Therefore, 


A 4 
M, = = (2) 
0 0 
(we take T/A as integer), As is known (see, for example, [1] 
— (3) 


where rg is the value of the second derivative rf for t= 0 (-r§ is the variance of the process x}). In the calculation 
of the function m{x) we will assume that the second derivative x? of the process x; remains constant during the whole 
overshoot, equal to its value at the beginning of the overshoot. — 


This is a natural assumption, as can be seen from the following considerations. It is obvious that the measure- 
ment cycle 4 must be chosen so that it can be expected that the majority of overshoots will be detected. This is pos- 
sible only in the case when the main part of the spectrum of the process consists of frequencies that are considerably 
lower than the frequency of measurement. Since only overshoots with durations not greater than 4 are considered, 
then the probability is small of a sharp change in the quantity x, or even in x? during the overshoot. In any case, 
this assumption becomes all the more valid when A — 0, 


With this assumption, the duration of an overshoot beginning at the zeroth instant is given by the expression 


(4) 


m (z) =\\ z) dydz (:<0,%<z), (5) 


where p(y,z) is the conditional distribution density of the vector (xj, xf ) under the condition that x= 5. When we 
carry out the calculations, we obtain 


—oo 


where 


% 
Therefore, 
ort 
x- 2 


where 


1@= CG 
—= for G<0 


To simplify this expression, we note that the quantity M, can be assumed to be small in comparison with T. We 
may therefore assume that 5 2 2, 


It can be verified that if the spectrum of the process is concentrated in a frequency band from 1% to 2% of the 
integration frequency, then the expression in the square brackets in (6) is close to unity (for example, if the frequen-' 
cies in the spectrum de not exceed 2 9/54 and the spectral density does not increase with increasing frequency, then 
this expression is between 0.8 and 1.2), If we note that the assumption that the process is stationary and has a normal 
distribution is usually valid only to a first approximation, and also that the statistical character of the process is not 
known accurately, it is obviously not usually sensible to use the complicated formula (6). It can, in fact, be replaced 


by the approximate expression 
M,=T 0.026 (8) 
If we also use the approximate relation 
=z, 


which is accurate to within 15% for x = 2; then for 7 =2 we obtain 


M ra\ & 
Example A, A process with spectral density 


54 
x 


and correlation function 
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and 
i r. 8 
p(w) = 
t 


M,T MJM, 
A B A B 
2 | 0.023 6.7-10-* 5.1-10-* 0.029 0.022 
3 |0.00135| 8.2-10-* 6.2-10-* 0,06 0.047 


The results of this example show that if the spectrum of the process is in a band with frequency about 1/5 of 
the measurement frequency, then the total duration of the undetected overshoots is on the average 2-@ of the total 
duration of all overshoots beyond the prescribed level. 


SUMMARY 


1. A method is obtained for calculating the fraction undetected of the total duration of overshoots of a con- 
tolled variable quantity beyond a permissible level when the time interval between successive measurements is 4, 


The mean sum of the durations of all overshoots and the mean sum of the duration of undetected overshoots 
occurring during a time interval of length T can be calculated for a measurement cycle 4 from the formulas (1) and 
(8). 


2. If, with sufficient accuracy, the process can be assumed to be stationary and to have a normal distribution, 
and if the variance of the process and its derivatives (rp,rg, r'')) are known with sufficient accuracy, then the calcu- 
lation of M, can be carried out with considerable accuracy by applying numerical- graphical methods of integration 


and the formulas (6) and (7). 
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THEORY OF A SINGLE-CORE MAGNETIC AMPLIFIER 
WITH RECTIFICATION OF THE LOAD CURRENT BY MEANS 
OF A SYMMETRICAL NONLINEAR RESISTANCE 


Dai Tse-Hsin 


(Moscow) 
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pp. 483-492, April, 1961 

Original article submitted June 21, 1960 


It is demonstrated that when a nonlinear resistance is connected in series with a single-core saturable 
reactor, the current in the ac circuit contains a dc component which is proportional to the dc magne- 
tizing signal, 


An expression is derived for the dc component of current in the ac circuit; this expression establishes 
a relationship between the dc component in the load circuit and in the input signal, 


We know that there are magnetic amplifiers which rectify the even voltage harmonics by means of symmetri- 


cal nonlinear resistances. Push-pull magnetic amplifiers basedon this principle have a lower sensitivity threshold of 
the order of 107!4-10" [1). 
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However, the theory for magnetic amplifiers of this type has not been developed sufficiently. Therefore, we 


have posed the problem of clarifying the physics involved in the operation of amplifiers designed according to the 
principle of combining ideal saturable reactors and symmetrical nonlinear resistances (Fig. 1). 


Fig, 1. 


In this paper we make the following assumptions. 


a) The characteristic for the symmetrical nonlinear resistance is approximated by 
three straight-line segments shown in Fig. 2: 


du 
r =F = 00° for |u|<|uol, (1) 


du 
r=%=0 for 


Such a volt-ampere characteristic can be obtained by connecting semiconductor 
diodes in a bridge circuit as shown in Fig. 3; here ordinary diodes are connected in the 
arms, and a reference silicon diode (stabilitron) of the type "D808-813" is connected in 
the diagonal. If we use silicon diodes as the rectifying diodes, then the volt-ampere 
characteristic for such a bridge is very close to the characteristic shown in Fig. 2, Even 
at a high temperature (of the order of 90° C) a high reverse resistance is maintained in 
diodes of this type. The forward resistance of the diodes, as well as the reverse dynamic 
resistance of the reference diode, can be treated as a resistance connected in series with 
the bridge. The reference diode does not affect the symmetry of the volt-ampere char- 


acteristic for the bridge, since the voltage across the reference diode connected in the diagonal always bucks the 
voltage applied to the bridge. Carborundum resistances are good as far as stability is concerned, but they have high 


losses, 


b) the magnetization curve is represented by three straight-line segments (Fig. 4): 


=const for |B\<|B,|, 


aB (2) 


3 c) The network is supplied with a sinusoidal voltage. 


Up 
Fig. 2. Fig. 3. 
S 
Ne 
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H. wf 
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Fig. 4. Fig. 5. 
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d) The current in the control circuit does not vary for variations in the ac circuit, This is achieved by con- 
necting a large ac impedance in the circuit (in particular, a large inductance). 


e) In this paper we shall study only the steady-state mode, 


f) A resistor is connected in series with the symmetrical nonlinear resistance and the faturable reactor in the 
load circuit. 


The network which is studied in this paper contains only one core, However, the results obtained here can 
be used for combining single-core networks to form multiple-core networks, 


The Operation of a Saturable Reactor with a Symmetrical Nonlinear Resistance 
For the network in Fig. 1, it is expedient to examine three modes of operation. 


Mode 1, This mode occurs when the magnetizing force produced by the maximum control current is less than 
or equal to the magnetizing force corresponding to the knee of the curve in Fig, 4 (i.e., We [|=] Wei! , Fig. 5). 


This mode is of greatest interest. It shall be demonstrated that there exists an approximately linear relation- 
ship between the dc component Ip of the current i, in the load circuit and the control current I... 


We shall study how the voltage U, across the saturable reactor, the flux © in the saturable reactor, the voltage 
across the nonlinear resistance U,,,and the current i; in the load current vary (in Fig. 6 ly = 0). 


When there is no signal during the interval 0 =t =t, the supply voltage U,,sin wt is less than the reference 
voltage Up across the stabilitron, and all of the supply voltage appears across the nonlinear symmetrical resistance. 


At the instant t= t, at which Up )sinwt,= Up, the current i, begins to increase according to the equation 
Umsin ot —U, = + iy, = + (3) 


where L, is the inductance of the winding W, of the unsaturated core, r, is the resistance in the load circuit. 
Keeping in mind the fact that wL, > 1, it is possible to neglect the voltage drop iyr;. Then from Eq, (3)we 
have 


W10*D ~ Um sin wt —U,. 


The flux varies according to the expression 


t 
AO = \(Umsin at —U,) dt. (4) 
Henceforth we shall assume that the values of U,, and Us are such that the quantity 4@ defined by formula(4) 
is equal to 4s. Then, if the flux varies from zero, we obtain ® = %, for t= t. 
The current i, varies according to the expression 


i, = (5) 


For t= t= t the voltage across the nonlinear resistance does not vary and remains equal to Up. 


At the instant t+ t, the supply voltage again becomes equal to Up, the voltage across the saturable reactor de- 
cease to zero, the flux reaches and value 5, and the current i, reaches the value I, (Figs. 5 and 6), 


For t> t, the voltage difference U,,sin wt-U, becomes negative, and the reference diode starts to block the 
load circuit, But a reduction of the flux from the value + ; produces an emf which maintains the reference diode 
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in a conducting state; therefore a positive current i; flows through the load circuit. This continues until the flux de- 
creases to the original value @= 0, 


We see that when a symmetrical nonlinear resistance is present there is a negative voltage across the saturable 
reactor, wherease the supply voltage is positive. 


For t> t, the process is repeated, but in the opposite direction. The average value of the current i, over the 
entire cycle is equal to zero, 


When a signal(-Ic) is present, where | Wd. |=|Wd,| in the time interval t= t= ty (Fig. 7), the process occurs 
in the same manner as when there is no signal,except for the fact that the flux ©, is displaced downward by the 
amount 


i 
De = (6) 


where L, is the inductance of the saturable reactor from the control circuit side. 


At the instant t= t,,when U,psinwt = Up, the flux begins to increase when the original value ,. At the in- 
stant t= t, the flux retums to the same original value. 


Thus when there is a negative signal -I,,it is possible to have a change in flux during the positive half-cycle 
of flux density which is equal to 


im =|, — =| + > (1) 


but the actual flux change A@ when the signal is present remains the same as when there is no signal. The variations 
in the fluxes with and without a signal are,respectively, equal to (Fig. 5) 
A® = 0,— = 90,—0=9,, A® = 0,— = (8) 
The shape of the variation in the curve i; for t, = t = t, (Figs. 6 and 7) is identical with and without a signal. 
For tg=t St, the current i, and the flux vary according to equation (3). The initial conditions are 


t=, i, = 0, 
It should be noted that now the maximum possible flux change in the negative direction is equal to 


At the instant t = ts the core is saturated. Determining ty in the equation U,,sinwt, - Up * 0, we find the in- 
sant ts from the equation 


— i,r;) dt W,10* (—®,) + Lek. (10) 


The negative current i; for ts=t = t, is determined by the expression 


For ts =t= t, there is no voltage across the saturable reactor, the flux does not change, and = -4g. 


At the instant t = t, the current i, reaches the value (-1, + 1,), and the saturable reactor emerges from its satu- 
ration state, The instant t, is determined from the equation 
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U,,sinet —U, 


sin ot, —U, =r, (—1,+1.) 0. 


a 


; 

& 


Fig. 6. 


If we neglect the voltage drop it;, then the latter expression becomes 


(12) 


We shall derive the expression for the dc component ly of the current i, in the load circuit. 
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When the saturable reactor is unsaturated, the magnetizing current i, in the load circuit is small and can be 
neglected, 
Then 


We shall express the current I, in terms of the network parameters, 


The integral of the voltage across the saturable reactor with respect to time for tg t = t, is equal to the 
change 46_,, (Fig. 7): 


t, t, t 


since for tg t <t the saturable reactor is saturated, /* U; dt = 0, and the voltage U,,, sin wt—Up = iyr; as shown in 
Fig. 7. ts 
Neglecting the voltage drop i,r;, the Sin wt —Uy)dt (Fig. 7) which corresponds to presence of the 
signal is equal to the integral /* (U,,, sin wt-U,)dt = —W,10-*®, (Fig. 6) which corresponds to absence of the signal. 
Therefore for Fig. 7 


sin ot —U,) dt = — W,10-*O,, 


(cos wt, — cos wt.) U, (ts — W,10°*,. 


From (14) we obtain 


U 
—, + Lele ) = — (cos wt, — cos wts) — (ts — ts). (16) 
Substituting (15) and (16) into (13) we find the final expression for the dc component of the current iy in the 
load circuit: 


Lele) (17) 


f 0.4nW, pi0-* 
Ie), (17a) 


where 1 ,, is the length of the magnetic path in the core. 


From Eq. (17) it is evident that the dc component of current in the load circuit is directly proportional to the 
signal I, and the frequency of the supply voltage, and is inversely proportional to the resistance r, in the load circuit. 
For an increase (in absolute value) of the signal in the negative direction from zero to I,,the magnetizing current i; 
in the positive direction remains constant; the average value of the current i, in the negative direction increases. The 
average value of the current i in the negative direction reaches the maximum negative value I; neg, max for a con- 
trol current equal to ~I,. Substituting —1, = ~1, into Eq. (17), we obtain the maximum dc component 


Tomax = Tuneg.max = ( (18) 


(cos wt, — cos wt,) — U, ty). (13) 
or 
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Mode 2. This mode occurs when the magnetizing force produced by the 
dc control current is greater than the magnetizing force corresponding to the knee 
on the magnetization curve (Fig. 5) and less than the maximum possible mag- 
netizing force produced by the current i; in the load circuit (that is, 


Wels| <|We le |< |W, cf. Fig. 8). 


For a further increase in the modulus of the negative signal above | —1,| , 
the current i, in the negative direction will remain constant. Its average value, 
in accordance with Eq. (18), is equal to Ijneg, max- However, as we shall demon- 


strate below, for I, > 2 ies current i, in positive direction becomes greater 
than for 0 <| 1,| <| 


For mode 2 we shall denote the control current I'>. 


We shall wee how the current i, in the positive direction varies for |I',| > 
> 


saturable reactor is saturated and the current ty St<t', (Fig. 8) is determined from 


: Until the instant t';,corresponding to the knee of the magnetization curve.the 
! the expression 


U_, sinat —U 
i= — (19) 


When at the instant t = t'; the condition 


U,sinet,—U, Wel.—Wel, 


(20) 


is satisfied,the saturable reactor emerges from the saturation state, and the flux 
rises from the value ~,. When the flux varies within the limits +4, to —@,, the 
current i, in the load circuit varies within the limits (W,I', + Wel,)/W; to 
(Wel. If |I',], then we can assume that while the flux varies 
Fig. 8. the current maintains the constant value Wcl'./W,. The assumption that Ig | 


<< | I',| is equivalent to approximating the magnetization curve by three straight- 
line segments (Fig. 9): 


From (20) we determine the instant t';: 


= + are sin| (U, + Wi ~ + are sin[ U, (22) 
For tg St St’, the current i, has its sign maintained constant by the emf induced by the flux change. 
e At the instant t = t', the flux returns to its original value ~ @,, and the current i; goes to zero. 
Ri The instant t', is determined by equality of the time integrals of the voltage across the saturable reactor: 
The fe 
on \ (Um sin —U, — dt + | Um sin ot —U,— iyr;) dt = 0, 
i, 

3) ve ts 


Wd 


(Um sin ot —U,— dt=0, 
t 


& 
i 
\ 
| 
ph ot 
| 
| 
d | 
| 
) 
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U 
(U, ) — (cos wt, — cos = 0. (23) 


Computing t’, and t', from (22) and (23), it is possible to obtain the average value for the current i, in the 
positive direction: 


won — ty) +e (cos wf, — cos wt,) (24) 


| 
max 
Fig. 9, Fig. 10, Fig. 11. 


In this mode the dc component ly of the current is equal to the sum of the currents determined by Eq, (18) and 
(24): 


From Fig. 8d, Eq. (18),(24), (25), and from the discussion above it is evident that for |I,| <|I,|< Par “nis =| 
1 


the current 1; pos, (for t;=t <t',) increases in magnitude in the positive direction as I, increases, while the current 
Ig max from (18) remains constant. Therefore the dc component ly from Eq. (25) decreases as the control current 
increases. 
Mode 3, This mode occurs for a further increase in the 
iy signal |, to the value | I, | > a 
c 


; under 
1 

these conditions the magnetizing force produced by the de 
control current I, in the load circuit. Thus the saturable 


t reactor cannot become desaturated. The current i, in the 
load circuit does not contain a dc component. When the 
=0 signal increases, the current i, for t= t= t', tends toward 
V7, the shape shown in Fig. 8d for tg=t= &. 
1 Toe (led <Us\) Thus modes 2 and 3 are not working modes. For 


! mode 1 the de component lp is proportional to the signal 
I, (for 0 < < [Welgl). 


: The control characteristic computed from Eqs. (17) 
: and (25) is shown in Fig. 10. The characteristic which is 
v taken experimentally for Hg = 0.05 oerst, 1, = 0.5 ma, is 


Fig. 12, shown in Fig. 11. It is evident that its shape is close to the 
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theoretical shape. From Fig. 11 it is evident that the characteristic is practically a straight line for |1,| < |1,| , and 
this coincides with the theoretical conclusion (cf. Fig. 10). 


Figure 12 showed the oscillogram for the current i; in the load circuit. The nature of the curve verifies the 
theoretical conclusions (cf. Fig. 7d). 


The experimental investigation was performed on a core with the dimensions 20 x 28 x 7 that was made of 
the material "79MN." The core was made of a ribbon 0.04 mm thick. The ac winding had 400 turns of 0.2 mm 
diameter wire. The control winding had 1000 turns of the same wire. As a symmetrical nonlinear resistance we used 
the network shown in Fig. 3 and a carborundum resistor (¢ 20 x 5). The supply voltage had a frequency of 500 cps. 
SUMMARY 


When a nonlinear symmetrical resistance is connected in series with a single-core saturable reactor the current 


in the ac circuit contains a dc component which is proportional to the magnetizing current for |I,|< |I,|. This 
makes it possible to obtain controlled rectification. 


The derived expression (17) provides a preliminary quantitative relationship between the dc component of the 
current in the ac circuit and the network parameters. 


The author expresses his appreciation to M, A, Rozenblat for his valuable remarks and advice, 
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A TRANSISTOR AMPLIFIER OPERATIONING 
IN A THYRATRON MODET 


M. V. Ol'shvang 


{Moscow ) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 493-500, April 1961 

Original article submitted July 22, 1960 


The paper studies the operation of a transistor amplifier in a switching mode in which the tran~- 
sistors are maintained fully conducting up until the instant that the collector current passes through 
zero [1]. 


The characteristics of the amplifier are clarified when a magnetic amplifier is used to control the 
circuit. Experimental results are cited. 


A switching mode of operation provides the following basic advantages in a transistor amplifier in comparison 
to the proportional mode: 


1) A high coefficient of utilization for the transistor—up to 100 and greater; 


* Original Russian Pagination. See C. B, translation. 
t From a report made to a seminar in the Institute of Automation and Remote Control of the Academy of Sciences 
USSR on contactless magnetic elements, May 25, 1960. 
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2) a high efficiency—up to 95-98%, 
3) a high temperature stability. 


As a consequence of this the use of the switching mode is especially favorable for use in power output stages. 
In control systems such an amplifier often operates into a reversive RL dc or ac load (collector and induction motors, 


* 


ad 


if 


Fig. 1. 


However, a switching mode involves a substantial shortcoming where RL loads are concerned: when the “switch* 
opens (the transistor is cut off), the breaking of the current through the inductance causes overvoltages to develop across 
the transistor, and the energy stored in the inductance is dissipated in the transistor. 


For the case of a nonreversive dc RL load, this shortcoming is eliminated by shunting the load with a reversed 
rectifier. The use of uncontrolled reversed rectifiers for reversive amplifiers and for amplifiers with an ac output is 
impossible. 


The breaking of the current by the transistor can also be avoided when the load current periodically drops to 
zero due to supply from the ac line. Then the transistor operates in the following mode: before the arrival of the 
control signal it is cut off; at the instant the control signal arrives it is fully conducting and remains in this state 
as long as the load current flows through it, irrespective of the magnitude of the control signal; after the current passes 
through zero the transistor is again cut off. Thus the transistor made in this network is analogous to the operating 
mode of a thyratron [1]. 


The Principle of Operation for the Circuit 


The amplifier circuit is shown in Fig. 1a. The transistors T, and T, are connected in common-emitter fashion 
with a positive feedback signal through the current transformer TFB, The diodes D, and D, block the transistor from 
any voltage with the nonworking polarity. 


The TFB has the following windings: the working winding w,, through which the load current i, flows, the 
feedback windings Wo, = Wa = W» through which the transistor base currents ip; and ig, flow, the windings w, and w, 
which are subjected to the blocking i, and control i, signals, respectively, through the ballast resistors Rj, and R... 


electromagnetic clutches, etc.). 
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The wave shapes of the currents and voltages are shown in Fig. 2, We shall study the operation of the net- 
work by assuming that the feedback transformer is ideal. During the positive halfcycle the voltage is applied to the 
transistor T, and the diode Dz. Assume that at the instant voltage is applied to T, the relationship 


> (1) 
is satisfied. 
Here the direction of the magnetizing force (mmf) 4,W, is such that the transistor T, is cut off, and the current 


flows through the winding we; i. e., if (1) is taken into account we have 
; 
Under these conditions the magnitude of the voltage which blocks T, is determined from the formula 
dB 
— €or = = = WoS (3) 


Here B is the flux density in the core of the TFB; S is the cross section of the core; t@ = Ry» + Tiny, where 
Ry, is the resistance of the winding we, and rjn., is the input resistance of T,. 


Assume that at the instant © = 6» a control signal appears whose mmf is directed opposite to the mmf of the 
blocking signal; assume further that its magnitude satisfies the inequality 


i> (4) 


Here the resultant mmf will change sign, and the flux in the core will begin to decrease (dB/dt< 0), Asa 
result the conditions €9; > 0, eg < 0 will be satisfied; i. e., T, will shift into the active region, The process of a 
further shift of T, toward full conduction proceeds to completion in avalanche fashion (until T, is in the region of 
collector current saturation), 


After T, has become conducting, the inequality 
wp > ive > ©) 
is satisfied. 
From the instant @y we have 
Ww— + — 0. (6) 
From (5) and (6) we obtain the expression for the base current 


“LY! w (6a) 


= = 
The transistor is fully conducting if 


Where By is the current gain in a common-emitter circuit when U,o; = 0, and the specified value of collector 
Current is igo}. 
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cross 
is 
asses 
421 


From (6a) and (7) we obtain the condition governing full conduction for the transistor: 


= 
(8) 


kp = —, (9) 


where 6 is current gain of the transistor in a common-emitter circuit for the specified collector current and a col- 
lector voltage Up), < 9. 
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| After the transistor shifts into the active region it is encompassed by positive feedback with a coefficient 
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After the transistor has been shifted into the active region by the control signal,the inequality B> By is 
valid in the entire region where the current i; increases. If (8) is satisfied then kg, > 1, and thus at the instant @, 
the avalanche process of transistor conduction occurs and the transistor shifts into the region of collector current 
saturation. 


Condition (8) is sufficient for assuring avalanche conduction of the triode after the control signal satisfying (4) 
has been applied, provided that the process of reaching full conduction in the transistor occurs sufficiently rapidly in 
comparison to a cycle of the supply frequency. The necessary and sufficient condition for assuring that the process of 
reaching transistor conduction is an avalanche process when the transistor is connected in a common-~emitter circuit 
with feedback through a current transformer is cited in [2] for the case where the dynamic characteristics of the tran- 
sistor and the core are taken into account. 


The time required to reach full conduction in the circuit whose parameters are cited below amounted to 30-40 
sec for a resistive load, and 5-10 psec for an RL load. The latter value is explained by the fact that the diode be- 
comes fully conducting at the very beginning of the next rise of the load current from zero to the maximum, since 
the time constant for the RL load appreciably exceeds the time constant for the process of reaching full conduction. 
Therefore,when the transistor operates into an RL load the power dissipated during commutation can be neglected, 
and the power dissipated in the transistor is in practice determined solely by the power dissipated in the region of col- 
lector current saturation. 


At the instant 9 = con, (8 = © + Ocon,) the current ico}.3 (icol.2) becomes equal to zero, The quantity 9... 
is determined by the instant at which the transistor becomes fully conducting and by the quantity 9%, = arc tan 


The presence of the switching mode and the control of the transistors through the transformer permits easy 
realization of parallel and series connection of transistors; ae as pape Se corresponding number of feedback 
windings are added to the TFB. 


The Control of the Amplifier 


Changing the power to the load is achieved by changing the conduction angle 9%, of the transistor. Therefore, 
in order to control the amplifier it is possible to use various phase-shift networks. The most expedient network for 
this purpose is a magnetic amplifier (MA). 


Figure 1a shows a variant of a control circuit using an MA whose supplied directiy from the line voltage. The 
amplifier in Fig. 1 shall henceforth be called a magnetic~transistor amplifier (MTA); the MA is its input stage, and 
the transistor section plus the TFB is the output stage. 


While the MA is unsaturated (the excitation interval) inequality (1) is satisfied; after saturation inequality (4) 
is satisfied (i.e., at the instant the MA is saturated the transistor becomes fully conducting). 


On the basis of studying the graphs in Fig. 2a it is possible to come to a conclusion concerning the input-output 
response for the MTA. For operation into a resistive load it is the same as the response of the input stage (the MA), 
since the excitation interval of the MA coincides with the intervals during which the transistor is blocked. For oper- 
ation into an RL load,the interval during which the transistor is in the blocked state proves to be less than the ex- 
citation interval of the MA due to the lag of the load current relative to the supply voltage (cf. Fig. 2a); therefore 
the input-output response is steeper. Even for an MA saturation angle equal to the load angle ¢,, the load current 
teaches a maximum value and then remains constant. 


Evidently the most expedient form of control signal is a pulse. Under these conditions the ‘next pulse must 
arrive only after the load current has dropped to zero. On the other hand, the blocking signal can be applied only 
during the interval for which the transistor is in the blocked state. 


These conditions can be realized if the control and blocking circuits are connected in parallel with the output 
stage. Then the voltage across these circuits will be applied synchronously with the voltage applied to the transistors. 
Such a network is shown in Fig. 1b. The transistors T,, Ts and Tg, T4 are connected in series, The voltage across the 
magnetic amplifier MA is applied through a step-down and isolating transformer Tr. This lowers the weight and di- 
mensions of the MA and permits the control and blocking signals to be applied to the same winding of the TFB. 


The waveshape diagrams for i, and i, are shown in Fig. 2b for this case. The excitation interval for the MA 
and the interval during which the transistor is in the blocked state coincide for this network irrespective of the nature 
of the load, Therefore, the input-output response of the amplifier shown in Fig. 1b also repeats the MA response ir- 


respective of the type of load. Since the response of the MA depends only slightly on the load angle, it follows that 
the resultant response of the MTA will also have this property. 


Figure 3 shows the input-output responses for an MTA with the transistors connected in series according to the 
network shown in Fig. 1b. The supply voltage is U = 80 v (U,, = 115 v), f = 400 cps; the load impedance is Z; = 
30 ohm. 


We shall cite the basic parameters for the network. The TFB core is made of 80NKh material with 20 rings; 
Wy = 12, We = 70 (n = 11.7), and w,, = 500. The MA core is identical with that of the TFB; w,, = 2400, w.., = 3000 
(Rem = 300 ohm). The transistors T;~T, are germanium transistors of the type P4B (without additional cooling fins), 
D,~—D, are germanium diodes of the type D302 ; D,—D, are germanium diodes of the type D2V . The MA is con- 
nected through a step-down transformer with a ratio 6/1. 


The characteristic 1 (Fig. 3) was taken for the network shown in Fig. 1b and is practically indepndent of the load 
angle. It is easy to see that it repeats the conventional response for a self-saturating magnetic amplifier. The input- 
output response for the network shown in Fig. 1a (with the transistors connected in series as for Fig. 1b) coincides with 
the characteristic 1 for gy, = 0. The characteristic 2 was obtained for gy, = 40°, It has a pronounced limiting of the 
output current and a large linear segment; this verifies the reasoning cited above. 


amp 
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Fig. 3. 


The average value of the current at the output of the MA in the network shown in the Fig. 1b did not exceed 
2.5 ma. Thus the power gain for the output stage in Fig. 1b is 5 - 10°; in Fig. 1a it is 1- 10°. Such a high gain for 
the output stage is achieved due to the fact that the control signal must merely bring the transistor into the active 


region. 


The over-all power gain of the MTA over the linear segment of the characteristic is approximately 2 - 10”, 
The efficiency of the amplifier is 95-98%, 


The network was tested under conditions where the ambient temperature was varied through the range from 
-50 to +60°C. Under these conditions the instability of the MTA characteristic was determined solely by the tem- 
perature drift of the MA characteristic. 


We shall compare MTA networks with synchronous and independent connections. The advantage of the net- 
work in Fig. 1a is the fact that for operation into an RL load the MTA corrects the characteristic of the MA in the 
region of maximum output (where the gain of the MA arops sharply), since the region of maximum output from the 
MA is not used for controlling the output stage. However, in view of the necessity of synchronizing the voltage ac- 
ross the control circuits with the voltage across the output stage the use of Fig. la is difficult for polyphase circuits. 
The network in Fig. 1b is more flexible and can be used in all those cases where an MA is usually used. The other 
advantage of the network in Fig. 1b is its higher power gain for the output stage, 


Special Operational Feature of the TFB 


The magnetic state of the TFB core varies along a hysteresis loop as shown in Fig. 2. Here the core must not 
saturate; i.e., the condition 


AB < 2B,, (10) 


must be satisfied, where B, is the saturation flux density. 


From (3), (6a),and (8) we obtain 
R 
AB = 


R, = Ry, + Rp, 
as a result of integration; here 1; is the average value of the load current; f is the frequency of the supply source; R,, 
is the resistance of the winding wo, and R,,, is the average value of the resistance rp over a half cycle. 


For the case of TFB saturation before the instant 6,, premature blocking of the transistor occurs due to dis- 
ruption of the feedback. 


An appreciable difference in the magnitudes of the resistances Rg, and Re in the feedback circuits leads to 
saturation of the TFB during that half cycle when the transistor operates with a larger distance Rg [this follows from 
Eq, (12a)]. An insignificant deviation in the values of Rg; and Re will have no effect due to the fact that inequality 
(8) is satisfied with some margin, 


The output stage of the MTA is practically inertialess. 


For a nonuniform and asymmetrical arrangement of the windings on the core of the TFB,the normal mode of 
operation of the MTA is violated; this is a consequence of the harmful effect of the leakage inductances. The leak- 
age inductance of the feedback winding causes a certain lag in the current ig with respect to the current i, in the 
primary winding during the interval for which the output stage is fully conducting, and with respect to , during the 
interval for which the transistor is in the blocked state. At the instant the primary current passes through zero, the 
current ig » 0; this current then becomes the magnetizing current, the flux begins to change in the opposite direction, 
and a positive voltage is applied to the emitter junction of the second transistor at that moment when it must be 
blocked. 


In view of this,the leakage inductance of the feedback winding relative to the windings w, leads to a growth 
of the no-load current, and the leakage inductance of the feedback winding relative to the warking winding may cause 
the appearance of step variations on the input-output response of the MTA, 


These phenomena are absent when the distribution of the windings on the toroid of TFB is uniform. In the net- 
work shown in Fig. la the no-load current did not exceed 10 ma for the parameters indicated above. 


Another method for eliminating these phenomena is to shunt one of the TFB windings witha capacitor whose 
magnitude (normalized to the winding wg) is 0.2-2 yf. However, a capacitor of several microfarads leads to an ap~ 
preciable retardation of the process involved in reaching full conduction and thus causes an increase in the losses. 


The effect of the leakage inductance is the smaller, the greater the magnetic permeability of the core ma~ 
terial. As tests have shown, it is possible to use either permalloy cores or cores made of steel £31C for the TFB. 


The feedback transformer can be treated as ideal if we choose the value of blocking current in such a way that 
during the greater part of the cycle the inequality i, » iy is satisfied, where i,, is the magnetizing current nor- 
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malized to the winding w,. In practice it is sufficient for the current }, to exceed the magnitude of the magnetiz- 
ing current by a factor of 5 to 10. 


SUMMARY 


In the proposed magnetic-transistor amplifier the transistor operates in a mode which is analogous to a thyratron 
mode: when the control signal is applied it becomes conducting via an avalanche process and remains fully con- 


ducting while current flows through it. This is achieved by introducing positive feedback through a current transformer, 


As a result the amplifier acquires the following advantages: 1) the possibility of operating into an RL load; 2)a 
small time required to reach full conduction (this time is'practically independent of the wave front of the control sig- 
nal, and thus a high use coefficient forthe transistors and high stability are achieved); 3) a large gain for the output 
stage and a low inertia. 


Based on the network analyzed above, we designed and successfully tested amplifiers with an output working 
into a reversive induction motor (two-phase and three-phase), A somewhat modified version of the network can be 
used for feeding the armature of the dc reversive motor. 


The author expresses his deep appreciation to R, A. Lipman who made a number of valuable suggestions when 
this paper was edited. 
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ON THE THEORY OF SELECTIVE RECTIFICATION OF EVEN 
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The paper studies one of the networks which is used for selective rectification of even voltage 
harmonics using symmetrical nonlinear electrical resistors. The application of this network in 
the output circuit of a "ferroprobe" (or a magnetic amplifier of the even harmonic type) is 
studied. A simple formula is derived which makes it possible to evaluate the sensitivity of 
such a “ferroprobe.” Experimental data are cited which verify the theoretical conclusions. 


In a number of papers [1-9] a study is made of selective rectifiers based on symmetrical nonlinear electrical 
resistors (SNER). In [7, 8] a study is made of the operation of SNER that are used for selective rectification of the 
even harmonics in a voltage obtained from the output (indicating) winding of a “ferroprobe." It is demonstrated that 
the presence of a sharp knee near a certain point on the volt-ampere characteristic of the SNER is of great significance 
in increasing the sensitivity of the network. In this regard the SNER network treated by R. Ya. Berkman [9] is very 

favorable. It consists of a Gretz bridge consisting of four rectifier elements connected to a common blocking voltage. 


* Original Russian Pagination. See C. B. translation 


at 
nce 


ge. 


The volt-ampere characteristics for such SNER is shown schematically in Fig. 1. Near the points U;, and —U}, whose 
position is regulated by the magnitude of the blocking voltage, there is the required sharp knee. 


R. Ya, Berkman uses an expansion of the measured voltage into a Fourier series; as a result the final formula 

contains the sum of the series and essentially makes it possible to draw merely the qualitative conclusion that the 
network under study achieves selective rectification of the even harmonics in the 
measured voltage. As far as quantitative computations are concerned, this for- 

l mula can be used for simple computation only in certain particular cases (for 

. example, the case where one of the even harmonics is isolated). We have derived 

a sufficiently exact and at the same time simple solution for the SNER network 
cited above which has been justified by experimental results. 


| 
Theory 


The network under study is shown in Fig. 2. We assume that: a) the capaci- 
Fig. 1. tances are sufficiently great so that the ac component of the potential difference 
across their plates can be neglected in the steady-state mode; 


b) the characteristics of the rectifier elements are identical, and in order to achieve additional simplification 
of the computations we assume that the reverse resistances of the rectifier elements are infinitely large; 


c) the supply voltage u,. does not contain even harmonics; 


d) the amplitude e,, of the measured voltage e is considerably less than the amplitude u,;, of the commutating 
voltage u,: 


(1) 


is formed. For |U,| = U, the bridge is blocked and i'y= i*g = i= 0, For |Ug| > Uj, (i.e., in the “transmission in- 
tervals" r' and 7”; cf. Fig. 3) the bridge is unblocked and 


” 


t=T 


= ¢. (2) 
If the measured voltage is e = 0, then in view of the symmetry of the network we have: 
i=0, ip = ig = io, 


” 
Up =U 


on 
ner, 
a 
e<u, =u || 
(ny and ng,respectively, denote the number of turns on windings I and I). « 
th 
‘ 
\/ 
Fig. 2. Fig. 3, . 
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For e » 0 we have 


(3) 


where I denotes the dc component of the current i flowing through the load resistance R;. From (1) it follows that 
i « ig. Choosing the conditional directions of the currents i'y, i"9,and i (cf. Fig. 2) for the time segments r‘ and r* 
when the bridge is unblocked, we shall obtain the following Kirchhoff equations for the symmetrical loops: 
“¢ + = =) + (+ 5) (4) 


FB + = 5) + (io — 5) 6) 


Here Uy = R,I is the potential difference across the load resistance, U; (ig + i/2) is the potential difference 
across the rectifying element, R,,;, Rhg), are ballast resistors, where R},,) includes one-half the resistance of the 
winding Il; the upper sign in front of U,/2 corresponds to the interval r', and the lower sign corresponds to the in- 
terval r". Since i « ig, it is possible to write approximately 


ty, (ig = Up ig) = tp, 6) 


where R, q (ig) is the differential resistance of the rectifying element. Adding (4) and (5) and taking (6) into ac- 
count, we obtain 


Uy = (io) + + bo (7) 


This equation determines the current ip in the intervals r' and r". Over the remaining portion of the cycle, 
ig = 0. From (7) it follows that a) the current ig has a waveshape in the form of peaks (cf. Fig. 3), and b) if uz 
does not contain even harmonics (i. e., uw(t) = —u (¢ +37), then ip also does not contain even harmonics; 


Taking (8) and (2) into account together with the rectifying effect of the bridge, we have 


(9) 


We shall integrate (7) over the limits of the interval r' while assuming u, = u,,,, sinwt. Taking into account 
the fact that Up, = Ug, cos w7/2 and that according to the generalized theorem of the mean we have 


ts ts 


dt = \ Ry (ia) dt = Fy, \ ig at 


where R, is a certain “average” value for the resistance of the rectifying element, and taking (9) into account, we 
obtain 


in view of the same symmetry. 
| It is evident that 
I 
AU, = 
2¢ 
Uy = 2Rp | indt. = 
t 
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This equation can easily be solved graphically. However, for practical computations it can be further sim- 
plified. For wt < 100° we have 


with an error that does not exceed 30%, 


Substituting this into (7) we obtain 


~0.15 (11) 


Formula (11) yields a value for the ratio r/T which is somewhat too high in comparison to the value obtained 
from (10). On the other hand, if we assume that the reverse resistance of the rectifying elements are infinitely large, 


this very fact causes a certain decrease in the theoretical value of r/T. These two errors thus often compensate each 
other. 


From (4), (5),and (6) we obtain the equation for i if we take into account the fact that AU}, =(R,/2)-1 and 
Up = RI: 


this €quation applies to the intervals r' and r". Over the remaining portion of the cycle i= 0. From (12) we obtain 


4 (eae 4. (\ +) idt) + (“bal + + R,) as) 
From the generalized theorem of the mean we have 


¢ ans Re Ted. 
Ry. dic) idt Rr dic) idt) = )= (14) 


t, 1 


where R, g is a certain “average” differential resistance of the rectifying element. We shall study the left-hand side 


t 
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of Eq.(13). The voltage e, generally speaking, consists of two components €ey and eoqq, where e,, denotes the 
sum of the even harmonics of the voltage e and has the property 


fey, (t= 5) (15) 


the quantity eogq tepresents the sum of the odd harmonics of the voltage e and has the property 


Codd (=~ (t + 5). 


From (16) and (2) it follows that 


4 
Coaadt = 0, (17) 
ty 
Using the substitution 
t 
oy. tt (18) 


and taking (2), (14), (15),and (17) into account, we obtain 
BE, 
Rea , Ry (19) 
(3 + 


2 


from (13). 


Thus the dc component of the current flowing through the load 
resistance R; proved to be proportional to the average value of the sum 
of the even harmonics in the measured voltage over the transmission 
interval +. 


The SNER network under study can be used, in particular, for 
rectifying the sum of the even harmonics in an emf excited in the out- 
put winding of a “ferroprobe” (or an even-harmonic magnetic am- 
plifier). In that case we have the following well-known relationship 
[8] for an operating mode close to no-load: 


d 
2nS @0) 


(In the MKSA system). Here Hy is the “controlling” field, y, is the 
differential permeability of the core, n is the number of turns on the 
indicator winding, S is the cross section area of the core. Substituting (20 into (18), we obtain 


= — tar) Ho. (21) 


Here yg , and yq., are the values of the differential permeability at the instants t; and &,respectively. In 
order to obtain the maximum proportionality coefficient between I and Hg (i.e., the maximum sentivity of the net- 
work) it is necessary for the “phase shift” between e,, and the commutating voltage u, = ngu~/n, to be in accor- 
dance with that shown in Fig. 3 (i.e., the instants t, and ty coincide), Under these conditions the magnitude of the 
interval must be confined to the limits delineated by and Tz. In that case 9, = Ud, 

_ therefore, 


| 


E,= — He (22) 


The indicated optimal phase shift between e., and u, can be achieved practically by using a phase shifter or 
simply by choosing the appropriate resistance for the excitation circuit. 


From (19) it follows, in particular, that the sensitivity of the network depends substantially on the ratio r/T 
which can be computed approximately from (11). The quantity R,. which appears in this formula can be evaluated 
from the characteristic of the rectifier element only in a very approximate fashion [it can often be assumed that 
R,, = Ry. (ig/2)]. Notwithstanding this, computations according to formula (11) will prove to be sufficiently exact 
since the ballast resistors Rj, are introduced into the network to compensate the effect of instability in the rectifier 
element characteristics on the operation of the network; they must thus be appreciably greater than R, . However, 
Rpg, Must not be too great in order to avoid lowering the network sensitivity. If we substitute ratio 1/T found 
from (11) into (19) and find the optimal value of Rj, from the condition governing maximum sensitivity (assuming 
that E, is independent of Rp), we obtain 


Ry opt= (23) 


However, it is necessary to take into account the fact that this relationship is valid only within restricted 
limits. In fact, Rj, must not be too great or too small, since in either of those cases the interval r may,in accor- 
dance with (11),prove to be greater than Tz orlessthan ry, From (11) and (19) it also follows that an increase in 
Rbal.a» Rbal.2-and Ry leads only to a lowering of the sensitivity. 


If the forward resistance of the rectifying element on the working segment of the characteristic and the value 
of the ballast resistance Rj.) are small compared to the load resistance, then the terms (1/2)R, qg, and (1/2)Rpai.4 
in the denominator of formula (19) can be neglected. Substituting (22) into (19) we have 


for this case. 


For approximate computations that give a rough value for the sensitivity of "ferroprobes" in which the net- 
work under study is used for rectifying the sum of the even harmonics of the emfs excited in the output winding, we 
obtain the following results from (24) if we assume that Ry), = 0 and take (23) and (11) into account: 


(25) 
Here f is the excitation frequency. 


For a magnetic amplifier with the same output circuit,the current gain can be computed approximately from 
the formula 


I T (26) 


V + 


where I is the control current, Mg is the number of turns on the control winding, and ! is the length of the magnetic 
path in the core. 


In concluding this section we shall note that, as subsequent detailed analysis shows, the proportionality between 
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I and Hg is still maintained for the case where the reverse resistances of the rectifying elements are finite (for the 
condition that their volt-ampere characteristics are identical), 


Experimental Data 


An experimental verification of the derived theoretical results was obtained using magnetic amplifiers of the 
even harmonic type with an excitation frequency of 500 cps. We used germanium diodes of the type DGTs-24 as 
the rectifying element. The capacitors had a value of 20 yf. 


A. The table cites a comparison of the ratio r/T for various Rj, computed according to formula (11) and 
measured directly by means of an oscilloscope. 


The agreement is satisfactory. 

B, It is best to perform the direct verification of formula (24) using magnetic amplifiers with ferrite cores for 
which the quantity pq max, d0es not depend on the excitation mode and can therefore be measured easily. It is 
necessary to take into account the fact that formula (24) was derived on the basis of assuming no-load operation. 
Therefore in performing the measurements the ballast resistor Rj,,; .is introduced into the network. Figure 4 shows 
the dependence of tq max, Computed according to formula (24) from the experimental data on Rj,,; 2. It is evident 
that the curve uq.max, = Hd.max.®aj.2) 28ymptotically approaches the value corresponding to no-load operation, 
For various R, the no-load values of ug max, obtained in this manner fluctuate within the limits of from 4.8 - 10~ 
to 5.1 - 10 (web/amp - m). Direct measurement by the simplest method yields wg max. = (4.7 + 0.3)107 
web/amp m. 


C. Further experiments were performed using toroidal cores made of 79NM permalloy. The mean diameter 
of the toroids was 3.0 cm, and the cross section area was 0.25 cm*, The number of turns on the excitation winding 
was 63. Formula (24) can be used for this case in order to perform an approximate computation of 44 4x. from the 
measured ratios I/Hy and r/T. The results of such a computation for various Rj, are cited in the table. The "valleys" 
in the function ug max, = Md.max.(8b) at the edges afe explained by the fact that for Rj, < 1 kohm the trans- 
mission interval proves to be greater than rT 2,whereas for R, > 4 kohm it is less than 1. 


D. In accordance with (21) the function H,/I= f(R,) must be a straight line. This is verified by experimental 
results (cf. Fig. 5 which shows a typical curve for Hg/I as a function of R; from among the series of curves which we 
obtained). According to (24) the segment Ry shown in Fig. 5 must equal R,,/2 for Rg), = 0. Figure 6 shows the ex~- 
perimentally obtained functions Ry = Ro(n) for various R;,. For n -> 0 (this corresponds approximately to no-load oper- 
ation) we have Ry = (1/2)R,. 


E. From Fig. 7, which shows the experimentally obtained curves for the sensitivity I/Hy as a function of R, 
for various load resistances R;, it is evident that the inaximum sensitivity corresponds to a value of Rj, that equals 
R,; i.e., relationship (23) is also verified experimentally. 


SUMMARY 


1, For the indicated network the dc component of the current flowing through the load resistance is propor- 
tional to the average value of the sum of the even harmonics in the measured voltage over the transmission interval. 
It also depends substantially both on the load resistance R;, and on the resistance R}, across which the blocking voltage 
is produced. 
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| mu (0.22 | 0.52 | 1.02 | 2.0 | 4.00 7.00 
| | 0.24 | 0.20 | 0.16 | 0.13 
0.31 | 0.26 | 0.2 0.49 | 0.10 | os 
0.020 | 0.000 | 0.053 0.088 | 0.053 | 0.049 


al. 
age 


2. An approximate computation of the sensitivity for a "ferroprobe" in which the network under study is used 
for selective rectification of the sum of the even harmonics in the emf excited in the output winding can be per- 
formed according to formula (25). The approximate computation of the current gains for a magnetic amplifier with 
the same output circuit can be performed according to formula (26). 
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A SIMPLIFIED METHOD FOR DETERMINING STABILITY 
SEGMENTS 


Vv. E. Bogin 


(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No, 4, 
pp. 509-510, April, 1961 

Original article submitted August 6, 1960 


As we know, in investigating the stability of linear automatic control systems the method of D~-subdivision is 
used for determining the set of parameter values for which the system remains stable. 


In investigating systems by this method it is necessary to plot a D-subdivision curve which is associated with 
rather cumbersome computations. In practice, it is often necessary to determine a stability segment (i.e., a seg~- 
ment of the positive abscissa semiaxis which located within the stability region. 


Stability segments can be determined using a simplified method that is a corollary of the D-subdivision method 
involving one parameter. 


As usual, it is assumed that the coefficients of the characteristic equation for the closed-loop system depends 


linearly on the investigated parameter w; i.e., G(p) = Q(p) + wR(p), where G(p) is the left-hand side of the character- 
istic equation for the closed-loop system. 


As we know, the equation for the D-subdivision curve is given by 


Q (ja) 


' Original Russian Pagination. See C. B. translation. 


After separating the real and imaginary parts we obtain 
w (@) = A; (@) + (@). 
An example of the D-subdivision curve is shown in the figure. 


From the figure we derive the following: 


1) the boundaries of the stability segments are determined by the values of A;(w),), where w, are the real 
positive roots of the equation A,(w) = 0; 


2) by computing A,(w,, + Aw) and Ag(w, + Aw) it is possible to determine the direction of the hatching (i.e., 
it is possible to isolate the segments which are “candidates” for the stability segments). 


Thus the problem of determining stability segments reduces to the following operations: 
1) determining the real positive roots of the equation A2(w) = 0; 

2) computing the values A,(w,); 

3) computing the values A,(w, + Aw), Ao(w, + Aw); 


4) verifying the system for stability for a value of the parameter w taken from the segment which is a can- 
didate for the stability segment. 


The reduction of the computations in the method cited above is obvious, since here there is no need to deter- 
mine the large number of values w = w(w) which must be found for accurate plotting of a D-subdivision curve. 


The number of computations required for an approximate determination of the roots w,, is small, since such an 
approximation can be rather rough. 


There exist effective methods for determining the roots of the equations, For an approximate determination of 
the real roots it is convenient to use the method based on the Viet formula and the Bernoulli method. The results 
obtained by these methods can be refined using linear interpolation, the iteration method, the Newton method etc. [1], 


For equations which are not of high degree (up to the sixth) it is possible to use the Kryzhanovskii iteration 
variable w is substituted for according to the formula 


o=§ +5, 


here b is the upper boundary for the positive roots and 9 is a new variable. 


The computation of the values for the polynomials and their derivatives is conveniently performed according 
to the Horner scheme. 


In conclusion it should be noted that the method cited above is not expedient if the curve for D-subdivision 
with respect to one parameter is used for judging the system performance etc. (cf., for example, [3]). 
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THE DESIGN OF AUTOMATIC CONTROL SYSTEMS (ACS) 


Yu. I, Borodin and V. N. Plotnikov 


(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 511-523, April,1961 

Original article submitted October 31, 1960 


The paper develops a method for designing automatic control systems without lag for the case of 
first and second order objects when four types of controllers are used, The method is generalized 


to include systems with lag. Nomograms are given for choosing the parameters of the desired 
transfer function and for determining the type of controller. 


A number of papers [2-6, 9-13] have been devoted to the adjustment of controllers and the synthesis of cor- 
rective networks in ACS. 


In this paper we present a method for choosing the controller parameters on the basis of the primary perfor- 
mance indices in the case of first and second order objects (for the case of no pure lag and for the case where lag is 
present), and for four types of controllers: a) with respect to deviation (D); b) with respect to the integral (I); c) 
with respect to both deviation and integral (D1); d) with respect to deviation, derivative and integral (DDI). The 


method under study is a concretization of the method for synthesizing corrective networks in ACS which was based 
on choosing the desired logarithmic amplitude response (LAR) [4]. 


We shall study the block diagram of the ACS (Fig. 1a). It consists of a controlled object with the transfer 
function KgWo and a controller with a transfer function K,W,,(s) whose type and parameters are unknown, 


The transfer function for many objects can be represented by one of the expressions cited in Table 1; the tran- 
sfer functions of the controller types indicated above are.cited in Table 2. 


As the primary dynamic indices characterizing the automatic stabilization system (Fig. 1b) we shall analyze: 
1. The magnitude of the maximum deviation x,,,,, caused by a unit stop perturbation. 


2. The transient response time T caused by the same perturbation. (The minium time after which the devia- 
tion of the controlled quantity differs from its steady~state value by an amount less than or equal to A. For static 
systems T is determined by the formula x, ~A< |x(T)| < xg, + 4, and for astatic systems by the formula |x(T)|<A. 
The quantity A is determined by the engineering conditions, Here A is equal to 5% of the urit step perturbation). 


3. The time t, during which the deviation reaches its maximum value. 


4, The relative dynamic deviation of the controlled quantity, which characterizes the shape of the transient 


response; it is given by 1, = Xpy/X;y (for a static system—the curve 1) or 1). = Xg2/Xmax (for an astatic system — 
curve 2), 


5. The static system error xg, which is characterized by the deviation of the controlled quantity from the 
nominal value after the transient response caused by a unit step perturbation has been completed. 


First we shall study the method of choosing the type of controller and its parameters for objects that have a 


tansfer function of the type (1) and (2); then we shall generalize the method to include systems whose object has pure 
lag. 


‘TABLE 1 


K, 
| (1) 


DPI 


Kg + + 1) (8) 
& 


1st order 2nd order 


I, DI, DDI I, DI, DDI 


(10) 1) ( 


Fig. 1. Block diagram for the ACS (a), and the performance indices (b). g(t) is the stan- 
dard input signal [g(t) = const], f (t) is the perturbation, ¢(t) is the controlled quantity, 
x(t) is the deviation of the controlled quantity, y(t) is the control signal. 


Nomograms for choosing the desired transfer function of the closed-loop system. These nomograms establish 
the relationship between the parameters of the transient response and the transfer function for the closed-loop ACS, 
The transfer functions Y(s) for the ACS corresponding to the objects and controllers under study are shown in Table 3. 


The nomograms were plotted using formulas which determine the relationship between the primary performance 
indices T, Xmax, t, and , and the parameters of the transfer function for the closed-loop system K, Tyand {, 
(Table 4), The formulas were derived using the inverse Laplace transform. 


The nomograms are plotted for the transfer function (10) of the closed-loop ACS in the coordinates 20 log K vs. 
£23 te /Ty; 1 vs. €3 (Fig. 2) for the transfer functions (11) and (12) they are plotted in the coordinates 20 log K/T, 
V8. €33 t./T2, 1 v8. Cy (Figs. 3 and 4), The curves in Fig. 4 are plotted only for the case T, = T,; however, they 
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can be used for an approximate determination of the performance indices for T,/T, #1. For fixed values T,/T, #1 
the curves can be plotted analogously. The allowed value of the deviation A for the controlled quantity when the 
formulas in Table 4 are derived is assumed equal to 5% If Ahas a different value, then the axis of the abscissas 
should be shifted by the amount 20 log (4/0.05). 


20 log db 09 8 8 9 
5 — + 
LAE 
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FV 
— Nes 
GS 0607 08 09 0 1S 
NZ 
to, 
@ 02 03 04 06 07 08 09 12 18 14 15 


Fig. 4. Nomogram for choosing the desired transfer function of a closed- 
loop third-order ACS (astatic controller), 


of the closed-loop ACS, 


Example 1 (Fig. 3). The following requirements are imposed on the transient response: x, = 0; X 3, = 0.4; 


T = 25 sec. 


We shall specify T/T, = 5. Then T,=5 sec. At the intersection of the curves T/T, = 5 and x,,,, = 0.4 we 
find the point A with the coordinates 20 log K/T,; = -2.5 db or K/T, = 0.75; ¢,= 0.6. From this we have K= 
= 5 + 0.75 = 3.75, Under these conditions Yq(s) = 3.758(5°s' + 2 - 0.6 - 5s +1), At the same time, t,=T,"1.2=6 


(point B) and 1 = xg2/Xmax = 0.1 (point C). 


The nomograms for determining the type of controller in choosing the desired transfer function for the closed- 
loop system when the transfer function for the object is taken into account. We shall study a method of plotting 
homgrams for a first order object and controller of the type DDI (a). Astatic controllers of the type DI (7) and I (6) 
shall be treated throughout as particular cases of (8) when the coefficients of the derivative and the deviation (or just 


We shall use an example to demonstrate how the nomgrams are used for selecting the desired transfer functions 
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of the derivative) are equal to zero. Setting the transfer function for the closed-loop system equal to the desired 
transfer function (11), we have: 


Ke 


In order to satisfy the latter equation it is necessary that the condition 


i T. 
be satisfied. 
Solving these equations simultaneously, we obtain: 
K K 


If r* = 0, the coefficient of the derivative is equal to zero; then we obtain the condition for a proportional- 
integral controller of the type DI: 
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5) 


6) 


18) 


For the case of a type I controller the coefficient of the deviation must also be equal to zero: 


% 0. (19) 
Writing Eqs. (18) and (19) in the form 
20 log 20l0g (20) 
and 
20 log 720108 21. (21) 


we plot the curves in the coordinates 20 log =——vs. ¢ (ig. 5) on the same scale as that used for the nomogram in 


Fig. 3 for various values of T,/Tg (the thin solid line) and Ky= 1. The think solid line corresponds to Eq, (21), If 

Ky * 1, then it is necessary to shift the axis of abscissas by 20 log Ky. In choosing the desired transfer function for the 
closed-loop system it is necessary to make use of the nomogram in Fig. 5 in conjunction with the nomogram in Fig. 
3, For this purpose it is necessary to superimpose one nomogram on the other in such a way that the distance between 
the axes of abscissas equals 20 log Kg, For convenience in using the nomograms it is desirable to plot one of them 
(for example, the nomogram used to choose the type of controller) on transparent paper. 


20 db 20105 
‘| 
5 SALTS 
0 
-5 0 
4-- 
lo 
2 138, 08 
Fig. 5. Nomogram for determining the type Fig. 6. Nomogram for determining the type 
of controller for a first order object without of controller for a second order object without 
lag. lag. 


The parameters of the desired transfer function in the case of a type DI controller should be determined from 
the thin solid curve; in the case of an I controller they are determined from the intersection of the thick solid line . 
with thin solid line. 


For a second order object and a type DDI controller (8) we set Y(s) and Yq_ (s) equal to one another and find: 
K K K 
= + = + Ko = Tilt + (22) 


The first condition (22) is necessary for the case where all three types of controllers are used (this condition 
Corresponds to coincidence of the high-frequency asymptotes of the desired logarithmic amplitude response (LAR) 
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and the LAR of the object; the second condition (r* = 0) is necessary for a controller of the type DI, and the third 

condition (2¢,T, = 0) is for a controller of the type Il. Solving Eq, (22) relative toK/T, and taking the logarithm, 

we obtain formulas from which curves in the coordinates 20 log ( K/T,K) vs. ¢, (Fig. 6) for various values of T,/T, 
(the broken-line curves) and {4 (the thin solid curve); the 
thick solid line corresponds to the latter equation, The 
nomograms were plotted for Ky = 1. 
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parameters of the desired transfer function for the case of a 


DDI controller are determined from the broken-line curves 
Ol IN (T/T); the parameters for the desired transfer function in the 

YET case of a DI controller are determined at the intersection of 
ins 


Pl 2 the broken-line curves(T,/T» = const) and the thin solid 
rT TT curves ({ 9 = const); the parameters of the desired transfer 
mn a function for the case of an I controller are determined at the 
+ UJ intersection of three curves; i.e., under these conditions the 
4 a dos latter equation in the system (22) is taken into account (the 
+ thick solid curve). 


After the desired transfer function has been chosen, the 
parameters for controllers of the type DDI, DI,and I can be 
“mM - determined from formulas (22) in the case of a first order 
object. For all of the indicated types of controllers the gain 


isKe= (28). 
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Fig. 7. An example showing the selection of the 
desired transfer function for the ACS while taking 
the type of controller into account. The curves 

of the nomogram in Fig. 4 are denoted by broken 
lines. 


We shall use an example to clarify how to use the 
nomograms to choose the desired transfer function while 
taking the type of controller into account. 


Example 2, The transfer function for the object is 


1.78 
KoWo(s) = 10s + 1° 


It is required to clarify the possibility of using various types of controllers, The transient response is subject to 
the following requirements: x,,,, =0.9, T= 80 sec. We shall superimpose the nomogram in Fig. 6 onto the nomo- 
gram in Fig. 4 (Fig. 7) by raising the axis of abscissas by an amount equal to 20 log Ky = 5 db. 


The parameters of the desired transfer function in the case of a DDI are determined at the intersection of curve 

= 0.9 with the straight line T/T, = const. If T;/T, = 1.5 (point A,), then T/T, = 9.5 and T = 143 sec. If Ty/Ty* 
= 1.25 (point Ag), then T = 81 sec. If T/T, = 1.0 (point As), then T/T, = 11 and T = 110 sec. Thus for a specified 
value Of Xmgx the transient response with the minimal duration will be that corresponding to the desired transfer 
function 


35s 
The transient response time for the case where a controller of the type DDI is used can be reduced if xmax is 
less. Thus, for example, if x, = 0.4, then T = 45 sec. 


tersection of the curves x,,,,, = 9.9 and {9 = 0.9 (point B) and the straight line T;/T» = 1.15 passing through this point. 
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in- 


Thus we find T, = 11.5 sec, T/T, = 7.5 and T = 86 sec. The desired transfer function is written as 


25.88 


Ya@= F + $1) 


whereas K/T; = 2.24; ¢4 = 0.5. 


The transient response time here is greater than the required value, and it is not possible to reduce it 
for the type of controller specified. From Fig. 7 it follows that a controller of the type I cannot be used in this case, 


since this would lead to a large deviation of the controlled quantity (the curve 2¢.- T, = 0 is above the curve Xm, = 
= 1,2). 


wT, 
40° 


Fig. 8, Families of logarithmic amplitude responses for a closed-loop 
ACS with lag. a) First order object; b) second order object, {» = 1.0, 


pure lag (3) and (4). The transfer function for an object with lag is written as 


(s) = (s) (24) 


where KgW,(s) is the transfer function for the object without lag;e"7$ is the transfer function for the lag section, 
The transfer function for the closed-loop ACS in that case is 
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and includes as multiplier the same lag-section transfer function as does the object. Therefore the desired transfer 
function can be written conditionally as 


Yq, (s)e"™, (26) 


KW, (s) 


The transient response in a system that has the transfer function Y‘g(s) differs from the transient response in a 
system that has a transfer function of the form (25) by a time displacement in the amount r. 


For systems whose transfer functions are determined by expressions (25) and (27) the characteristic equations are 
identical. 


The equation has an innumerable set of poles. Among these poles it is possible to isolate the principal pair of 
complex poles which determine the nature of the transient response [9, 11]. 


First we shall study the method of choosing the parameters for the static controller (5). The characteristic equa- 
tion for the system in the case of the object (3) is given by 


Tos +1 + KoK,e~* =0, (28) 
and in the case of the object (4) by 


T2s? + 2oT os + 1 =0. (29) 


Figure 8 shows families of the LAR for a closed-loop system with a relative lag T = 1/T, = 0.3 and first order 
(a) and second order (b) objects (¢»= 1.0) when the principal pair of poles is taken into account. The LAR L,, | Yl, 
Lnl¥el » Lin! Ys]. Ly,!¥4] are plotted for the transfer function (27) as the values KgK,, = 0.58; 1.0; 2.16; 4.6. 


Figure 8 shows the lines for point M (the point M corresponds to the ordinate of the asymptotic LAR for w, = 
= 1/T) and the point N (the point N corresponds to the ordinate of the real LAR for w,; = 1/T,). The families of 
LAR which are cited make it possible to select the controller parameters, We shall clarify this using an example. 


Example 3 (cf. Fig. 8b), The transfer function for the object is KgWo(s) = 1.2e"/c2"s*+ 2 - 128+ 1) It is 
required to determine K, such that the performance indices become X,,= 0.25; Xgy $0.4; T $60 sec. Then 

20 log (K/Kg) =-13.6 (the low-frequency asymptote of the desired LAR). The point My at which this straight line 
interesects the curve M determines the conjugate frequency of the desired LAR.T, = T/2 = 6, The segment M,N, = 
20 log 2¢4 = -8 db, whence ¢4 = 0.20, K, = Kg — K/KgK = 1.2 — 0.25/1,2 + 0.25 = 3,17. From the nomogram in Fig. 


2a we find that the chosen desired transfer function corresponds to a transient response with the following performance 
indices: x,,,, = 0.4; T/T, = 8; T = 48 + 3.6 = 51.6 sec. 


Families of LAR similar to those shown in Fig. 8 can be plotted for systems with controllers of the type I, DI, 
and DDI; they can then be used to choose parameters that insure the required performance indices. However, it is 


simpler to generalize the method treated above for ACS with no lag and astatic controllers to include systems with 
lag. 


For the case of a first order object with lag (3) and a controller of the type DDI (8) the transfer function for the 
closed-loop ACS is determined by the expression 


Y,(s) =Y,(s)e-™. (30) 


c 


Y.(s)= 


| 


The desired transfer function Y' ,(s) [cf. (26), (27)] will be written as 


Ks 
T+ 


Yq (32) 


Requiring that the roots of the characteristic equation (31) must be equal to the desired roots, we determined 
the roots of the characteristic equation (32) by substituting them successively into (31); solving the resulting equations 
relative to and we find 


where fy and fg are determined by the formulas in Table 5 as a function of ¢;. Formulas (33) differ from the ana- 
logous formulas (17) for a system without lag in that they have the multipliers fy; and f,; (corrections to the lag), For 
convenience in determining the corrections, the formulas in Table 5 can be used to plot the curves fy; = f(¢) and {= 
= f(¢) for various r/T, and T,/T». 


For the case of a second order object with lag (4) and a controller of the type DDI (8) we have Y,(s) = Y',(s)e"™* 


1 
Y.(9) K x 
Ke + KP + s+ + ie 


Ks 
(Tys + 1) + 


(s) = (35) 


Originating from the requirements that the roots of the characteristic equation for the ACS (34) must be equal 
to the roots of the desired characteristic equation, we determine the roots of the characteristic equation (35), sub- 
stitute them into (34), and obtain the following solutions by solving the resulting equations for tr, and {¢; 


K 
_ T?— “Ko T ofa, 20:71 hea. (36) 


An additional condition is imposed on the transfer function; this condition is based on the presence of a third 
pole 


(37) 


t 


la = 


The formulas are written in the same form as that used for an object with no lag. The corrections fy and fy 
for the lag are determined from the expressions in Table 6, 


For convenience in determining the corrections from the formulas in Table 6 we plot the curves f,; = f(¢ 4), 
fez = f(¢ 1),and fgg = f(¢ for various values of r/To, T,/T and T,/T3. 


| 


34) 


37) 


Setting the right-hand sides of Eqs. (33) and (36) to zero and converting (33) to the form 


K K (38) 
7; = 7,2)» = 
and (36) and (37) to the form 


20 = 20) tn. 
20 log KT; = Dlog 


we plot the nomograms in Figs. 9 and 10 in the coordinates 20 log K/KgT, vs. ¢4 in order to determine the type of 
controller in an ACS with lag for fixed values of T,/T», €9 (for a second order object),and 1/T,g= 0.3, They are 
analagous to the nomograms in Figs. 5 and 6 for determining the type of controller in an ACS with no lag. 


Formulas (33) and (36) make it possible to choose the control parameters in such a way that the roots of the 
closed-loop ACS which determine the nature of the transient response correspond to the desired ones. However, in 
this case for the performance indices of the ACS with lag to be close to the desired indices of the ACS without lag 
(the transient response time differs by the amount 1), it is necessary that the remaining poles of the transfer func- 
tion for the ACS with lag lie considerably further away from the imaginary axes than the desired poles. Therefore 
there are limitations placed on the choice of the desired transfer function for the system with lag; these are deter- 
mined by the relationship between the object parameters and the desired transfer function for the ACS, The critical 
relationship between these parameters is determined by the minimal distance between the desired and the remaining 
poles for which it is still possible to estimate the transient response using only the principal poles, We adopt the phase 


stability margin for the open-loop system as a criterion which makes it possible to estimate the location of the poles 
of the transfer function of the system with pure lag. 


Under these conditions the formulas derived above are valid for those parameter values which cause the phase 
stability margin of the system with lag to be no less than the phase stability margin of the system with no lag. 


In order to determine the critical values of the parameters we computed the phase stability margins for the 
system with and without lag for various parameter values; we then compared the transfer functions of these systems 
for equal stability margins and found the relationships corresponding to the critical values of the parameters; i. e., 
K/KoT, = f(¢ 4;T;/To; €9) . The curves for the critical values were plotted in the coordinates 20 log K/KyT, and ¢, 
for fixed values of r/T, and T,/T 9; they were plotted on a nomogram for determining the type of controller in the 
system with lag (the dash-dot curves in Figs. 9 and 10), It should be noted that in Fig. 10 the curves for the additional 
conditions and the critical values coincide for the values T,/Ty = 0.5 and 1,0, Making use of the nomograms in Figs. 
9 and 10, it is possible to choose the desired transfer function for the system while taking the type of controller into 
account. Here we take into account the limitations governing the choice of controller parameters (these limitations 


were discussed above). After choosing the desired transfer function the controller parameters can be found from for- 
mulas (23), (33) and (36), 


We shall analyze an example which clarifies the method presented above. 
Example 4 (Fig. 11). The transfer function for the object is 


KWo (s) = $2.0.15-100 1° 


The following requirements are imposed on the transient response: Xmax=90.5; T=30 sec. We must deter- 
mine the type of controller and the numerical values of its parameters. We superimpose the nomogram in Fig. 10 
onto the nomogram in Fig. 4 and shift the axis of abscissas of the former by the amount 20 log Ky = 4 db. 


The desired transfer function is determined by the coordinates of the point A which is chosen on the intersec~- 
tion of the curve corresponding to the additional condition for T,/T» = 0.5 and the curve x,,,, = 0.5. 
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The desired transfer function is written as 


7.5s 
(5s + 1) (208° + 7.48 +1)’ 


Yq (8) = 


since T; = 0.5Ty = 5 sec; 20 log(K/T,) = 3.5 db; K = 7.5; 1 = 0.74; Tg = Ty. The parameters of the transient response 
afe Xmax = 0.5; T = 5 - 5.5 = 27.5 sec. The corrections fiz, fgg are found from the formulas in Table 6 (or from the 
graph): fy, = -1.02; fg = -1.5. The parameters of the controller are determined from formulas (23) and (36): 


K, =1/K = 0,132; 
The transfer function for the controller is 


KW, 1)/s. 


SUMMARY 


The proposed method permits the choice of the simplest type of controller and the determination of its para- 
meters either analytically or from the logarithmic frequency responses on the basis of specified performance indices. 
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DETERMINING THE OPTIMAL TRANSFER FUNCTION 
FOR "WHITE NOISE" IN A SERVOSYSTEM FROM 
THE FIGURE OF MERIT AND THE TRANSIENT 
RESPONSE DURATION FOR THE SYSTEM 


A. Kul'chii 


(Leningrad) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 524-529, April,1961 

Original article submitted July 24, 1960 


The paper establishes a relationship between the coefficients of the systematic error and the time 
required for the transient response when the minimum mean square error is ensured. Formulas are 
derived for determining the optimal transfer function for "white noise” in cases where the first 
two or three error coefficients are taken into account. 


A number of papers have been devoted to developing engineering methods for obtaining optimal block diagrams 
and dynamic~-accuracy estimates for automatic control systems. In [1], for example, a study is made of the problem 
of finding the optimal values of certain parameters for a servosystem with a specified structure on the basis of assur- 
ing the minimum mean square error (MSE) under conditions where it is equal to the square of the systematic error. In 
[2] the method of determining the optimal pulse transient response is generalized to include the case where the dynam- 
ic accuracy with which the regular useful signal is reproduced is determined by a priori specified coefficients of a cer- 
tain polynomial time function which characterizes this accuracy (the error coefficients C;). The sampling time T 
in [2] is chosen to be the transient response time for a unit step input. The amplitude-frequency and phase~-frequence 
responses are approximated semianalytically by means of the corresponding desired responses. 


In [3] K. I. Kurakin solves the problem of obtaining the responses for systems that are physically reliable and 
deviate the least amount the optimal system; this was achieved using an analytical method for approximating 
the transfer function @(p, e  *) by means of a fractionally rational Padeh function. The deviations from the optimal 
conditions are in this case also reduced to a situation where the error coefficients do not equal zero; i.e., the steady- 
state error with which the regular useful signal is reproduced is chosen equal to some allowable value. 


The satisfaction of the optimal conditions therefore assures a minimum MSE for a specified value of the 
systematic error. Thus, in order to determine the optimal transfer function and its corresponding MSE it is necessary 
to know the coefficients C; of the systematic error and the transient response duration T. 


In a number of cases the useful signal is a slowly varying time function; then the analysis can be limited to 
just the first two or three error coefficients, The coefficient Cy, which is a quantity that is the inverse of the figure 
of merit for the system with respect to speed, is usually specified; the error coefficients C, and C, in the cases in- 
dicated above have a comparatively small effect on the magnitude of the systematic error, and therefore it is ex- 
pedient to determine them from the conditions which govern the assurance of the minimum MSE. 


Based on the results of [3],we study the problem of determining the optimal transfer function for “white noise" 
in this paper; we achieve this by finding those values of the coefficients C, and C, for which the MSE is minimal. 
We consider C,; and T to be the specified quantities. 


In order words, the purpose of the paper is to determine the relationship between the accuracy of the system 
for random and regular inputs, and the time required for the transient response to be completed. 


By Original Equations 


We shall cite the expressions for the optimal transfer function and MSE of the servosystem in the case when: 
a) a stationary random noise of the “white noise" type with the spectral density Gpy (w) = €* = const acts at the in- 


put of the system, and b) the first two error coefficients satisfying the dynamic accuracy requirements are given to- 
gether with the transient response duration T [3]. 


The expression for the MSE in [3] becomes 
8C, 2C, 400,C, 
Here T is the duration of the transient response; C, and C, are the coefficients of the systematic error 6, 
dB (t Cc 
where 8(t) is the useful signal at the system input. 
_ From here on we shall operate with the relative value of the MSE which can be written in the form of the 
equation 
The equations for the optimal transfer function of the servosystem in the closed-loop pp (p) and open-loop 
Woprt,{P) states are given in [3] by means of the following expressions: 
grams 
lem T? 5 C T 2 
+ + pert 
T? 5 C. T 2 
nd For the system to be a stable and minimal-phase system in both the closed-loop and open-loop states the error . 
Pod coefficients Cy and Cy must satisfy the following inequalities: 
C,>0, (5) 
(6) 
sary 
(7) 
(8) 
ure 
n= For the case where it is necessary to take into account the first three error coefficients C,, C,, C,,the expres- 
™ sion for the optimal transfer function will be given by the equation [3] 
ise" 9 7C, 14C 3T? 14 C T 2C, 
m 
Under these conditions the relationships between the error coefficients must be determined by the set of in- 
equalities (6) and the inequalities 
7 97 3T 
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We shall now give the expression for the relative value of the mean square error for the optimal system when 
"white noise" appears at its input; this expression was derived on the basis of (9) using the formulas cited in [4]: 


2. The Regions of Allowable Values for C, 


From the inequalities (5)-(8) it is possible to plot the regions of allowable values for the error coefficient of 
the second derivative of the useful signal C, as a function of C, for specified T; these regions permit us to perform 
a direct estimate of the applicability of one or another value of C, from the point of view of the requirements im- 
posed on the optimal transfer function of the system. The indicated region for certain values of T are shown in Fig. 
1. Each of the regions is bounded from the left by the axis of ordinates (C, = 0), and above and below by straight 
lines whose equations are C, = 0.2T?—TC, and CG, = —TC,, respectively. On the right the region is bounded by the 
straight line Cy = 0.5 T. 


Knowing C, and T, it is thus possible immediately to determine the allowable limits for the values of C). 
However, this does not provide an answer to the question of whether or not we have chosen the optimal values of 
C, from the point of view of the magnitude of the MSE, 


The solution of this problem shall be sought, as indicated above, on the basis of the proposition that the error 
coefficients of the second and third derivatives have only a small effect on the magnitude of the systematic error; 
therefore they must be determined in such a way as to assure the minimum MSE, 


3. Determining the Optimal Transfer Function with Respect to the Figure of Merit 
and the Transient Response Duration 


Assuming that the magnitude of the figure of merit with respect to speed (1/C,) and the duration T of the 
transient response are specified, we find the expressions for the error coefficients which cause the relative value of the 
MSE tobe a minimum, For this purpose we use known rules to determine the minimum for the function (2) by assum- 
ing that the value of C, in it is variable; then we determine the minimum of the function (12), assuming that the co- 
efficients C, and C, are the variables. Having performed these operations, we obtain the following results. The func- 
tion (2) has a minimum for 


~ 
ai 
0 
| 
02 
48) 
~ — 
AS Sec; 
| 
Fig. 1. 
452 


It is evident that (13) satisfies both inequality (7) and inequality (8), Substituting (13) into (2), (3),and (4), we 
obtain the formulas for the relative value of the MSE and for the transfer function of the optimal system in the closed- 


loop and open-loop states in cases where it is possible to treat only the first two error coefficients. The enumerated 
formulas are written as: 


3c? 


Fig. 2. 


An analysis of the function (12) shows that it has a minimum for the values of the error coefficients C, and Cs 
which are defined by the equations 
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If we substitute (17) and (18) into (12) and (9), then we obtain formulas for determining the MSE and the transfer 
function for the system in the open-loop state in the case where it is necessary to take into account the first three 
coefficients of the systematic error C,, C,,and C, in finding the optimal parameters: 


196 (19) 
8uT 


T? T? 


The transfer function for the optimal system in the open-loop state is of the following form for the case under 
study: 
T? ling 1 


(p) = (21) 
Wope 
 + + + 720,? +! 


The transfer function (21) can be reduced to a form which is convenient for taking the logarithm; this makes 
it possible to plot the logarithmic frequency responses directly from a formula. Such a formula is written as follows: 


Vope (P) = 


(22) 


From expressions (13) and (14) we have plotted graphs in Fig. 2 which illustrate the dependence of the error 
coefficient C, (the solid line) and the relative MSE nj, (the broken lines) on the error coefficient of the first deriva~ 
tive of the useful signal for certain values of the transient response duration. 
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20) 


(21) 


4. Examples 


We shall study two examples of determining the optimal system from formula (16). Assume that it is necessary 
to obtain the characteristics of an optimal system for which it is sufficient to treat only the first two derivatives of the 
useful signal. The following quantities are specified: a) T =4sec, 
Cy = 0,025 sec ; b) T = 1.2 sec , Cy = 0.005 sec. Substituting the speci- 
fied values of transient response duration T and the error coefficient 


Cy into (16), we obtain the following expressions for the optimal transfer 
function: 


40 (1.82p + 1) (0.146p + 1) 


200 (0.045 + 1) (0.006p + 1) 
Wy = p (OAip + 1)(0.014p + 1) 


Using the transfer functions (23), we plot the logarithmic am- 
Fig. 4, plitude-frequency and phase~frequency responses (Fig. 3); curves 1 and 
3 apply to the system a, whereas curves 2 and 3 apply to the system b; 
A (w) is the amplitude response, and ¢(w) is the phase response. The 
shape of these characteristics in the low-frequency and intermediate frequency regions satisfies all the requirements 


that govern automatic control systems of high quality; in the high frequency region the characteristics must be cor- 
rected. 


Figure 4 shows the graphs for the transient response of the systems under study; these responses were obtained 
using an electronic integrator. Curve 1 corresponds to case a. and curve 2 corresponds to case b. 


SUMMARY 


1. Expressions are derived which make it possible to determine the optimal transfer functions when “white 


noise"appears at the input of servosystems that have specified values for the figure of merit and the transient response 
duration which assure a minimum MSE, 


2. The resulting formulas are convenient to use in engineering practice when using control systems. 
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ON THE PROPERTIES OF SINGLE-LOOP CONTROL SYSTEMS 


L. G. Sobolov 


(Leningrad) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 


pp. 530-535, April,1961 
Original article submitted November 18, 1960 


The stability conditions for certain types of single-loop control systems are given in the form of 
an interrelationship between the critical gain and the parameters of individual sections; this 
makes it possible to establish certain properties of these systems. 


Single-loop systems shall be defined as systems that consist of directional sections that are connected in series 
in a ring. We shall study systems that consist of aperiodic sections of a control circuit that controls an aperiodic ob- 
ject (case "A"), or an integrating object (case "B"). In both cases we shall also study so-called extreme systesm in 
which all the time constants of the control circuit are equal. 


A. The characteristic equation for the system is 


II (Ti + 1) + K =0, (1) 
i=o 


n 
where p= d/dt, K=]]X, i is the over-all system gain, T; for i > 0 are the time constants for the sections of the con- 
trol circuit, T; for i = 0 is the time constant of the object. 
For K = 0 the Mikhailov root locus for the open-loop system is 


D (jo) = 


n n 
where R=][ R,, = >) 0; +1, 0, arctan o7;, 
i=o i=0 


The maximum gain corresponding to the stability boundary is equal in magnitude to the radius-vector K,, = R 
of the root locus at its first intersection of the axis of abscissas (i.e., at 9'= #), We also have 


tan = 7, (é = 0,4,..., a) (3) 


for the condition 9'< 1/2 (since we are seeking the first intersection of the axis of abscissas by the root locus which 
occurs at the finite frequency w,). In solving (3) we also use the well-known trigonometric expression which can be 
reduced to the following form for n terms: 


4,j,....m=1 
(i<j<...m) 


uy => oT. 


| (2) 
Ca... 
arctan u,©arctan (4) 
i=1 C Cuf—... 
i=1 i=1 t=1 
where n n 
where uj" = > Um, 
| | 
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Cu=> Up uf = |] | 
i=1 


i=1 


n 
We shall also assume that (€ «= 1. On the basis of (4) we take into account the relationship 


tan >) 0; = —tan@, (5) 


i=1 


(which follows from the expression = and obtain 


i=1 


a" WA! 


/ 


i=0 i=1 f=} 


Fig. 1, 


Determining the minimal root w, from (6), we find the critical gain of the open-loop system. 


| 


Thus, for example, for a third order system ( n = 2) we have 


By performing an investigation of the extremum for (7) for T'; = var, but T'; = const we find 


T,=VT,, (8) 


i.e., the worst case from the point of view of stability isa square-law relationship between the ratios of the time con- 
stants for the control circuit sections and the object time constant. 


> 


Fig. 2. Graphical solution of the characteristic equation for an 
extreme single-loop system. 


An investigation of the extremum for T'; = var and T'; = var leads to the known [1, 2] relationships T', = T'; = 1 
corresponding to the minimum-minimorum K,,. The relationship K,, = K,,(T', T'z) is shown in Fig. 1. 


Note the interesting property of third order systems: 


(1+T,)? 
(Koy = (Ke =2 (9) 


For an investigation of the properties of extreme systems it is convenient to use the critical nonuniformity a,, 
of the control circuit. It is not difficult to show that 


| | = cos” cos (n§), (10) 


where 9 is the solution of the equation 
tan 9 = — Ttan(nQ) (11) 
for the condition 9 < #/2. Here T is the ratio between the time constant of the control circuit section and the object 


time constant. 


In the coordinates aT the curve for the critical nonuniformity of the open-loop system represents a family of 
single~parameter curves (the parameter is n) emanating from the origin. The maximum for a curve of any order can 
be found from the condition 


8 


= 
T 


where 
| 
6 
| 
2 
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=0, (12) 


which leads to the expression 


tan § = —tan(n§). (13) 


A comparison of (11) with (13) shows that the maximum for a curve of any order occurs for T = 1. 
Solving (10) and (11), we obtain the following result 


for a third order system (n = 2): 
Ker 
20 
5, 0055 and for a fourth order system (n = 3) we obtain 
00s— 


The nonuniformity of the closed control loop can be 


determined from the known relationship 


a 
N —S aN The maximum for the curve 6, of an nth order 
system also occurs for T = 1. 
— The expression for the conditions governing aperiodic 
\ foie. er stability of the extreme system can easily be obtained if 
6 Nei ‘= the characteristic equation for the system is written as 
46 where 
t 


3 It is obvious that the boundary for aperiodic stability 
nt corresponds to the case shown in Fig. 2. The ordinates corre- 
sponding to the points of tangency between y; and yg charac- 


terize the minimal nonuniformity of the open-loop circuit corresponding to the boundary of aperiodic stability. There- 
fore 


(n + T 


In extreme systems aperiodic stability is possible only for systems of order no higher than the third. 
B. In analogcus fashion it is possible to derive the expression 


hed 
Koy = Ton |] V +1 (18) 


ire 


n 
T T] (19) 


ni=1 


|= 


where 9;; and w, are found from the system of equations 


n 
2 >) tanQ,—07, (i= 1,2,...,n) (20) 


for the condition 9; < 1/2. 
The system (20)is reduced to the equation 


for the specified conditions. 


Determining the minimum roots (21) and making use of (18), we obtain 
Keg = Fit (22) 


for third order systems. 
The function Kcr = To/T:) is shown in Fig. 3. 
For extreme systems it is more convenient to use (19) and (21); on this basis we find 
cos" 


Sor |= (23) 


Value of the 
0.5 |1.125] 1.76 


2.39 | 3.03 | 3.67 | 4.30) 4.94) 


It is evident from (23) that the boundary stability curves in the coordinates §T consist of rays emanating from 
the origin with slopes depending only on n; the magnitudes of these slopes are equal to the values of the derivatives 
of a,, for T = 0 and for the corresponding value of n. The values of the slopes are given in the table as a function 
of the order of the equation (i.e., as a function of n + 1). 


The expression for the aperiodic stability curve for the extreme system is given by 
Sap = AT, (24) 
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Aperiodic stability is possible for systems of order no higher than the third. 


Expression (23) can also be recommended for determining the minimum achievable nonuniformity for non- 
extreme single-loop systems of arbitrary order. 


In conclusion, note that the transient responses for the unit step perturbation in an extreme single-loop system 
of nth order are completely determined by the two dimensionless quantities T and K (or T and 8); this makes it 
possible in principle for formulate transient response tables for the system. 
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ON AUTOMATIZING THE INTRODUCTION OF CERTAIN 
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Original article submitted November 29, 1960 


The paper studies a method for the automatic transformation of functions specified in the form of 
graphs, oscillograms, or other visible images into digital (tabular) form. Here the possibility of 
direct utilization of the resulting data in digital computers is assured for the purpose of solving 
many problems (for example, for finding interpolation polynomials, for performing numerical in- 
gration and differentiation, for solving a number of differential equation, etc.). 


In solving one or another type of problem on automatic digital computers it is necessary to reduce the data 
introduced into the machine to tabular form. If this data is initially in the form of graphs, oscillograms, photographs, 
etc., then the procedure of converting it to tabular form by means of a human operator occupies a great deal of time; 
this may be impermissible because of the current tendency toward increasing the speed of computation (cf. [1, 2]). 


Here it is expedient to study the possibility of automatizing the conversion of data from graphical to tabular 
form. 


Principle of Performance of the Transformation 


The graph (oscillogram) of a function may be continuous or quantized; a table is of necessity always quantized. 
This makes it possible to use a type of “retina” system for “seeing” the image; this system is made up of a matrix of 
individual light detecting elements (receptors). For functions which are given in rectangular coordinates it is con 
venient to use a rectangular matrix. In such a matrix the distance between the centers of two neighboring receptors 
on one line (or in one column) can serve simultaneously as a unit length and the measure of the resolving power. 
Based on this approach, the intervals between neighboring values in the table can contain only an integer number of 
unit lengths, and the minimum length for such an interval is equal to one unit length. Because of the fact that the 
functions to be processed f(x, y) = 0 may in the general case be functions that are not single~valued relative to x, 
it is expedient to use a scanning method in which the reference points are taken not along the x axis of the graph, 
but along the line proper. The table formed in such a case as a result of “reading” the graph corresponds to the ex- 
pressions xj = f;(s;), yj = (sj); i.e., the table "writes" the curve in the form of a pair of natural equations. In prac- 
tice this means that quantities referring to adjacent reference points on the curves are written in adjacent cells of 
the table, whereas in taking reference points along the x axis this proposition could be violated. For the engineering 
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realization of such a scanning method it is necessary to solve the problem of incorporating into the conversion device 
the capacity to establish whether or not two specified activated receptors correspond to neighboring (along the curve) 
reference points s; and 


- 
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In order to be specific we shall choose the minimum interval between reference points s;,, —8;; i.e., this 
distance will be chosen equal to the distance between the centers of two adjacent receptors. It can be demonstrated 
that the reasoning which follows below for the case of the minimum interval between reference points remains valid 
even for larger intervals, provided that these intervals are chosen according to definite rules and the curves satisfy 
certain conditions (the Deeny conditions) that are satisfied in the majority of cases [3, 4]. In that case a satisfactory 
method of solving the problem is one which might be called the method of formulating chains of numbers. The 
essence of this method is the following. 


All possible pairs of numbers of adjacent retina receptors are introduced into the memory section of the unit. 
We shall define two receptors as neighboring ones if their centers can be connected by a line segment which passes 
only through areas occupied by these two receptors. It is evident that the number of different pairs of numbers for 
adjacent receptors is equal to the number of different segments of a straight line connecting the centers of adjacent 
receptors; such a method for counting pairs makes it possible to assume pairs of number (i, j) and (j, i) to be iden- 
tical. 


The signals arriving from the retina which contain information on the numbers of the activated receptors are 
compared with the pairs of numbers for the adjacent receptors which are stored in the memory. Only those pairs for 
which both numbers correspond to activated receptors in the specified image are passed through for further processing. 
The formulation of the number chains is accomplished in the following manner. 


Any of the activated pairs is examined. Assume that the left number in that pair is taken, and among the re- 
maining activated pairs we find a pair which also contains a number. The second number for this pair which we 
have found is written at the left next to the original pair. Then we take the right number of the original pair and 
repeat the same process. The second number in the pair that we have found is written next to the original pair on 
the right. After that the entire procedure described above is performed for the left and right extreme numbers of 
the resulting chain of four numbers, etc. The process can be stopped in two ways: 1) we come to the end of all the 
activated pairs; 2) we find no pairs containing the extreme numbers in the chain among the remaining activated | 
pairs. The case where the process is stopped for the first reason corresponds to a graph containing one curve with no 
discontinuities. In that case we may encounter two variants: a) if in comparing the extreme numbers of the chain 
these numbers prove to be identical, the curve is closed and when necessary a signal is given to that effect; b) in the 
contrary case the curve is not closed. If the process has been stopped for the second reason, the chain counter re- 
cords the termination of one chain; then any of the remaining activated pairs is taken and the next chain of num- 
bers is developed from it. If activated pairs still remain, then the process is repeated, etc. Note that after ter- 
mination of each chain it is also possible to check for closure of the curve. 


For the case where the graph of the curve forms an intersection, the number of the receptor corresponding to 
the intersection point will belong to more than two activated pairs. The programming for the unit in this case will 
depend on the particular features of the investigated material. If there is no need to pay any special attention to 
the act of intersection, we can limit the process to finding only one activated pair containing the required number 
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and continue the formulation of the chain, However, if it is necessary to fix both the fact of intersection and the 
intersection point, it is necessary to check all the activated pairs each time in order to choose from them those which 
contain the specified number. If we find there is more than one such pair (not including the pair which is already 
incorporated in the chain), then the intersection point is memorized and one of these pairs (any one) is entered into 
the chain which is being formulated, The other pairs are entered into other chains. Then it is possible to compare 
the intersection points which are fixed in the memory unit for all the chains and to establish the number of inter- 
section points and the number of intersecting chains at each point. 


Each number which is incorporated in the chain contains the coordinates x; and y; of the corresponding ith 
receptor in coded form. Thus the result of recording each chain actually represents the sought-for tabular recording 
of the function. 


We shall make a brief examination of methods for achieving an engineering realization of the principles pre- 
sented above. The memorization of the pairs of numbers can be achieved without a memory unit in the conventional 
sense. In cases where the required memorizations of the information remain permanently the same, the memory 
function may be performed by the network structure proper. In fact, any network contains in its structure the in- 
formation on the rules according to which the signal is transformed in the network. 


In our case it is possible to use a well-known network for choosing the chain; this network is often used for 
converting numbers from the binary system to some other system (cf. [5]). An example of a portion of such a net- 
work for four receptors is shown in the figure (a). Activation of a receptor in this network corresponds to a positive 
voltage at its output. For the curve shown in Fig. 6 a positive signal appears at the output 3, 4 (i.e., the pair 3. 4 
is activated and is transmitted for use in the operative memory). It is easy to see that all of the remaining processes 
described above can also be accomplished using universal and widely used computer elements: coincidence net- 
works, noncoincidence networks, shift registers, etc. It can also be demonstrated that when necessary the memory 
network (a portion of which is shown in Fig. b) can easily be subdivided into individual parts that “service” the 
corresponding sections of the retina; moreover, in a number of cases these retina sections can be “serviced” alter- 
nately by one partial network. 


APPENDIX 
1. The number of pairs of adjacent receptors can be computed by an obvious method. For a quadratic retina 
consisting of A? receptors the number F of adjacent receptor pairs is equal to 


B = 2(A—1)(2A—4). 
Thus for a retina consisting of 100 x 100 receptors, B 4-104 


2. The method for indicating the receptor numbers, Using the binary number system for denoting the numbers 
(coordinates) of the receptors, it is possible to make use of two \aethods for denoting the numbers. 


A, Assume the retina contains $ = A® receptors. Then in order to denote any receptor number in the binary 
system (i.e., to denote any number from 1 to S) it is sufficient to use binary numbers containing logsS places. If 
the curve contains L unit lengths, then its notation based on this method will require L log, S binary places. 


B. Assume that S-place numbers are used to denote lines on this same retina; then the notation for any curve 
will contains S binary places, irrespective of its length L. 


It is evident that the use of method "A“* for method "B" depends on the relationship between the average length 
of the investigated curve and the dimensions of the retina: 


ifS > Lay logeS (i.e., Lay < S/iog, S), then we should use method "A"; 
if S< L,y logeS (ie., Lay > S/log,S), then we should use method “B." 


For the case where the image contains not only a line but also areas which must be taken into account, the 
method "B" will as a rule be preferable. 
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The paper describes the principles of operation for a wire frequency resonator that is used for wa 
industrial telemetering. The basic design relationships are cited. The results of experimen- om 
tal investigation are cited along with the circuit diagram for the wire frequency resonator. 


One of the most reliable forms of information used in telemetering is the frequency of an ac current whose 
magnitude depends on the value of the measured parameter. During the last few years a large number of instru- 
ments have been designed according to the frequency principle for the purpose of telemetering resistive and reactive 
power, voltage, current. and other quantities. However, a shortcoming of such instruments is the low accuracy with 
which the measured quantity is converted into the frequency of the ac current. The RC and LC oscillators which are 
used as transducers to convert the parameter into a frequency have appreciable conversion errors which are caused 
by variations in the temperature and humidity of the ambient medium, atmospheric pressure, and supply voltage. The 
electrical properties of the elements in the tank circuit of such oscillators are unstable in time. Therefore frequency 
methods of telemetering are at present used in those cases where the communication channels employed can pass 
only quantized information without substantial distortion. Examples of such communication channels are high-volt- 
age power lines, radio-communication channels, and conventional overhead or cable lines whose length amounts to 
tens or hundreds of kilometers. 


The growing requirements imposed on the accuracy of industrial measurements, coupled with the develop- 
ment of digital control machines for which the most convenient form of introducing information is the quantized 
form (pulse-code, digital-pulse, frequency), have increased the interest in transducers which convert industrial para- 
meters into frequency. 

The design of such accurate and reliable transducers will make it possible to realize the advantages of frequency 


methods over telemetering methods based on identifying features related to intensity in using communication channels 
with a high noise level. 


The conversion of a frequency into readings of a secondary instrument using digital indicators and printers, as 
well as the conversion of frequency into codes which can be used to introduce information into the control machine, a 
do not at present lead to any particular difficulties. 


Below we shall describe certain results of investigations performed on a wire frequency resonator designed to 
convert a parameter into frequency; these investigations were performed by the authors during a process of design- 
ing a digital telemetering system. 


The wire frequency resonator (Fig. 1) consists of the metallic wire 1 with one end is fixed to a support 2, The 
other end of the wire is connected to the kinematic element of the metering instrument and is not coupled rigidly 
to the frame of the resonator. The force F which stretches the wire is determined by the value of the measured para- 
meter. The proper frequency fy of the wire is determined by the relationship [1] 


where 1? is the length of the wire, p is the linear density of the wire. 


3 
Fig. 1. Fig. 2. Fig. 3. Fig. 4. 


From [1] it is evident that for a resonator in which? and p are constant quantities,the proper vibration fre- 
quency of the wire is determined by the force F. In order to excite vibrations in the wire the latter is placed in a 
transverse magnetic field produced by the permanent magnet 4; the wire is coupled to an electronic amplifier in 
such a way that the system forms an oscillator with self-excitation in which the resonator is the wire resonator. 


k =kiloohm 


The vibrations of the wire are described by a second-order partial derivative equation for small amplitudes; this 
equation can be reduced to an ordinary differential equation on the assumption that: 1) the vibrations of the wire are 
sinusoidal, and 2) the distribution of the magnetic flux density along the wire is sinusoidal. 


As a result of these simplifications it is possible to develop an equivalent circuit for the vibrating wire; this 
equivalent LC tank circuit is shown in Fig. 2., and its parameters are determined by the following equations [3]: 
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(2) 


2p (3) 
C= 


where B is the maximum value of the flux density, k is the coefficient of air friction. 


It is evident from the equivalent circuit that when the wire is connected directly into the feedback loop of the 
amplifier a positive feedback independent of frequency arises due to the resistance r of the wire. In order to elimi- 
nate this feedback the wire is connected into the arm of a resistance bridge made up of the resistors Ry = Rp and R, = 
=r (Fig. 3). The bridge is balanced at frequencies far from the resonant frequency. This circuit is more or less 
analogous to the corresponding circuit for the connection of a quartz oscillator [4]. Just as in the quartz oscillator 
circuit the selection of the magnitudes of the bridge resistors R,, Rg,and the normalized input resistance of the am- 
plifier is made in such a way as to avoid any substantial lowering of the Q for the wire. In order to reduce the varia- 
of the phase shift when the temperature changes, the resistance R, is made of the same material as the wire. The Q 
of the freely vibrating wire can be determined from expressions (3) and (4): 


F 
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The linear dependence of the Q for the wire on frequency is verified experimentally (Fig. 4). This data is 
cited for a tungsten wired = 0.05 mm in diameter and / = 20 mm in length vibrating in air, However, relationship 
(5) is not completely exact since it does not take into account the viscosity losses in the material of the wire and 
the friction at its points of attachment (terminals). It should be noted that one of the basic factors affecting the 
constancy of the vibration frequency for the wire is the method by which it is fixed to its supports. The microscopic 
unevenesses which are present on the surfaces of the terminals may lead to an appreciable lowering of the Q due to 
friction of the sections of the wire that are not clamped down. It is evident from (5) that it is advantageous to use 
wire materials with a maximum value for the ratio F/p. The best material for such wire is tungsten. 


In principle the wire resonator has a nonlinear characteristic which is determined by relationship (1). Since 
the Q decreases when the frequency is lowered, a certain initial force must be applied to the wire in order to assure 
. stable vibrations. Correspondingly, the resonator will have an initial frequency for the zero value of the parameter. 
However, in order to obtain beats the frequency corresponding to the zero parameter value should be taken as the 
maximum; then the frequency is decreased as the parameter increases. Under these conditions the frequency must 


either vary over wide limits relative to the zero frequency, or the frequency deviation must be several percent (depend- 


ing on the requirements imposed on the resonator). In the first case the nonlinearity is eliminated by the nonlinearity 
of the primary pick-up or by using profiled cams. 


For a small frequency deviation the nonlinearity § can be determined using expansion into a series 


where Fy is the initial tension, AF is the change in tension which is proportional to the measured parameter, 


for AF << Fy. 


An appreciable decrease in the nonlinearity occurs when beats produced by the frequencies of two wires Af = 
= fh. 
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The proper frequencies f, and f, are determined by an identical initial tension Fy for both wires and by the dif- 
ferential change in tension SF. The nonlinearity is determined from (6) and is equal to 


1 (AF \2 
which yields an error of 0.03% of the tension variation AF for a frequency deviation of 2 5 % 


We made an all-around experimental investigation of a wire resonator with the following basic parameters: 
1 = 20 mm, d = 0.05 mm, B = 500 gauss. 


The working frequency variation is assumed to be from 2 to 6 kc; however, during the process of the actual 
experiment the frequency varied from 1.5 to 7.5 kc. 


The stability of the frequency for fixed tension produced by means of a weight was monitored over a period of 
10 days for t = 20° 5°. 


During that time the maximum frequency deviation was less than 0.2% of the maximum working frequency of 
6 kc. The temperature dependence of the frequency was measured in the range from 20 to 85° C and was 0.09% per 


10 °C. In principle the resonator can be used right up to the highest temperatures whose limits are determined by 
the materials with which the device is constructed. 


Such a high temperature stability can be achieved if the wire resonator is coupled to the primary metering 
instrument by means of a spring, a membrance, or some other elastic element whose rigidity is very much less than 
the rigidity of the stretched wire. 


The error due to the elongation of the wire when the temperature changes is determined by the coefficient of 
linear expansion for the wire material; it is appreciably less than the over-all temperature error of a wire resonator. 


The use of tungsten wire assures stability of the length and diameter of the wire with time. 


A photograph of the resonator with its cover removed is shown in Fig. 5; the dimensions of the resonator are 
32 x 32 x 44mm. 


The two-stage amplifier is designed using transistors of the type P13 and P14 (Fig. 6). The matching trans- 
formers at the input and output of the amplifier use armored cores made of ferrite F-1000 20mm in diameter. For 
temperature variation in the range from +20 to + 80° the amplifier maintains its stability; this is achieved by stabiliz- 
ing the operating mode of the first stage [5]. The power drawn from the supply circuit is 15 Mw for u = 6.3 v and 
can be lowered appreciably by using a diode amplitude limiting circuit. A variation of the supply voltage within the 
limits +50 % does not affect the accuracy of the resonator. 


We measured the frequency by comparing it with the frequency of a standard quartz oscillator. The com- 
parison error did not exceed 4 0.1% 
SUMMARY 


The high Q of the described wire resonator, coupled with the large range of possible frequency variation, makes 
it possible to use it for the accurate telemetering of various industrial parameters. The parameters can be simply con- 
verted into a digital code. 


Depending on the required speed of response and passband widths of the communication channels, the resonator 
can be designed for any frequency band from 100 cps to 10 kc. The resonator is small. 


The high reliability of the resonator is assured by the absence of moveable elements which wear, by its simple 
construction, and by the use of transistor amplifiers, 
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The paper determines the frequency responses of carbon pile voltage regulators for dc and ac 
generators. Based on an analysis of the transient responses performed in a linear approximation, 
expressions are derived for the transfer functions of the regulator and for its corresponding fre- 
quency responses. A method is developed for the experimental determination of the regulator 


frequency responses in a closed-loop generator-regulator system. Results are cited for an ex- 
perimental investigation of one of the standard manufactured voltage regulators. 


1. Introduction 


Notwithstanding a number of shortcomings of carbon pile voltage regulators which are due to the comparatively 
rapid wear and low heat resistance of the carbon pile, such regulators have been widely used for the voltage regula- 


tion of dc and ac generators, and for converters and rectifiers [1, 2]. Nowadays,carbon pile regulators are being re~ 
placed by more reliable voltage regulators based on self-saturating reactors and semiconductor devices. However, in 
a number cases where the cooling conditions are not the limiting factor carbon pile voltage regulators will still re- 
main in use. 


Based on the above, a study of the frequency responses of carbon pile voltage regulators and of the engineering 
method for determining those responses is of definite interest. 


The transient responses in carbon pile voltage regulators have been studied in a linear approximation by a num- 
ber of authors [1,3, 4 etc. ]. However, in a majority of cases the study is a simplified one: it does not take into 
account the electromagnetic transient responses [3], or it does not consider the effect of the load equalizer [1, 4]. 
The problem of an experimental determination of the frequency responses for carbon pile regulators has been treated, 
as far as we know, only in one paper [4]. However, the method proposed in this paper for determining these res- 
ponses is so complex that it cannot be applied in engineering practice. A special feature of this method is the fact 
that the responses are taken for an open-loop regulation circuit, and therefore a complex electronic device of the 
“hard” feedback regulator type is necessary for stabilizing the initial mode of operation of the regulator. 


As will be demonstrated below, for obtaining the responses of a closed-loop regulation system the stabilization 
of the initial mode of operation is achieved by the regulator proper, and no auxiliary units are required. 


A theoretical investigation of the transient responses of a carbon pile regulator shall be performed for a typical 
voltage and load regulation circuit for parallel operation of dc generators (Fig. 1). The connection of the carbon pile 
voltage regulator for ac generators differs from the circuit shown only in that the regulated ac voltage is first recti- 
fied and then applied to the working winding W,, of the electromagnet. 
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The signal from the reactive load equalizer is usually introduced into the circuit of the working winding at the 
input of the rectifier [1]; this simplifies the analysis of the system. In a number of cases generator excitation cirrent 
equalizers are used for parallel operation of ac generators instead of reactive load current equalizers. In that case 
the equalizer circuit does not differ in principle from the circuit shown in Fig. 1. It is evident from Fig. 1 that the 
electromagnet of the regulator (usually of the armored type) is equipped with two windings—the working winding W,, 
and the equalizing winding Weq . The equalizing circuits of the regulator are connected in parallel for parallel 
operation of several generators; this assures load regulation of the generators in accordance with the method of imagi- 
nary statism. 
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Fig. 1. Circuit diagram for the carbon pile voltage regulator. 1) 
generator; 2) carbon pile voltage regulator (a is the electromagnet, 
b is the carbon column); 3) stabilizing transformer; 4) calibrated 
resistors. FW is the excitation winding of the generator; R, is a 
remote resistor for setting the regulator; Reg, Regr,and Ret are the 
resistances of th: windings on the stabilizing transformer. 


Besides these two windings, simple carbon pile voltage regulators are usually equipped with several more 
auxiliary windings which include a temperature compensation winding and a corrective winding (Fig. 2). 


The corrective winding is not shown on the original circuit diagram (Fig. 1), since,due to the presence of the 
rectifier D, (Fig. 2) in the circuit of this winding,it operates only for large excitation currents (for small values of 


the carbon pile resistance R,,). However, these modes of operation for the voltage regulating system are not of great 
interest from the point of view of stability, since in studying the transient responses the effect of the corrective wind- 


ing can be neglected. 
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Regulator circuits usually include a multiplier resistor R, which introduces “hard” negative feedback. In the 
RUG~-82 circuit (Fig. 2) this resistor is connected through the rectifier D; in such a way that for large excitation 
currents (i.e., for small generator speeds) no current flows through the resistor; this eliminates any additional lowering 
of the generator voltage in such a mode due to the presence of “hard” feedback. 


equalizer 
bus i 


Fig. 2. Electrical circuit diagram of a type RUG-82 regulator com- 
bined with a TS-8 stabilizing tansformer: a) complete circuit; b) 
section of the circuit which applies to the working winding. G is the 
generator; Weor is the corrective winding of the regulator; Wy is the 
temperature compensation winding; TS-8 is the stabilizing trans- 
former; BS-18 000 is a calibrated resistor; VS-20 is the remote 
for setting the regulator (setting the regulated voltage); Dy and D, are 
selenium rectifiers. 


The voltage regulator circuit includes stabilizing transformer which in the general case has four windings (Fig. 
1): Wy, Wg, Ws, and Wg they are connected as shown in Fig. 1. 


2. The Equation for the Transient Responses of the Regulator 


For a theoretical investigation of the regulator transient responses we shall make the following basic original 
assumptions. 


1, The deviation of the variable quantities from their steady-state values are assumed to be small; this corre~ 
sponds to the problem of developing a method for determining the frequency responses of the controllers for small 
sinusoidal perturbations. 


2. The leakage fluxes of the windings on the regulator electromagnet and stabilizing transformer are neglected. 


The resistance R,, of the carbon pile depends, as we know, on its compressive force; the latter is determined 
by the magnitude of the air gap § between the armature and core of the electromagnet. Thus, R,, = f (8). 


For small increments of the variables it is possible to write 


AR, =—Kphd, 
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The equation of motion for the electromagnet armature when small deviations of the variables from their 
steady~state values occur is written in the following form in operator notation for zero initial conditions: 


(mp* + Kap + K,) 40 + KgA®Dy=0, (2) 


oF, 
where K,=Ky, +X, ~ wh is the tangent of the slope angle for the mechanical characteristic of the 
OF, 
regulator = = f(6) at the point Ky = | , F, is the traction force of the electro- 
magnet; ,, is the magnetic flux in the electromagnet core; ®, = f(6, aWw)= iwWy + leqWeqtpWry: K is the co~ 
efficient of viscous friction; m is the over-all reduced mass of the moving sections of the electromagnet, and Aé is 
the air-gap increment. 


From equation (2) we obtain 


M (Pp) 


M (p) = Tp + p+4, 


The transient responses in the windings of the electromagnet can be described by the equations 


ok 


+ y— = (ut ty — igh), 
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(5) 


where u is the generator voltage (the regulated voltage); u,, is the voltage across the terminals of the working winding 
circuit of the regulator (cf. Fig. 2); u,, is the voltage across the excitation winding of the VA yi Nyaa 
age in the equalizing circuit which is equal to the potential difference between the equalizing bus and upper 
terminal of the calibrated resistor (cf. Fig. 1); Rwy, Req are the resistances in the circuits of the working and equaliz- 
ing windings (taking into account the resistances of the stabilizing twansformer windings); @g_ = f(awg,) is the mag- 
netic flux in the stabilizing transformer, aWg; = iW; + iestW2 ~ieqWs ~iwwe 


Varying the system of equations (5) and solving it relative to AS wy, we find the sought for equation for the 
electromagnetic transient responses in the regulator: 
w. 
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where 
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When the equalizing circuit is cut out, we have i., = 0. 
Setting Reg = e in (6) we,find 
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Fi (p) =1 + Tate) + Ey 
Ty — = Se 2+ Tot 


From Eqs. (6)and (9) it is evident that the resistance R, introduces "hard" feedback with respect to the voltage 
Un into the circuit along with "soft" stabilizing feedback; this worsens the accuracy with which the generator volt- 


age is regulated in steady-state modes. 
Therefore the feedback coefficients Kg, = Ry/R, in standard manufactured regulators is made very small (0.013- . 
0.03). Thus in investigating the transient responses the feedback through the resistor Rg can be neglected. 


Solving the system of equations (1), (3),and (6) relative to AR... , we obtain the general equation for the 


regulator (for Ry/R, = 0): 
AR, + G,,Au,, + G,Ai ® (10) 
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The functions Gy, (P), Geq(P), Gr(p).and G,(p) are the transfer functions for the regulator with respect to the 
corresponding input signals. 

In analogous fashion, for i.g = 0 (Req = =) we obtain 

AR, =G,,Au+G,,,Au,, + G,,hi,, (12) 
3 (P 1 wm wy 
(0) =F, (p) M (Pp) — Ke (A + 4 + PT = sts 

The functions Gw; (p), Gm (p)and Gay (p) are the transfer functions for the regulator with respect to the corre- 
sponding input signals for an open equalizer circuit. 

Replacing the operator p by jw in the expressions for the transfer functions, we obtain an analytical expression 
for the corresponding frequency responses of the regulator. In particular, for ieg = 0 and no stabilizing tansformer 
we have 

Guy Go) = — Ke) + + Te] + — (14) 

In order to use the analytical expression for practical computations it is necessary to know the regulator para- 
meters. However, the experimental determination of these parameters is undoubtedly a much more complex problem 
than the experimental determination of the frequency responses proper. 

It should be noted that experimentally it is sufficient to determine the frequency responses G,{jw) and Gwdjw) 
since the remaining responses can be computed if the indicated frequency responses of the regulator and the stabilizing 
transformer are known, Thus, for example, from (13) we obtain 

Cet, = pMx PT sta 
Rost (1 + PT st.) ~ PT (16) 

The frequency response Gg;+(jw) is determined as the ratio between the complex amplitude of the emf across 
the terminals of the winding W, of the stabilizing transformer and the complex amplitude of the sinusoidal com~- 
ponent of voltage across the terminals of the primary winding Wz when the remaining windings are energized: 

(io) = for i,=0, iw=9, eq =9, u, = + AU, sin wf. 

Since the stabilizing transformer is not saturated, it follows that in determining the frequency response we 
can assume Uno = 0. Having determined the time constant Tg from the amplitude- and phase~frequency responses 
of the transformer, and knowing the frequency response of G,,(jw) and Gw;(jw) we can easily obtain all the other 
frequency responses of the regulator. 

3. Experimental Determination of Frequency Responses for the Regulator 

The simplest method for experimentally determining the frequency responses of the carbon pile regulator con~- 
sists of measuring its steady-state response (the amplitudes and phases of the fluctuations in the carbon pile resistance) 
to constant-amplitude sinusoidal input signals at a number of fixed frequencies. 
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Fig. 3. Circuit diagram for obtaining the frequency responses of a carbon pile voltage regulator when os - 
cillations are excited from the regulator side: 1) generator; 2) regulator; 3) stabilizing transformer; 4) 
low-frequency signal generator; Rg, f4y. Ishg are noninductive series resistors; R, is the remote resistor 
for setting the regulator. 


Fig. 4. Circuit diagram for obtaining the frequency responses of a carbon pile voltage regulator when os- 
cillations are excited from the generator side: 1) generator; 2) regulator; 3) stabilizing transformer; 4) 
low-frequency signal generator; 5) rotary amplifier. 


Here we must keep in mind the fact that since the carbon pile regulator is a nonlinear device its parameters 
(ws Kp, Kg. Kyy etc.) and frequency responses depend to a substantial extent on the initial unperturbed state (69, 
Reta). 


Therefore in order to determine the frequency responses it is necessary to create a mode in which small forced 
oscillations of the electromagnet armature occur relative to the intial state (6 9). 


The high values of gain (for small AK, = 1—K,), and the more so for astatic adjustment of the regulator (Ke * 
= 1), such a mode can be assured when feedback is present between the carbon pile resistance and the voltage across 
the working winding of the regulator; i.e., it can only be achieved in a closed-loop voltage regulator circuit. Taking 
the frequency responses of the regulator as an individual element outside of any regulation system proves to be im- 
possible, since the smallest variations in the intial voltage across the working iiRion (Ug) leads to full compression 
or expansion of the carbon pile. 


In the paper by Mills [4] a special electronic regulation circuit was used to obtain the frequency responses of 
the regulator; this special circuit assured stabilization of the initial value of the carbon pile resistance (Roy) and 4 
mode in which small oscillations in the carbon pile regulator occurred. 


However, the use of such a method is unjustifiably complex, since the mode of small forced oscillations can 
easily be produced in the closed-loop generator-regulator system without using cumbersom auxiliary electronic 
devices, Under these conditions the stabilization of the initial regulator parameters is assured by the regulator itself. 
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The excitation of small sinusoidal oscillations in the voltage regulation system can be achieved in two ways: 
a) from the regulator side (Fig. 3), and b) from the generator side (Fig. 4). 


In the first case a small series resistor R, is connected into the over-all regulator winding circuit; the low- 
frequency signal generator 4 is connected to the series resistor through an isolating capacitor C. The signal supplied 


by the generator causes the regulation system to enter a mode of forced small oscillations about the specified steady - 
state mode. 


In order to determine the frequency responses of the regulator it is necessary to maintain a constant amplitude 
for the voltage Atreg across the regulator at each frequency, and to measure the amplitudes and phases (relative 


to Atreg) of the sinusoidal components of the voltage drops across the carbon pile (Afi, ) and across the resistor 
Ig, connected in series with the carbon pile (Aly,). 


For small sinusoidal oscillations in the system we have (the switch S2 is closed, and S1 is open): 


AR, = Gy(jw) Aureg 
Ad, = AR, Ing + co, 


Depending on the position of the switches S1 and S2, we can thus determine the various characteristics of the 
regulator. 


1) If the switch S1 is open and the switch S2 is closed, we obtain the characteristic Gw(jw). 
2) If Sl and S2 are open, we obtain Gw;(jw). 
3) If S1 is closed and S2 is open, we obtain G'w,(jw). 


It must be noted that the indicated circuit can also be used for obtaining the frequency responses of the entire 
voltage regulation system as a whole. 


In that case the equalizing circuit of the regulator and the windings W, and W, of the stabilizing transformer 


are cut out as shown in Fig. 1. The regulator equation for this circuit can be written in the following form in terms 
of small increments: 


= Gw(Auger—Aujn) + Gc Aue +G,,Au,, + (19) 


In obtaining the frequency responses for the entire system the amplitude of the voltage across the series resistor 
(4u;,) is maintained constant; the series resistance R, under these conditions must be made small (1 to 2% of Rwo) 


in order that its connection into the regulator circuit should not cause a noticeable variation of the initial steady- 
state operating mode of the generator (Up). 


Small oscillations can also be excited in the system when the frequency responses of the regulator and of the 
entire system are tkane even when the series resistance is not connected into the regulator circuit; this is done by 


connecting a low-frequency signal generator 4 to the remote resistor Rr which is included in the regulator circuit 
(cf. the broken line in Fig. 3). 


= a7) 
From this we have 
R (18) 
ed 
e* 
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It can be demonstrated that for such regulator circuit its equation for small increments of the variables is 
written as (for ugen = Ureg = U): 


ARe (p) = Gyp) bu — Gyhuin+ Ge duc +6,,Au,, + 6,Ai,, (20) 


R, 
Au, = Aix. R. = — Au 4 
in Sh r "she she 


Cw i + 


Taking into account the fact that T; = T, [cf. (11)] and 


1+ 
1+ pT:~ 


=i, 


(22) 


Eq. (20) can be written as 


= Gy(Au— Aujn,) + Gy Aug +G,,Au,, + G,Ai,, 


i 
(24) 


It should be noted that such a method of determining the frequency responses when the winding WT is present 
in the regulator is less accurate than the method described above, since Eq. (22) is approximate. 


For excitation of oscillations in the system from the generator side,an ac source of low frequency is applied to 
the generator excitation winding. A rotary amplifier of the type ADU-1300 was used as such a current source at 
was fed from a low-frequency generator, 

and the second was fed from a storage battery. At frequencies 
above 40-50 cps an analagous generator was used for supply- 
ing the generator excitation winding; in this case low-fre- 


Cw quency modulation of the auxiliary generator voltage achieved 
f= 2Scpp 6° by connecting the equalizer winding of the generator regulate 
20} | 260 to a low-frequency signal generator through an intermediate 
electronic amplifier. 
4 180 In all remaining respects the method of obtaining the 
6 t regulator frequency responses does not differ from that des- 
cribed above. 
4) 120 
ol Rj In order to measure the low-frequency voltages and 
. x their phases (relative to a certain reference voltage) we used 
. a? a phase detecting voltmeter of the type VP-253; the low- 
’ frequency master generator consisted of an infra-low~-fre- 
a 7 | quency signal generator of the type OS-103 (a “solartron”). 
bint Figure 5 shows the frequency responses G'w,(jw) deter- 
Fig. 5. Logarithmic amplitude- and phase~fre- mined in this fashion for an RUG-82 regulator when the 
quency responses for an RUG-82 carbon pile equalizer winding is connected and no stabilizing transformer 
voltage regulator for an open equalizer winding is used. 


G'wa- 


For this experiment the effective value for the voltage 
was Atreg = 50 mv, and Reog = 11.8 ohm. 


From Fig. 5 it is evident that in the frequency range from 10 to 100 cps the amplitude G'y, decreases with an 
increase in frequency, and the phase 9 decreases monotonically. The initial segement of the amplitude response is 
approximated by two straight lines with a junction frequency fj = 25 


c 
is 


(21) 
(23) 
where 
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At a frequency of approximately 120 cps the amplitude response has a clearly defined maximum. The shape 
of the frequency responses shown in the figure is in full agreement with expression (14). 
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The paper describes methods for designing functional potentiometers using unprofiled cases. A 
graphical method for computing functional potentiometers is given. 


Functional potentiometers are widely used in instrumentation, and in automatic and computer devices. At 
the present time profiled frames are often used to obtain the required functional dependence of the output voltage 
Upyt On the position of the movable potentiometer arm. Below we describe methods for designing functional potentio- 
meters with unprofiled frames and a resistance that is uniformly distributed with respect to length. Such frames are 
more convenient than the profiled ones to manufacture; the winding of the resistance on them and the control of the 
winding quality can be achieved more simply. 


The required functional relationship can be achieved by giving the movable arm a special shape. We shall 
call this method of obtaining the required functional relationship the shaped-slider method. The principle of opera- 
tion for such a potentiometer shall be clarified using the example of a linear displacement converter (Fig. la). The 


number 1 on the diagram indicates the resistance wound on an unprofiled frame in this particular case it may be a 
theochord in the form of a wire. 


Assume that the slider 2 slides along the rheochord made of wire shaped in the form of a flat curve. If the 
equation for the curve which describes the shape of the slider is y = f(x) in rectangular coordinates and the rheochord 


is shaped in the form of a straight line segment parallel to the y axis, then when the slider is shifted by the amount 
x is the output voltage will be expressed by the formula 


i 
Yout= "in T = "in (1) 
net 
where ujp is the rheochord supply voltage applied to terminals 3 and 4; 2 is the length of the rheochord. 
ge 


The positive direction of displacement for the slider is in this case opposite to the positive direction of the x 


axis. In obtaining formula (1) it is assumed (as it is throughout) that the potentiometer was unloaded, Thus, for ob- 
taining the dependence of u,,, on x in the form 


where c; = const, the equation for the shape of the arm must have a similar form. 


The functional relationship between u,,,, and x can be obtained in the same way if the rheochord is shaped in 
the form of a flat curve and the slider is shaped in the form of a straight line segment. We shall call this method the 
shaped rheochord method. It is not difficult to see that in this case 


out= Min (@), (2) 


where s(x) is the length of the rheochord arc enclosed between the point 0' and the slider; c, = const. 


Fig. 1. 


The special feature of these potentiometers is the fact that: a) for a rheochord which is not closed and two 
supply terminals it is possible to obtain only a monotonically varying quantity up, as a function of x, since s(x) for 
a monotonic variation of x can either only increase or only decrease; b) for obtaining the same functional relation- 
ship between up, and x it is possible to use various shapes of the rheochord curve. For example, the rheochord shape 
can be replaced at its end (cf. Fig. 1b) by the broken-line curve which is the mirror image of the segment 0° 0™ that 
is parallel to the x axis. However, all the rheochord curve shapes for which the same function uo, (x) is obtained 
must satisfy Eq. (2). 


It is possible to design a functional potentiometer using a combination of the two methods described above. 
We shall call this the combined method. 


All of these methods can be extended to apply to potentiometers with angular displacement of the slider. 
For example, in order to obtain the required dependence of the output voltage on the angle of rotation a of the 
slider when the first of the methods above is used, the slider must be given the shape of a curve whose equation in 
polar coordinates is similar to the output function required under the condition that the rheochord must be a straight 
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line segment passing through the pole of the polar coordinates (as does the axis of rotation for the slider). For a sine 
potentiometer, for example, the equation for the circle » = 2R sin « must be obtained, where p is the radius-vector, 


a is the polar angle, R is the radius of the slider circle. Schematically a sine potentiometer with a shaped slider is 
shown in Fig. 1c. 


When a shaped rheochord is used in potentiometers which convert angular displacement the same function of 
the output voltage can also be obtained for different shapes of the rheochord curve. For example, in order to obtain 


a linear dependence of upyz on x it is possible to use the following curve shapes (Fig. 1d): the semicircle a; the semi- 
circle b; the quarter of a circle c, and the curve fe d. 


The equations for the shapes of the rheochord curves that produce the same output function uoy,; must satisfy 
the equation 


@) 
in 


where cs = const. For example, in order to obtain the function 


Hout = %in sin 


the rheochord, in accordance with Eq. (3), must be given the shape of a cardidoid 


The arrangement of such a potentiometer is shown in Fig. 1c. Here, just as in the sine potentiometer cited 
above, it is necessary to isolate the point 0 of the rheochord from the slider for all a except a= 0. In both cases 


this is achieved by cutting a notch in the slider. In Fig. 1c this notch is denoted by the letter c. The notch permits 
us to avoid breaking the circuit when the slider passes the value a = 0. In order to achieve this in potentiometers 
with a shaped rheochord it is sometime sufficient to replace the specified rheochord curve shape with another curve 
which does not pass through the axis of rotation of the slider (cf. for example, the broken line in Fig. 1c). 


Other conditions being equal, a shaped rheochord has a higher resistance than the rheochord of a potentiometer 
with a shaped slider; in some cases this is a positive factor, The shortcoming of a shaped rheochord consists in the 


possibility of dissipating the uniformity of the winding at points with a large curvature. 


A common shortcoming of potentiometers with a shaped slider and a shaped rheochord is the fact that the angle 
at which the slider intersects the rheochord varies depending on the position of the slider. This latter fact can lead 
in practice to a situation where the number of turns of the winding resistance shorted by the slider does not remain 
constant; when the rheochord is in the form of a wire this can lead to a situation where the length of the rheochord 
section occupied by the slider contact is not constant, and acertain error is produced. Therefore it is desirable that 
the angle of intersection of the rheochord with the slider be as close as possible to a right angle. This can sometimes 
be achieved using the combined method. 


There is a certain difficulty in finding functions which satisfy Eq, (3). It is possible to find the curve for a shaped 
theochord using the graphical method of computation. For this purpose n points are marked on the developed length 
of the rheochord (segment 1 in Fig. a); the distances of these points from one of the rheochord ends are,respectively, 
equal to 


@) 


where ¢ (a) are the values of the unknown function which are proportional to Uy; for a = a; k= 1, 2,3,.... 
Then the rays Ry, Re, ..., Ry, ..., R are drawn from the poles; for the combined method the curves Ky, Kg, ..., Kk, «+ 
» Kp, Corresponding to the shape of the slider in the positions ... , , are drawn (cf. Fig. 2). 


By reproducing “knees” in the rheochord at the division points it is possible to achieve such a deformation of the 
theochord that the above-mentioned points lie on the corresponding rays or on the corresponding curves Ry, when the 
combined method is used. The resulting broken line will be the approximate shape of the sought for curve. 


Figure 2 shows an example of designing a sine potentiometer with the limits of rotation for the slider ranging 
from a = 0 to a = 9/2. The curve a corresponds to a potentiometer rheochord designed according to the combined 
method, and curve b corresponds to a potentiometer with a shaped rheochord. For such a graphical computation at 
the division points the output function will acquire the specified value exactly, and in the interval between those 
points it will be subject to an error caused by the process of replacing the chord arc. This error will depend on the 
number of division points and on the specific shape of the curve. 


In conclusion it should be noted that for a wire rheochord it is theoretically possible to use the methods cited 
in this paper to obtain an absolutely exact required relationship between the output voltage and the position of the 
slider. 
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REVIEWS 


REVIEW OF THE BOOK BY A. A, SHEVYAKOV 
TITLED "'THE AUTOMATION OF AVIATION POWER 
INSTALLATIONS" (DEFENSE PRESS, 1960) 


B. N. Petrov, B. S. Voronkov, and 


N. P. Kolpakova 


Translated from Avtomatika i Telemekhanika, Vol. 22, No. 4, 
pp. 556-557, April,1961 


The book being reviewed was written in accordance with a program for a course in aviation power installations 
and was accepted as a text for aviation colleges in the "Control Systems” specialization. The basic material for this 


text was drawn from a lecture course given by the author for a number of years in the Ordzhonikidze Moscow Aviation 
Institute. 


The author quite justifiably notes the present gap in university courses between the approach to engine design 
and to problems of developing automatic control systems for engines, Because of this the book makes an attempt to 


bring these two problems closer together and to study the basic characteristics of power installations from a unified 
point of view. 


The paper studies systems for controlling aircraft pisto. _ngines and almost all types of gas turbine engines — 
single-shaft and two-shaft turbojet engines, turbopropeller jet engines, two-loop turbojet engines, single~shaft and 
two-shaft turbopropellar engines. In terms of method the material in the book is arranged in such a way that in the 
beginning a study is made of control objects when the basic information on their performance characteristics is cited; 
the equations of motion are derived (in a linear approximation), and the expressions for the coefficients in the equa- 
tions of motion are derived in terms of the basic parameters which determine the operating mode of the engine. 


Based on this conclusions are drawn concerning the properties of the engine as a controlled object under all conditions 
of aircraft flight. 


Then a study is made of the possible laws for controlling the engines; the basic requirements which must be 
imposed on the control systems are formulated, and the problems involving the dynamics of these systems are studied. 


For the majority of the engine types examples are given for performing computations with determine the values 
of the coefficients in the linearized equations of motion and analyze the dynamics of the control system. The prob- 


lems involving the dynamics of control systems are studied using specific examples of control networks that are used 
in practice. 


The existing educational plans in aviation colleges in the “Control Systems” and “Instrumentation” specialities 
are laid out in such a way that almost no attention is devoted to problems of engine theory. Moreover, in order to 
understand the essence of this subject a minimum preliminary knowledge in this field is necessary. In this regard the 
book under review makes an attempt to approach the subject in a broader fashion; the book thus quite justifiably in- 
cludes sections devoted to basic information on the performance characteristics of engines. 


An important part of the book is devoted to problems involving the dynamics of control systems for turbojet 
engines with a forcing chamber, since this type of engine is one of the most common used at present. Here a number 
of possible methods for controlling the forcing loop are studied: regulating the pressure before and after the com- 
pressor, regulating the degree of gas expansion onto the turbine, regulating the gas temperature. 


For each of the control methods the dynamics problem is solved in such a way that the evaluation of the system 
is performed according to the transient responses caused by the most substantial perturbations. The same section clari- 
fies problems associated with the process by which the engine enters into a forcing mode. Here studies are also made 
of several possible methods for controlling the engine when the forcing loop controller acts either on the fuel ex~- 
Penditure in the forcing chamber or on the valve of reactive jet. Under these conditions we take into account different 


values for the time shift relative to the beginning of controller operation and the instant at which the forcing mode 
begins. The transient responses treated in the book permit a proper evaluation of the possibilities of each system. 


A significant place in the book is reserved for problems involving the dynamics of control systems for turbo- 
propeller engines; these problems also apply to the basic types of gas turbine engines which are widely used at present, 
These problems are studied for single-shaft motors with one and two propellers (with a differential reducer), or for 
two-shaft motors with one propeller. Here a study is made of those control laws for which the number of propeller 
revolutions is maintained at the specified number by means of controllers which affect the angles at which the propel- 
ler blades are set; the specified gas temperature is maintained using controllers which operate in closed and open 
loops and affect the fuel discharge. The problems associated with the effect of the stabilizing devices in these con- 
troliers are clarified along with the effect of the engine torque limiter on the properties of the control system. Data 
on the transient responses makes it possible correctly to evaluate the advantages and disadvantages of each system 
which is studied. 


An important section of the book is the chapter dealing with problems of autonomous control of gas turbine 
engines. With respect to single~shaft turbopropeller jet engines and turbopropeller engines with one propeller the 
book demonstrates the possibilities of such control systems and compares them with nonautonomous systems. Here 
a study is made of conditions governing complete and partial autonomy, and the advantage of using autonomous 
control systems is demonstrated. 


One of the very important problems which to a great extent determines the basic properties of any aviation 
power installation is the acceleration characteristic for the engine. The narrow range of operating stability for the 
compresser and the combustion chamber, as well as the regorous limits with respect to the maximum allowable gas 
temperature at the turbine input, impose especially difficult conditions on the correct selection of parameters for 
the acceleration system in the engine. A separate chapter is devoted to this subject, where a study is made of the 
behavior of an engine acceleration system using a phase plane (since a nonlinear problem is being solved). The 
method given for computing the acceleration systems of turbojet engines also makes it possible to select the con- 
troller parameters in such a way that the necessary margin is achieved with respect to the maximum allowable values 
of gas temperature and the pumping characteristics of the compressor. 


- The book compares favorably with others in terms of its completeness, the use of specific examples, and the 
way in which the problems dealt with coincide with the curriculum. 


Besides systematizing previous known material, the book offers a number of new results, including an analysis 
of the dynamics of two-shaft and two-loop turbojet engines, single-shaft turbopropeller jet engines, turbopropeller 
engines with a differential reducer and coaxial variable-pitch propellers, autonomous systems for turbojet engines 
and turbopropeller engines, and acceleration systems for turbojet engines. Investigations of the transient responses 
are made for all of these case, and the corresponding estimates are given. 


The book includes certain clearly identified material which extends beyond the range of the student curriculum; 
this is done to give the book a scope beyond that of a textbook. 


The book is adequately illustrated and is written in easily understandable language. 


The book presents extensive material on the classification of aviation engines as controlled objects. However 
(and this is one of its most substantial shortcomings) the book does not provide any comparison and classification of 
the automatic control systems for the indicated engines. The author has not made a critical analysis of various types 
of control systems. Critical remarks should also be made on the way engineering design problems are treated; this 
applied particularly to problems of designing the automation equipment. The book essentially contains no descrip- 
tion of the design and technological solutions, nor does it include a survey of typical engineering specifications and 
conditions for aviation engine automation systems; problems of improving the control reliability are also not dealt 
with. 


The “auxiliary” aggregates of the control systems—the fuel supply, the starting system for the engine, the safety 
features, etc.—are treated in a very cursory way. New problems in the theory of aviation engine control are not treat- 
ed at all (problems of optimization and a consideration of discontinuous perturbations). 


As a whole this book is an original text which deals in easily accessible form with the problems set forth in it; 
its publication is very timely. 


